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PREFACE 


1 


In  the  aftermath  of  the  1939—45  war,  as  a  more  scientific  approach  began  to  supplement  cut-and*try  as  a  basis  for  the 
design  of  parachutes,  W.D.Brown  first  established  a  sound  aerodynamic  foundation  for  his  subject  in  his  book,  ‘Parachutes'. 
Roles  which  were  being  fulfilled  by  the  parachutes  which  he  described  in  it  were  firstly  man-carrying  for  both  life-saving  and 
military  applications;  secondly  weapon  and  store  dropping;  and  thirdly  deceleration  of  aircraft.  Looking  back  over  the  35 
years  that  nave  elapsed  since  this  publication,  although  their  form  has  been  greatly  expanded,  all  of  these  are  still  essential 
tasks  for  parachutes  to  fulfil.  However,  additional  ones  also  are  now  required.  For  example,  crews  of  spacecraft  as  well  as 
aircraft  make  encapsulated  descents  using  parachutes  as  essential  parts  ot  their  total  escape  systems  and  in  addition  to  drag 
production  contemporary  parachutes  may  also  need  to  demonstrate  sifeificant  lift-producing  abilities  or  the  capability  of 
rotation  at  chosen  spin  rates  during  their  descent.  Not  ony  have  parachufps  been  used  on  a  number  of  occasions  to  assist  re¬ 
entry  into  the  earth's  atmosphere,  they  have  also  been  employed  in  landing  essential  instruments  on  to  other  planets  as  well 
as  on  to  the  planet  earth. 

They  have  also  had  to  find  their  place  within  the  much  wider  classification  of  aerodynamic  decelcrators,  in  which  they 
share  a  primarily  decelerative  task  with  bailoons  and  various  metallic  non-inflatable  devices.  In  his  glossary  in  1 95 1 ,  Brown 
defined  a  parachute  as  ‘an  umbrella-shaped  device  to  produce  drag,  commonly  used  to  reduce  the  rate  of  descent  of  a  falling 
body'.  A  more  precise  and  limiting  definition  is  now  customary.  Within  this  broader  classification  of  aerodynamic 
deceierators  parachutes  are  now  considered  to  be  that  class  of  drag-producing  bodies  whose  essential  characteristics  include 
their  flexibility,  their  inflated  shape  being  dependent  on  the  flow  field  which  surrounds  them. 


A  greater  precision  in  definition  epitomises  changes  taking  place  in  many  engineering  fields,  not  only  in  that  of 
parachute  aerodynamics.  With  the  availability  in  most  design  offices  of  powerful  mainframe  computers  and  in  many 
situations  of  relatively  inexpensive  mini-  and  micro-computers,  new  possibilities  exist  for  engineers  to  establish  the  relevant 
basic  relationships  and  to  develop  their  subsequent  solutions.  For  these  task  to  be  performed  adequately,  basic  principles 
must  be  appreciated  and  agreed  sign  conventions  implemented. 

It  is  to  meet  all  these  kinds  of  need  that  ‘  The  Aerodynamics  of  Parachutes'  has  been  written.  In  the  subject  of  parachute 
aerodynamics  this  AGARDograph  is  envisaged  as  a  direct  descendant  of  Brown's  'Parachutes',  for  it  has  the  same  emphasis, 
that  of  selecting  the  principal  aerodynamic  characteristics  of  parachutes  and  the  various  known  factors  which  affect  these 
characteristics'.  It  takes  into  account  not  only  many  of  the  subsequent  publications  which  have  been  summarised  in  the  1 963 
and  1978  United  States  Air  Force  Parachute  Design  Guides,  but  also  the  proceedings  of  the  American  Institute  of 
Aeronautics  and  Astronautics  Aerodynamic  Decelerator  Conferences  which  have  been  held  every  two  and  a  half  years,  the 
Helmut  G-Heinrich  Decelerator  Systems  Engineering  Short  Courses  which  took  place  in  1 983, 1 985  and  1 987  and  ‘  The 
Parachute  Recovery  System  Design  Manual',  which  will  shortly  be  used  by  the  United  States  Naval  Weapons  Center. 

It  has  been  anticipated  that  its  main  readers  will  be  recent  engineering  graduates  entering  research  establishments, 
parachute  companies  or  related  industries.  In  its  preparation  some  appreciation  on  the  part  of  the  reader  of  basic  mechanics, 
elementary  fluid  mechanics  and  the  principles  of  computing  has  been  assumed. 

Apart  from  Karl-Friedrich  Doherr  and  my  own  research  associates,  too  many  other  individuals  have  contributed  their 
components,  criticisms  and  suggestions  for  me  to  mention  their  names  individually.  I  can  only  hope  that  they  will  recognise 
their  invaluable  contributions  in  the  publication  which  has  resulted  from  all  of  our  efforts. 


David  Cockrell 
Leicester  —  1987 


Au  lendemain  de  la  guerre  de  1 939—45,  alors  qu’un  demarche  plus  scientifique  commcnfait  a  completer  les  methodes 
empiriques  du  genre  "on  decoupe  et  on  essaie”  comme  base  de  la  conception  des  parachutes,  W.D.Brown  fut  le  premier  a 
elaborer  une  theorie  aerodynamique  saine  pour  ce  qui  faisait  l’objet  de  son  livre  “Les  parachutes'.  Les  fonctions  remplies  par 
les  rnateriels  qu’il  y  decrivait  etaient  d’abord  I’emport  des  hommes  dans  le  double  but  de  la  sauvegarde  de  la  vie  humainc  ct 
des  applications  militaires,  en  deuxieme  lieu  le  largage  d’armes  et  d'approvisionnements,  cn  troisieme  et  dernier  lieu  la 
deceleration  des  avions  it  l’atterrissage.  Si  on  etudie  les  35  annecs  qui  se  sont  ecoulees  depuis  cctte  publication,  tous  ccs 
cmplois  sont  encore  essentiellement  ceux  que  Ton  attribue  aux  parachutes  meme  si  teur  aspect  exterieur  s'est  bcaucoup 
diversifie.  Neanmoins  il  en  faut  maintenant  quelques  model  es  supplemental  res.  Par  exemple,  les  equipages  des  vehicules  de 
I'espace  comme  ceux  des  avions  font  des  descentes  enferme  dans  des  capsules  equipees  de  parachutes  qui  constituent  la 
partie  essenticlle  de  l'ensetnble  de  leur  systeme  d’evacuation;  et  en  plus  de  la  production  d'engins  bases  sur  la  trainee,  il  peu; 
aussi  etre  demande  aux  parachutes  modem  es  de  posseder  des  qualites  de  portance  ou  d'aptitude  a  toumer  sur  eux-memes  a 
une  vitesse  donnee  au  cours  de  leur  descente.  Des  parachutes  ont  etc  utilises  non  settlement  en  de  nombreuses  occasions 
pour  fad  liter  la  rentrie  dans  I’atmosphere  terrestre,  mais  aussi  comme  instruments  essentiels  d’atterrisage,  quit  s’aguse  de  sc 
poser  sur  not  re  globe  ou  sur  cfautre  planetcs. 

II  a  egalcment  fallu  leur  trouver  un  crcneau  dans  la  dasse  bcaucoup  plus  vaste  de  decelerateurs  aerodynamiques  ou  its 
rempUssent,  concurrement  avec  les  ballons  et  divers  dispositifs  metaliiques  non  gonflabks,  une  tache  qui  consistc 
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cssentiellemcnt  a  assurer  un  freinage.  Dans  son  glossaire  de  195 1 ,  Brown  a  defini  le  parachute  comme  “un  dispositif  en 
forme  de  paraptuie,  foumissant  une  trainee,  utilise  communement  pour  diminuer  la  vites.se  de  desccnte  d'un  corps  qui 
toinbe".  La  definition  courante  actuelle  est  plus  precise  et  plus  limitative.  Au  sein  de  la  categorie  tres  vaste  des  dcceierateurs 
aerodynamiqucs,  les  parachutes  sont  generatcurs  de  trainee  dont  une  des  caracteristiques  esscntiellcs  est  la  souplssse 
d’emploi,  car  une  fois  gonfles,  leur  forme  depend  de  l’ecoulernent  de  Pair  autour  d'eux. 

Cette  meilleure  precision  dans  la  definition  resume  bien  revolution  qui  s’est  faite  dans  beaucoup  de  secteurs  de 
I'industrie,  et  pas  seulement  dans  Paerodynamique  des  parachutes.  Grace  k  la  presence  de  puissants  ordinateurs  centraux 
dans  la  plupart  des  bureaux  d'etude  et,  dans  beaucoup  de  cas,  &  l’utilisation  de  mini-  et  de  micro-ordinateurs  relativcment 
peu  couteux,  les  ingenieurs  disposent  de  possibility  nouvelles  pour  etablir  les  fonnules  de  base  voulues  et  devclcpper  les 
solutions  qui  en  decoulent  en  vue  de  la  conception.  Pour  que  ces  taches  puissent  etre  executees  efficacement,  it  faut  bien 
cemer  les  principes  de  base  et  appliquer  les  conventions  de  symbolique  convenues. 

Ccst  pour  satisfaire  tous  ces  besoins  que  1‘ouvrage  “Aerodynamique  des  parachutes"  a  etc  ecrit.  Dans  le  domaine  dont 
il  port  le  titre.  cet  "AGARDographe"  est  considere  comme  la  suite  directe  du  livre  “Les  parachutes”  de  Brown,  car  i!  insiste 
sur  le  meme  theme:  “les  principaies  caracteristiques  des  parachutes  et  les  divers  facteurs  connus  qui  les  affectent”.  II  tient 
compte,  non  seulement  du  grand  nombre  de  publications  posteri cures  qui  ont  ete  repertoriees  et  resumees  dans  les  “Guides 
de  conception  des  parachutes”  de  l’armee  de  Pair  americaine  de  1 963  a  1 978,  mais  cgalement  des  comptes  rendus  de  debats 
dy  conferences  de  PInstitut  americain  d’Acronautique  et  d'Astronautique  (“American  Institute  of  Aeronautiy  and 
Astronautiy")  sur  la  deceleration  aerodynamique  qui  se  sont  tenuy  tous  ly  deux  ans  et  demi,  dy  cours  techniquy  abreges 
de  Helmut  Giieinrich  sur  la  technique  dy  systemy  de  deceleration  qu‘il  a  donnes  en  1 983, 1 985  et  1 987.  et  du  “Manuel  de 
conception  d’un  projet  de  recuperation  par  parachute”  qui  sera  prochainement  public  par  le  Centre  des  armcments  navals 
americain  (“U.S.  Naval  Weapons  Center"). 

On  a  prevu  que  ly  principaux  lecteurs  de  cet  ouvrage  sentient  ly  ingenieurs  fraichement  diplomes  qui  sont  sur  ly 
point  d'entrer  ly  etablissements  de  recherche,  ly  society  de  fabrication  de  parachuty  ou  ly  industriy  qui  leur  sont 
associey.  Pour  preparer  sa  presentation,  on  a  suppose  que  le  lecteur  possedait  quelqucs  connaissancy  de  mecanique 
fondamentale,  dy  notions  elementairy  de  mecanique  dy  fiuidy  et  des  principy  de  Pinformatique. 

En  dehors  de  Karl-Friedrich  Doherr  et  de  my  propry  associes  en  matiere  de  recherche,  ly  autry  personny  qui  ont 
apportc  la  contribution  de  leurs  connaissancy,  de  leurs  critiquy  et  de  leurs  suggy  tions  sont  trops  nombreusy  pour  que  je 
puisse  |y  citer  touty  individuellement.  Pypere  seulement  qu'ellcs  pourront  reconnaltrc  au  passage  ly  apports  inytimable,; 
qu’clly  ont  fails  a  cette  publication  qui  yt  Paboutissement  de  tous  nos  efforts  conjugues. 

David  Cockrell 
Leicy  ter  — 1987 
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1. INTRODUCTION 


l.t  FUNCTION  OF  A  PARACHUTE 

The  word  'parachute'  is  derived  from  the  French  words  parare ,  meaning  to  shield,  or  to  ward  off  anti  chute , 
meaning  a  fall.  Thus  the  word  parachute  means  any  natural  or  artificial  contrivance  which  serves  to  check  a  fall  through 
the  air  and  thereby  support  in  the  air  sorae  load  or  store.  But  for  present  purposes  this  is  too  wide  a  definition.  Within 
thli  AGARDograph  the  parachute  will  be  considered  to  be  one  subdivision  within  a  much  larger  classification  of 
aerodynamic  decelerators.  Ibrahim14  defines  the  latter  as  ‘devices  whose  primary  function  is  to  maximise  the  drag  of 
systems  of  which  they  form  a  part’  and  parachutes  as  ‘flexible,elastic  bodies  whose  inflated  shapes  are  dependent  on  the 
flow  conditions*.  Because  of  the  structural  form  of  the  rigid  decelerator,  the  improper  primary  function  of  the  flexible 
wing  and  the  invariability  of  the  ballute’s  inflated  shape  (a  balhue  is  a  cross  between  a  balloon  and  a  parachute,  used  as 
a  high-speed  aerodynamic  decelerator)  these  devices  are  excluded  from  consideration  as  parachutes  here,  except  as 
possible  illustrations  o',  certain  aerodynamic  principles. 

Although  the  primary  function  of  the  parachute  is  the  maximisation  of  system  drag,  among  possible  secondary 
rbles  may  be  the  provision  of  a  horizontal  velocity  component  (termed  drive )  to  the  system  or  even  to  be  the  means 
whereby  this  system  acquires  a  lift  force.  Thus  the  distinction  between  parachutes  and  flexible  wings  may  not  be  as 
clear  cut  as  it  first  appears. 

From  their  initial  man-carrying  role,  parachutes  have  been  required  to  undertake  a  much  wider  range  of  tasks.  As 
well  as  men,  stores  are  often  decelerated  by  parachutes.  These  include  weaponry  as  well  as  other  supplies.  Parachutes 
can  also  be  used  to  decelerate  rapidly-moving  vehicles,  such  as  aircraft  and  motor-cars.  Sometimes,  instead  of  an 
individual  ejected  from  a  malfunctioning  aircraft,  an  entire  crew-carrying  module  is  decelerated  by  a  parachute.  As  Fig. 
1.1  indicates,  parachutes  must  be  capable  of  operation  over  a  wide  range  of  velocities  and  hence  dynamic  pressures,  in 
very  different  environments.  For  example,  extra-terrestrial  applications  may  be  required  when  instruments  have  to  be 
landed  on  other  planets.  Even  on  the  earth's  surface,  as  natural  aerodynamic  decelerators  such  as  dandelion  seeds  bear 
witness,  some  applications  of  parachutes  call  for  unusual  designs. 
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PIC.  1.1  PARACHUTE  PERFORMANCE  ENVELOPE 
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1.2  AERODYNAMICS:  ONE  AMONG  MANY  DESIGN  CRITERIA 

The  aerodynamics  of  parachutes  is  only  one  aspect  of  their  multi-faceted  design.  Their  flexible  structure  must  be  of 
light  weight,  yet  strong  enough  to  withstand  the  high  loads  imposed  during  inflation.  Their  deployment  must  be  both 
simple  and  reliable.  They  must  be  relatively  easy  to  manufacture  so  that  their  cost  is  small  compared  with  that  of  the 
store  which  they  decelerate  and  then  deliver.  Parachute  aerodynamics  is  concerned  with  the  mechanism  of  inflation,  the 
determination  of  flight  mechanics,  such  as  the  nue  of  system  descent  in  a  given  application,  the  drive  and  the  lifting 
behaviour  of  the  parachute.  Included  are  system  stability  characteristics  such  as  equilibrium  angles,  frequency  of 
oscillation  and  the  oscillation  damping  rates.  To  predict  satisfactorily  all  these  characteristics  the  necessary 
mathematical  relationships  must  first  be  developed  and  relevant  data  acquired  through  appropriate  experimental 
programmes.  Having  beat  acquired,  these  data  must  then  be  made  generally  available.  The  Parachute  Design  and 
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Performance  Data  Banku,  established  at  the  United  States  Air  Force  Flight  Dynamics  Laboratory  in  the  period  1970- 
73,  was  a  significant  step  in  data  dissemination  and  this  ACARDograph  has  been  written  to  further  assist  data 
formalisation,  application  and  dissemination. 

1.3  BRIEF  HISTORICAL  SURVEY  (REFERENCES  1.3  TO  1.11) 

Although  evidence  of  parachute-like  devices  to  lower  both  animals  and  humans  from  high  towers  exists  in  Chinese 
archives  from  as  early  as  tbe  12th-century  and  in  1514  sketches  of  parachutes  were  made  by  Leonardo  da  Vinci,  the 
first  authenticated  parachute  descent  was  not  made  until  October  22 1797,  when  Andrf-Jacques  Ganerin  jumped  from  a 
balloon  over  Paris.  Early  users  of  parachutes  were  stunt  men,  descending  from  towers  or  from  tethered  balloons  for  the 
entertainment  of  spectators.  By  the  early  nineteenth  century,  exhibition  parachute  descents  from  balloons  were  being 
made  all  over  the  world  and  in  this  capacity,  in  1808  a  parachute  first  saved  a  human  life. 

At  the  outbreak  of  the  1914  world  war,  participants  on  both  sides  ’tudied  activities  behind  the  enemy  lines  by 
stationing  observers  in  baskets,  slung  beneath  tethered  hydrogen-filled  talloons.  As  these  balloons  exploded  if  they 
were  hit  by  machine-gun  fire,  on  the  approach  of  enemy  aircraft  the  observers  would  bail  out,  using  for  this  purpose 
cotton  parachutes,  some  9-11  metres  in  diameter  which  were  tethered  to  their  baskets.  A  large  number  of  observers' 
lives  were  thus  saved,  there  being  407  successful  parachute  descents  in  France  by  members  of  the  British  Ballocn 
Wings  alone  and  a  further  125  by  members  of  the  United  States  Forces.  For  French  observers  a  parachute  system  was 
evolved  in  which  the  entire  basket  was  retrieved,  becoming  the  first  recorded  'encapsulated'  retrieval  by  pararhute. 
Similar  observational  practices  were  carried  out  by  the  German  Army,  with  a  corresponding  high  success  rate  in 
observer  retrieval 

When  aircraft  entered  the  war  the  use  of  parachutes  by  aircrew  was  delayed.  The  probable,  reason  for  not  using  them 
initially  was  because  of  difficulties  in  egress  from  aircraft  cockpits.  Since  it  took  valuable  time  to  put  the  parachute 
harness  on  and  to  extract  the  parachute  from  its  container  which  was  fixed  to  the  aircraft,  often  there  was  insufficient 
altitude  remaining  for  the  parachute  to  fully  open.  But  eventually,  die  prime  importance  of  pilots'  lives  was 
recognised,  probably  because  of  the  considerable  investment  in  training  cost  and  time  which  they  represented.  The  first 
recorded  saving  of  life  from  an  aircraft  by  a  parachute  was  in  1916.  At  that  time  parachutes  were  opened  by  static  lines 
attached  to  the  aircraft,  the  opening  being  delayed  until  the  parachutist  was  well  clear  of  the  machine. 

By  this  stage  in  the  war,  individual  aviators  on  the  German  side  were  equipping  themselves  with  appropriately 
modified  Heinecke  parachutes  which  had  originally  been  intended  for  balloon  observers.  As  these  parachutes  opened 
they  lifted  the  aviators  clear  of  their  cockpits.  Following  the  development  of  'packaged'  or  pack  parachutes  by  Charles 
Broadwick  and  others  in  the  early  190D's,  all  combatants  rapidly  made  the  necessary  developments  in  materials  and  in 
parachute  packing.  On  April  28  1919,  Leslie  L.  Irvin  made  the  first  free  parachute  descent,  from  1500  ft  above  the 
ground.  By  this  lime,  parachutes  were  in  regular  use  for  the  dropping  of  flares  and  in  1918  they  were  often  the  means 
by  which  spies  wet?  infiltrated  behind  enemy  lines. 

The  first  parachute  designed  for  military  personnel  was  standardised  in  1924.  After  that  time,  first  in  the  United 
States  and  later  in  Great  Britain,  the  use  of  parachutes  became  compulsory  for  aircrew.  By  about  1930  the  Soviet 
Army  had  begun  to  equip  and  train  some  of  its  units  for  airborne  operations,  using  parachutes.  Corresponding  German 
units  were  deployed  in  Holland  and  Belgium  during  the  early  stages  of  the  1939-45  war. 

By  this  stage  in  many  countries  a  systematic  testing  and  development  programme  had  become  essential.  There  was 
an  over-riding  need  for  reliability,  thus  for  a  better  appreciation  of  parachute  materials  characteristics,  of  structural 
strengths,  opening  factors,  drag  characteristics  and  stability  behaviour.  Research  took  place  in  many  places  but 
increi  -ingly  in  the  United  Kingdom  and  in  Germany.  By  the  outbreak  of  the  second  world  war  in  1939  there  was 
considerable  experience  in  using  parachutes  for  weapon  stabilisation,  required  both  for  impact  attitude  and  the  need  to 
obviate  high  g-loading  in  the  direction  normal  to  that  of  the  weapon  axis,  in  tbe  dropping  of  supplies  by  parachute  And 
in  paratrooping.  During  that  war  there  was  considerable  development  in  ail  these  applications  as  well  as  in  the  aircraft 
deceierator  rNe,  made  necessary  through  both  the  advent  of  dive  bombing  and  tbe  rapid  deceleration  on  landing  required 
by  some  fighter  aircraft.  During  tbe  1930's  tbe  needs  for  high  aircraft  deceleration  led  to  the  development  of  ribbon 
parachutes  by  Georg  Madelung.  At  the  high  speeds  which  were  necessary  such  parachutes  were  able  to  provide  the 
required  low  opening  shock  loads  and  also  exhibit  stability  in  pitch. 

These  various  applications  were  demanding  differing  parachute  characteristics,  for  example  a  low  degree  of  parachute 
stability  tolerable  to  a  member  of  aircrew  making  an  emergency  escape  from  his  aircraft  would  be  quite  unacceptable  to 
a  regular  parachutist  such  as  a  paratrooper,  or  for  an  aircraft  deceierator  system.  For  such  ejector  systems  knowledge  of 
the  relevant  parameters  influencing  parachute  inflation  became  essential  so  that  satisfactory  predictions  of  the  time 
taken  foe  inflation  and  tbe  corresponding  forces  which  were  developed  could  be  achieved.  Using  parachutes,  guided 
missiles,  such  as  the  V.l  and  V.2,  as  well  as  missile  components  were  successfully  recovered  in  1944  and  the  earliest 
ejector  seat  deceleration  was  made  by  parachute  about  1944-6,  the  idea  for  so  doing  originating  in  Sweden. 

During  the  1939-45  war,  at  the  various  research  establishments  parachute  sections  were  established.  For  example, 
at  the  Royal  Aircraft  Establishment  under  WU.Brown,  the  British  Parachute  Section  was  established  in  1942.  After 
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this  war  was  over,  when  T-FJohns,  a  member  of  this  Section,  published  the  Report  "Parachute  Design",  he  stated  that 
most  usual  requirements  for  parachutes  were:* 

(i) .  that  they  will  invariably  inflate; 

(ii) .  that  they  will  develop  specified  drag  forces  at  particular  descent  speeds; 

(iii) .  that  they  will  be  sufficiently  strong  to  withstand  opening  at  speeds  which  are  usually  higher  than  their 

descent  speeds,  and 

(iv) .  that  they  will  give  specified  degrees  of  stability  to  the  payloads  to  which  they  are  attached. 

In  this  AGARDograph,  the  aerodynamic  aspects  of  requirement  (ii)  are  discussed  in  chapter  2,  those  of  (i)  in  chapter 

5  and  (iv)  in  chapter  4,  The  aerodynamic  aspects  of  requirement  (iii)  are  considered  in  chapters  2, 3  and  4. 

_  After  the  end  of  the  second  world  war,  in  the  United  States  the  military  engagements  in  Viet-Nam  and  elsewhere 

i  stimulated  more  parachute  research  into  gliding  parachutes  such  as  the  ram-air  inflated  textile  wing  originally  proposed 

i  in  1961  by  Jalbert,  the  emergency  escape  of  aircrew  and  the  airborne  delivery  of  personnel,  stores  and  weapons  as  well 

*  as  into  aircraft  retardation  and  vehicle  recovery  over  a  wide  dynamic  pressure  range.  Deceleration  through  the 

deployment  of  a  series  of  parachute  canopies  in  a  number  of  separate  stages  became  commonplace.  By  the  1960's  the 
ribbon  parachutes  developed  for  this  purpose  were  used  for  the  deceleration  of  the  United  States  astronauts  returning 
from  die  Moon  in  the  Mercury,  Gemini  and  Apollo  spacecraft  as  well  as  those  in  the  Soviet  Union's  Vostok  and 
Soyuz  space  vehicles.  Yuri  Gagarin  safely  landed  Vostok  I  by  parachute  in  April  1961  and  in  February  1962  John 
Glenn  used  a  ringsail  ribbon  parachute  to  land  a  Mercury  spacecraft  In  July  1976,  using  parachutes,  successful 
landings  of  the  first  of  two  Viking  spacecraft  was  made  on  the  planet  Mars. 
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2.STEADY-STATE  AERODYNAMICS 

Although  it  is  more  logical  to  begin  by  describing  the  deployment  and  the  inflation  of  the  parachute  canopy, 
proceed  with  its  deceleration  and  then  to  consider  its  behaviour  in  the  steady  state,  basic  aerodynamic  concepts  for 
parachutes  are  most  readily  introduced  by  first  considering  their  fully-deployed  steady  descent. 

2.1  SOME  DEFINITIONS.  RELEVANT  DIMENSIONLESS  PARAMETERS 

In  flight  mechanics  the  characteristic  forward  direction  of  vehicles  in  motion  is  first  determined.  Then  a  common 
procedure  is  to  establish  orthogonal  sets  of  axes  which  are  fixed  in  these  vehicles,  passing  through  an  origin  which 
itself  is  fixed  in  the  vehicle.  Such  axes  are  referred  to  as  body  axes.  Customarily  the  axis  O-x  is  so  positioned  that  it 
points  in  that  characteristic  forward  direction.  As  Fig.  2.1  illustrates,  the  angle  of  attach,  a  4s  the  angle  measured 
between  the  component  in  the  Oxz  plane  of  the  resultant  airflow  V»  in  that  forward  direction  and  the  body  axis  O-x. 
The  body  axes  are  right-handed  in  direction. 

As  Terns  and  Symbols  for  Flight  Dynamics11  makes  clear,  in  aeronautical  parlance  the  term  angle  of  incidence  is 
no  longer  an  acceptable  alternative  for  angle  of  attack,  a . 

Like  any  other  immersed  body,  when  a  parachute  moves  through  a  fluid  a  resultant  aerodynamic  force  R  is 
developed  on  it.  Experiments  can  be  devised  by  which  to  measure  both  the  magnitude  of  this  force  and  its  moment 
about  any  specified  location.  Though  these  measurements  are  sufficient  to  define  the  line  of  action  of  this  resultant 
force  they  do  not  determine  the  precise  position  of  the  centre  of  pressure,  a  specific  point  on  that  line  of  action  at 
which  the  resultant  aerodynamic  force  can  be  considered  to  act.  In  most  aerodynamic  applications  the  location  of  the 
centre  of  pressure  is  determined  by  convention.  Thus,  for  a  section  of  aircraft  wing  section  or  for  a  gliding  parachute 
the  centre  of  pressure  position  is  defined  as  being  at  the  intersection  of  the  resultant  aerodynamic  force  line  of  action 
with  the  chord  line  of  the  aerofoil  section  which  constitutes  the  wing  or  the  gliding  parachute.  A  gliding  parachute, 
such  as  that  with  a  ram-air  canopy  described  in  Appendix  2B,  is  one  which  is  capable  of  imparting  a  horizontal 
component  of  velocity  or  drive  to  the  parachute  and  its  payload.  Momentarily  it  is  possible  for  a  system  comprising 
a  gliding  parachute  and  payload  to  develop  a  resultant  lift  force.  In  contrast  is  the  conventional  parachute,  possessing 
solely  a  drag-generating  rNe.  A  number  of  conventional  parachute  canopies  are  illustrated  in  Appendix  2A.  In  some 
situations  this  distinction  between  these  two  types  of  parachute  canopies  becomes  artificial,  since  conventional 
parachutes  become  gliding  parachutes  if  appropriate  panels  are  removed  from  the  canopy.  When  this  occurs  either 
definition  of  parachute  could  be  adopted,  whichever  is  the  more  convenient 

In  the  physical  appreciation  of  parachute  behaviour,  such  as  when  formulating  and  solving  equations  of  motion,  it 
is  sometimes  desirable  (though  not  essential)  to  know  the  centre  of  pressure  location.  For  conventional  parachutes 
the  centre  of  pressure  position  is  defined  to  be  at  the  the  intersection  of  the  line  of  action  of  the  resultant  aerodynamic 
force  with  the  parachute  axis  of  symmetry. 
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D*  =  NS./jt)]1"  (24) 
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Since  the  nominal  total  surface  area  of  the  canopy  is  not  always  clearly  defined  there  can  be  confusion  over  the 
magnitude  of  the  parachute’s  nominal  diameter  D,.  To  avoid  this  confusion  the  term  constructed  diameter,  D,  is 
sometimes  used.  The  Recovery  Systems  Design  Guide,  reference  1.9,  defines  the  constructed  diameter  of  a  parachute 
canopy  as  ‘the  distance  measured  along  the  radial  seam  between  points  where  the  maximum  width  of  opposing  gores 
intersects  that  radial  seam’. 

The  Mach  number  Ma  equals  V^a ,  where  a  is  the  local  speed  of  sound  in  the  undisturbed  fluid.  The  ways  in 
which  these  and  other  aerodynamic  force  and  moment  coefficients  vary  with  the  angle  of  attack,  Reynolds  number  and 
Mach  number  are  described  in  Section  22. 

For  gliding  parachutes  the  aerodynamic  reaction  R  is  usually  expressed  in  terms  of  its  two  components,  lift  L  and 
drag  D  respectively  perpendicular  to  and  parallel  to  the  resultant  airflow,  as  illustrated  in  fig.2.2  The  positive 
direction  of  lift  is  in  the  opposite  sense  to  the  weight  of  the  system  and  the  positive  direction  of  drag  is  in  the 
opposite  sense  to  the  parachute’s  resultant  velocity.  Like  tangential  and  normal  force  components,  lift  and  drag  forces 
are  similarly  expressed  in  terms  of  non-dimensional  coefficients  as  Ct.aiKi  C,,: 


cl  = 

and  Co  = 


L 


Ws° 

(2.5) 

D 

2PVRS» 

(2.6) 

The  non-dimensional  force  coefficients  Q.  and  Q,  can  be  expressed  as  functions  of  the  angle  of  attack,  the 
Reynolds  number  and  the  Mach  number  in  exactly  the  same  way  as  were  the  force  coefficients  C,.  and  CN  in 
equation  23. 

To  avoid  any  confusion  in  determining  the  sense  of  the  lift  component  force  it  is  preferable  to  confine  the  use  of 
tangential  and  normal  force  components  to  conventional  parachutes  while  reserving  lift  and  drag  components  for 
gliding  parachutes.  When  considering  flight  mechanics  however,  the  component  of  aerodynamic  force  in  the  direction 
of  the  relative  airflow,  i.e.  the  drag  component  is  often  required,  so  transformation  is  necessary  from  one  set  of 
aerodynamic  force  components  to  the  other.  The  relationships  shown  in  Fig.2.3  are: 


Conventional 
Parachi 


Fig.  2.3  RELATIONSHIPS  BETWEEN 
THE  AERODYNAMIC  COEFFICIENTS 
CL  AND  C„  AND  THE  COEFFICIENTS 
CN  AND  CT 


CB  = 

CT  coa  a 

♦  C„  sin  a 

(2.7) 

anti 

Ci  * 

C„  cos  a 

•  CT  tin  a 

(2.8) 

Correspondingly, 

CT  « 

CB  coa  a 

•  CL  sin  a 

(2.9) 

▲ 


Vf 
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and  C*  «  Cfc  cot  o  ♦  CB  tin  o  (2.10) 


And  since  l/pV},  a  *  (2.11) 

where  the  subscripts  o  and  E  respectively  imply  measurement  at  sea  level  and  the  equivalent  air  speed,  the 
expressions  for  the  aerodynamic  force  components  in  equations  2,1, 2.2, 25  and  2.6  can  be  written  in  terms  of  the 
equivalent  air  speed  V,.  The  latter  is  the  appropriately-corrected  speed  be  recorded  by  an  air  speed  indicator. 


2.2  SOME  STEADY-STATE  AERODYNAMIC  CHARACTERISTICS 

To  illustrate  the  functional  relationships  expressed  in  equation  23,  some  typical  aerodynamic  characteristics  of 
various  canopy  shapes  are  now  considered  A  list  of  the  most  common  parachute  canopy  shapes,  together  with  a  brief 
account  of  their  aerodynamic  characteristics,  is  appended  to  Section  2 
2.2.1  Shape  of  Parachute  Canopy 

As  explained  in  the  Introduction,  the  process  of  parachute  design  is  inevitably  one  of  making  compromise 
decisions.  The  shape  of  the  parachute  canopy  is  determined  by  considering  all  the  idles  which  the  parachute  may  be 
required  to  fulfil.  Some  of  the  factors  which  will  influence  the  choice  of  design  are  outlined  below. 

2.2.1.1  Opening  Characteristics  The  speeds  at  which  parachute  canopies  are  required  to  deploy  and  to  inflate 
strongly  influence  the  maximum  structural  loads  which  they  must  be  designed  to  withstand.  On  strength 
considerations,  <f  inflation  is  required  at  high  equivalent  air  speeds,  ribbon  parachute  canopies,  such  as  that  shown  in 
Fig.2.4,  are  almost  exclusively  chosen. 


2.21.2  Is  ‘Drive*  Required?  Drive,  defined  in  Section  21,  may  be  a  strageric  requirement,  as  it  is  for  airborne 
forces  parachutes,  or  it  may  be  undesirable,  as  would  be  the  case  for  dropping  stores  by  parachute  into  a  confined  zone. 
As  will  be  outlined  in  Section  23  X  a  parachute  canopy  possessing  drive  is  likely  to  be  strongly  statically  stable  in 
pitch. 

When  drive  is  required,  a  gliding  parachute  ctnopy  must  be  adopted.  But  where  drive  would  be  an  undesirable 
characteristic  a  conventional  parachute  canopy  is  tsed  instead. 

Since  many  gliding  parachutes  are  effectively  inflatable  low  aspect  ratio  wings,  their  aerodynamic  characteristics 
vary  with  the  angle  of  attack  and  wing  aspect  ratio  in  a  manner  typical  of  these  ratio  wings.  Gliding  parachutes  ate 
limited  by  control  considerations  to  a  maximum  aspect  ratio  of  about  3:1,  giving  a  gentle  stall  and  a  correspondingly 
slow  increase  in  drag  coefficient  The  ratio  of  lift  to  drag  is  low;  at  the  present  state  of  the  art  about  3:1  is 
characteristic  but  higher  ratios  are  attainable  with  the  more  advanced  design  of  swept-wing  closed-cell  ram-air  gliding 
pasrachutcs  described  in  Section  9.5.1.  The  characteristic  variation  of  lift  coefficient  with  angle  of  attack  fur  ram-air 
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gliding  parachute  canopies  having  aspect  ratios  varying  from  1:1  to  3:1  is  shown  in  Fig.  2.5.  These  experimental 
results  were  quoted  by  Ungardu  from  an  earlier  report  by  Nicolaides*4. 


Fig.  IS  LIFT  COEFFICIENT  VARIATION  WITH  ANGLE  OF  ATTACK  FOR 
RAM-AIR  GLIDING  PARACHUTES 


Fig.  2.6  DRAG  COEFFICIENT  VARIATION  WITH  ANGLE  OF 
ATTACK  FOR  RAM-AIR  GLIDING  PARACHUTES 


In  Fig.  2.6,  die  corresponding  drag  coefficient  variation  with  angle  of  attack  and  aspect  ratio  is  shown  for  iam-air 
gliding  parachutes. 

22.13  What  Value  of  Cj.  is  Required?  It  will  be  shown  in  Section  2 33  that  for  a  given  canopy  size  and 
payload  mass  the  rate  of  descent  of  a  parachute  decreases  as  the  magnitude  of  the  tangential  force  coefficient  Cj 
increases.  In  general,  as  high  a  value  as  is  pmcticable  is  desirable  for  Cr.  At  zero  angle  of  attack,  Cr  is  equal  to  the 
drag  coefficient  Q,.  Some  typical  characteristics  for  the  variation  of  Cp  and  Cj  with  angle  of  attack  ate  shown  in 
Figs.2.6  and  2.7.  They  are  seen  to  be  dependent  on  the  porosity  of  the  parachute  canopy,  a  property  which  will  be 
discussed  in  Section  2.4.  Whereas  in  Fig£6  characteristics  for  a  typical  gliding  parachute  have  been  given,  in 
Fig£.7  they  are  shown  for  flat  circular  parachute  canopies,  so  called  because  when  these  canopies  are  spread  out  flat 
on  a  plane  surface  they  are  circular  is  shape 


Fig.  2.7  TANGENTIAL  FORCE  COEFFICIENT 
VARIATION  WITH  TANGLE  OF  ATTACK 


Imoorous  Canopy 

Medium  Porosity  0.40 ' 
Canopy 

Veiy  Porous  Canopy  0.20  ■ 


-  FLAT  CIRCULAR  CANOPIES 
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2.2. 1.4  What  Cr  Characteristic  is  Required?  It  will  be  shown  in  Section  23.1  that  the  attitude  of  the  canopy 
when  in  equilibrium  is  detennined  by  the  angle  of  attack  at  which  C*  is  equal  to  zero.  Further,  the  condition  for  a 


•  FLAT  CIRCULAR  CANOPIES  ref  2.2  . 


■Sill 

gill 


■ 

i 

■mill 

I 

i 

rmiil 

nsmsm 

♦*  Imporous  Canopy  "•*  Medium  Porosity  •"  Very  Porous  Canopy 
Canopy 


Very  Porous  Canopy 


Fig.  2.8  NORMAL  FORCE  COEFFICIENT  VARIATION  WITH 
ANGLE  OF  ATTACK 

^  j 


parachute  to  be  statically  stable  in  pitch  is  shown  in  Section  232  to  be  that  when  in  equilibrium,  d<yda  must  be 
positive.  Almost  any  pnchute  canopy  will  descend  stably  if  only  it  is  made  sufficiently  porous  for  example,  ribbon 
patachute  canopies,  being  highly  porous,  characteristically  display  strong  static  stability  in  pitch. 

Typical  characteristics  showing  the  variation  of  C*  with  angle  of  attack  for  flat  circular  parachute  canopies  are 
shown  in  Fig.2.8.  They,  too,  are  strongly  porosity  dependent.  For  these  and  other  parachute  canopy  shapes,  as 
porosity  is  increased  not  only  does  the  tangential  force  coefficient,  Cr  markedly  decrease,  resulting  for  a  given  canopy 
shape  and  size  in  » increased  descent  velocity,  but  as  very  porous  canopies  inSa*  they  may  exb&it  spudding,  defined 
by  Brown'*  as  a  tendency  for  the  open  canopy  to  coBapie  to  a  fora  in  which  the  open  diameter  lies  between  one- 
third  and  one-quaner  of  the  fuUy<open  diameter;  she  shape  of  the  collapsed  parachute  canopy  then  resembling  that  of  a 
squid.  The  phenomenon  of  squidding  is  further  discussed  in  Section  52. 
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2.2.2.  Variation  of  Aerodynamic  Coefficients  with  Reynolds  Number 

Published  dam  on  the  variation  of  parachute  aerodynamic  coefficients  with  Reynolds  nianber  are  almost  entirely 
concerned  with  drag  coefficient  variation.  Provided  that  the  Reynolds  number,  baaed  on  parachute  canopy  nominal 
diameter,  is  greater  than  about  10*,  Figs.  2.9a  and  2.9b  indicate  that  iiule  variation  in  dreg  coefficient  with  Reynolds 
number  occurs.  By  inference,  at  these  Reynolds  numbers  little  variation  of  other  aerodynamic  force  coefficients  with 
Reynolds  number  is  anticipated.  In  R*2.9a,  comparing  corresponding  data  for  other  bluff  bodies,  Le.  a  sphere  and  a 
circular  disc  normal  to  the  flow,  k  is  evident  that  at  high  Reynolds  numbers  the  boundary  layers  detach  and  near  the 
leading  edge  of  these  bodies  flow  separation  occurs. 


Since  for  these  bluff  bodies  there  is  little  or  no  boundary  layer  variation  with  Reynolds  number,  the  resulting 
variation  m  aerodynamic  coefficients  with  Reynolds  number  is  also  small.  In  Figs.  2.9a  and  2.9b  cross  parachute 
canopies  canopies  are  referred  to.  These  are  manufactured  from  two  rectangular  strips  of  material  joined  together  to 
produce  a  cross  or  a  cruciform  shape,  appearing  as  shown  in  Fig.  2.10. 
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As  Fig,2.9b  indicates,  for  Reynolds  numbers  above  10*.  smalt  variations  with  Reynolds  number  in  the  drag 
coefficient  of  parachutes  can  occur,  but  these  variations  are  of  little  practical  significax*.  However,  ftom  free  flight 
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2.2.3  Variation  of  Aerodynamic  Coefficients  with  Mach  Number 

With  bluff  bodies  such  ms  parachute  canopies  substantial  variations  in  drag  coefficients  with  Mach  number  are  to 
be  anticipated.  This  variation  is  strongly  dependent  on  both  the  canopy  shape  and  the  local  Reynolds  number.  Fig. 
2.11a  is  generally  illustrative  of  the  substantial  drag  coefficient  dependence  on  Mach  number,  a  dependence  which 
develops  in  subsonic  flows  as  Fig.  2.11b,  describing  the  characteristics  of  a  parachute  cluster,  illustrates.  Hyperflo 
parachute  canopies,  whose  characteristics  are  shown  in  Fig.  2.11a,  were  flexible  ribbon  concepts  which  were 
specifically  designed  for  supersonic  operation  by  the  Cook  Research  Laboratories. 


In  order  to  establish  high  Mach  number  aerodynamic  characteristics  of  parachute  canopies,  it  is  essential  to  use 
either  ve^r  large  high-speed  wind  tunnel  facilities  to  obviate  blockage  effects  or  else  to  flight  test  full-scale  canopies. 
The  blockage  constraint  of  parachute  models  in  wind  tunnels  is  discussed  in  Section  6.3. 


*  »  \T 

V  ' 


2.3  STEADY-STATE  FLIGHT  MECHANICS 

2.3.1  Equilibrium 

2.3.1.1  The  Conventional  Parachute  Bv  resolving  and  taking  the  moments  of  the  external  forces  acting  on  the 
conventional  parachute  shown  in  Fig.2.12  the  conditions  for  its  equilibrium  can  be  established.  As  the  figure  shows, 
the  angle  0  is  the  inclination  of  the  parachute  axis  of  symmetry,  defined  in  Section  2.1,  with  the  vertical.  Initially, 
the  drag  of  the  payload  is  considered  as  negligible  compared  with  that  of  the  canopy. 


O  Origin  of  co-ordinate  system 
G  System  centroid 
cp  Canopy  centre  of  pressure 
Vg  Resultant  velocity  at  origin 
N  Normal  aerodynamic  force 
component 

T  Tangential  aerodynamic  force 
component 
a  Angle  of  attack 
8  Inclination  angle  of  the 
parachute 


positive 


V* 


Fig.  2.12  EQUILIBRIUM  CONDITIONS 
FOR  CONVENTIONAL  PARACHUTES 


Resolving  in  the  Ox  direction: 


Resolving  in  the  Oz  direction: 


mg  coa  0 


mg  sin  0  - 


Taking  moments  about  the  system  centroid,  G: 


M0  =  -N(x,-*«)  =  0  •  (2.14) 

Then  since  (x,-xj  is  non-zero,  for  equilibrium  the  normal  aerodynamic  force  component.  N  must  be  zero. 

Equation  2.13  then  shows  that  since  mg  is  necessarily  non-zero,  at  equilibrium  the  angle  0  must  be  zero.  Under 
this  condition  equation  2.12  shows  that  the  tangential  force  component,  T  equals  mg.  Hence,  when  it  is  in 
equilibrium  a  conventional  parachute  descends  with  its  axis  of  symmetry  vertical,  at  such  an  angle  of  attack,  a  that  no 
normal  aerodynamic  force  component  is  developed  on  it 

When  the  drag  of  the  payload  is  not  negligible  when  compared  with  that  of  the  canopy,  at  equilibrium  the  normal 
force  component  is  small  and  positive.  Under  these  conditions,  a  stable  canopy  descends  at  a  small  positive  angle  of 
attack,  dekribing  a  coning  motion  with  the  semi-apex  angle  of  the  cone  equal  to  this  angle  of  attack.  In  order  to 
minimise  this  coning  motion  the  drag  of  the  payload  must  be  small  compared  with  that  of  the  canopy. 

2.3. 1.2  The  Gliding  Parachute  For  exattly  the  same  reason  as  for  a  conventional  parachute,  if  the  drag  of  the 
payload  attached  to  a  gliding  parachute  is  neglected  then  when  uus  parachute  is  in  equilibrium  it  descends  so  that  its 
axis  Oz  is  vertical  This  axis  is  drawn  through  both  the  canopy  centre  of  pressure,  cp,  located  on  the  aerofoil  section 
chord  line  at  about  the  quarter-chord  position  and  the  parachute-payload  system  centroid,  G.  From  the  origin  O, 
selected  on  the  axis  Oz,  the  axis  Ox  extends  at  right  angles  to  Oz  in  the  plane  of  symmetry  and  in  the  same  sense  as 
that  of  the  parachute's  resultant  velocity  V*. 
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Fig.  2.13  has  been  drawn  in  (he  equilibrium  position.  In  this  figure  the  line  OP  has  been  drawn  through  the 
origin  O  and  perpendicular  to  the  aerofoil  chord  line.  The  adjustable  rigging  angle  $  lies  between  OP  and  the  axis  Oz. 

In  reference  2.1  the  symbol  y  is  used  for  the  angle  of  climb  of  an  aircraft.  Throughout  The  Aerodynamics  of 
Parachutes  the  angle  of  descent,  y4  will  be  adopted  instead,  where 


Y<  =  -T-  (2.15) 

Using  the  symbols  shown  in  Fig.2.13  and  applying  equilibrium  conditions  corresponding  to  those  adopted  for 
conventional  parachutes  in  equations  2.12  to  2.14,  including  neglecting  the  drag  of  the  payload,  by  resolving  forces  in 
the  Ox  direction: 


L  sin  y4  D  cosy*  =  0  (2.16) 

by  resolving  forces  in  the  Oz  direction: 

ug  -  L  cos  y4  •  D  sin  Ye  =  P  (2.17) 

and  by  taking  moments  about  the  centroid, G: 

M0  =  -  (L  sin  y4  -  D  cos  y4  M*,-  *,)=  0.  (2.18) 


O  Origin  of  co-ordinate  system 
G  System  centroid 
cp  System  centre  of  pressure 
V  Resultant  velocity  at  origin 
L  Lift  aerodynamic  force 
component 

D  Drag  aerodynamic  force 
component 
a  Angle  of  attack 
♦  Rigging  angle 
Y-  Angle  of  descent 


perpendicultr 
tc  chord  line 

Fig.  2.13  EQUILIBRIUM  CONDITIONS 
mg  FOR  GLIDING  PARACHUTES 
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From  the  inset  diagram  in  Fig.  2.13: 


tan  y, 

=  VI L  = 

1  • 

(2.19a) 

L/D 

and 

cos  Ye 

*  L/mg  . 

(2.19b) 

From  the  conditions  established  in  Fig.  2.13,  at  equilibrium  the  angle  of  descent  x-is  related  to  the  angle  of  attack 
a  by: 


'it 


a  +  4 


(2.20) 


IS 


If  the  drag  of  the  payload  is  not  neglected,  then  at  equilibrium  the  Oz  axis  of  the  gliding  parachute  is  inclined  at  an 
angled  to  the  vertical,  so  that 


a  +  f  «  9  +  .  (2.21) 

By  adding  the  payload  drat  to  the  forces  shown  in  Fig.  2.13  it  can  readily  be  seen  that  the  condition  for 
equilibrium  is  that  the  axis  Oz  is  inclined  at  a  small  negative  angle  9  to  the  vertical  As  this  angle  is.  in  general 
small  equation  2.20  is  approximately  valid  even  when  it  is  necessary  to  take  the  drag  of  the  payload  into 
consideration.  Thus; 


y4  -  a  +  $  .  (2.22) 

Now,  through  equation  2.19,  y4is  a  function  of  the  angle  of  attack,  a .  Equation  2.22  demonstrates  that  the 
purpose  of  allowing  the  rigging  angle  $  to  vary  is  to  adjust  the  angle  of  attack  at  which  the  gliding  parachute  flies, 
making  possible  a  range  of  equilibrium  ratios  of  lift  to  drag  at  which  the  parachute  descends. 

2.3.2  Static  Stability 

Whether  or  not  a  system  is  stable  is  determined  by  its  response  to  small  displacements  from  its  equilibrium 
position.  Static  stability  is  solely  concerned  with  the  directions  of  the  moments  which  are  developed  on  such  a 
displaced  system,  a  statically  stable  system  being  one  which,  as  a  consequence  of  a  small  displacement  from 
equilibrium,  develops  a  moment  in  a  direction  which  would  restore  the  system  to  equilibrium.  This  concept  has 
nothing  to  say  about  the  equilibrium  of  forces  which  act  on  the  system  after  displacement  or  about  the  frequency  and 
attenuation  or  amplification  of  any  resulting  oscillations,  nevertheless,  it  is  a  very  valuable  concept  in  a  number  of 
fields  including  aircraft  dynamics  and  in  Section  4.3.2  it  is  shown  to  have  particular  merit  when  applied  to  parachute 
dynamic  pitching  motion.  Since  a  system  which  exhibits  static  stability  about  one  axis  need  not  necessarily  do  so 
about  any  other,  it  is  important  to  define  the  axis  about  which  the  stability  erf  a  system  is  being  considered. 


The  concepts  of  static  stability  and  instability  are  illustrated  in  Kig.2.14  by  reference  to  a  pendulum,  hanging  on 
the  end  of  a  light,  straight  rod.  When  this  supporting  rod  hangs  vertically  the  pendulum’s  suspended  mass  is  in 
equilibrium  Work  must  be  done  to  displace  it  from  this  equilibrium  position  and  to  incline  the  supporting  string 
through  an  angle  A9.  Hiving  made  such  a  small  displacement,  the  pendulum  develops  a  moment  AM  about  the 
suspension  point  0.  This  moment  is  in  the  direction  to  restore  the  system  to  the  equilibrium  state  and  is  opposed  to 
that  of  the  displacement.  The  static  stability  of  the  system  is  characterised  by  'he  relationships  that: 

(0.  in  equilibrium,  the  moment  about  O: 


M.(9)  >  0 


(233) 


00-  and  after  a  small  displacement: 


WW 


dM^dO  <  0  .  (2.24'. 

However,  if  the  pendulum  were  inverted,  with  the  supporting  rod  vertical  the  system  is  also  in  equili  brium.  But, 
in  allowing  the  supporting  rod  to  deflect,  work  is  done  by  the  system  on  the  rod.  Having  deflected  through  a  small 
angle  A0  the  moment  then  developed  about  the  suspension  point  0  is  in  the  same  direction  as  that  of  the 
displacement.  Hence,  for  this  invertea  jAndulum  dM^d0  is  greater  than  0. 

In  exactly  the  same  way,  the  criteria  required  for  the  equilibrium  and  the  static  stability  of  a  parachute  in  pitch  are, 
from  equations  2.12  to  2.14  and  2.16  to  2.18,  that 

M0(a)  =  0  (225) 

and  dMc/3a  <  0  .  (2.26) 

By  specifying  moments  about  the  parachute  centroid  it  is  ensured  that  the  moment  consequent  upon  a  small 
disturbance  from  the  equilibrium  state  is  wholly  aerodynamic.  Its  magnitude  and  its  sign  can  then  be  determined 
readily,  solely  from  steady-state  aerodynamic  tests. 

For  a  conventional  parachute,  the  necessary  conditions  for  equilibrium  and  static  stability  in  pitch  are  given  from 
equations  2.14  and  225  as: 


together  with 


M0  (a)  =  -N(x,.x,)  =  0  . 

3M  0/3a  <  0 


or,  from  equation  227: 


dN/da 


Whereas,  for  a  gliding  parachute,  from  equations  2.18  and  225: 

M0  =  -  [L  sin  (a  +  $)-D  cos  (a  +  6)]z0  =  0  .  (230) 

together  with  dMG/3a  <  0  .  (231) 

Some  typical  steady-state  aerodynamic  pitching  moment  characteristics  for  conventional  parachutes  have  been 
taken  from  the  Recover  Systems  Design  Guide1  •*  and  are  shown  in  Fig.  2.15. 

f  Pitching  Moment  Coefficient  Cm(,  ^ 
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The  aerodynamic  pitching  moment  coefficient  Cm  is  defined  in  a  similar  way  to  the  force  coefficients  in 
equations  2.1, 2.2, 23  and  2.6  by: 


C.  «  M  (2.32) 

‘/jpV’SJD, 


the  nominal  diameter  J),  being  defined  in  equation  2.4. 

Since  pitching  moments  are  the  products  of  the  normal  force,  N  and  an  appropriate  moment  arm,  when  tabulating 
experimental  results  it  is  important  to  specify  precisely  the  location  of  the  axis  about  which  the  pitching  moment  has 
been  measured  .  As  equation  2.25  indicates,  the  system  centroid,  G  is  often  an  important  location  about  which 
pitching  moments  are  required  but  as  this  is  dependent  on  the  mass  of  the  payload  which  is  suspended  from  the  canopy 
its  location  may  not  be  known  with  precision  at  the  time  of  performing  aerodynamic  experiments  on  the  parachute 
canopy.  In  much  of  the  established  pitching  moment  coefficient  data  for  parachute  canopies  the  point  about  which 
this  moment  has  been  measured  is  unspecified.  Often  it  is  the  suspension  line  confluence  point  which,  for  all 
practical  purposes,  can  be  considered  to  be  the  system  centroid.  This  lack  of  precision  can  make  these  data  unreliable. 
In  his  experimental  data  Doherr110  has  overcome  this  problem  by  specifying  the  canopy  centre  of  pressure  location 
relative  to  the  canopy  hem  line  at  the  skirt  periphery. 

The  parachute  canopies  whose  characteristics  are  shown  in  Fig.2.15  are  all  in  equilibrium  at  zero  angle  of  attack. 
However,  of  the  five  canopies  illustrated  only  those  that  are  cross-shaped  are  statically  stable  at  this  angle,  the  other 
three  canopies  illustrated  exhibiting  simultaneous  equilibrium  and  static  stability  at  certain  positive  angles  of  attack. 
For  example,  the  flat  circular  canopy  is  in  equilibrium  and  is  statically  stable  in  pitch  at  an  angle  of  attack  of  about 
20  degrees.  These  pitching  moment  characteristics  are  skew-symmetric  about  zero  angle  of  attack,  thus  the  fiat 
circular  canopy  described  is  also  in  equilibrium  and  is  statically  stable  at  -  20  degrees  angle  of  attack. 

Both  of  these  solutions  are  within  the  Oxz  plane.  The  flat  circular  canopy  exhibits  the  required  equilibrium  and 
stability  characteristics  whenever  it  flies  with  its  axis  of  symmetry  vertical  and  the  resultant  relative  airflow  lies  on 
the  surface  of  a  cone  whose  axis  is  the  canopy  axis  of  symmetry  and  whose  semi-apex  angle  is  20  degrees.  Since  no 
particular  direction  is  preferred,  during  descent  the  parachute  oscillates  through  approximately  ±  20  degrees.  In  the 
literature  such  a  parachute  is  referred  to  as  an  unstable  parachute . 

On  the  other  hand,  the  cross  (or  cruciform)  parachute  canopy  of  3.8:1  arm  ratio  shown  in  Fig.  2.10,  is  referred  to 
as  a  stable  parachute .  It  is  in  equilibrium  and  is  statically  stable  in  pitch  at  zero  angle  of  attack,  thus  any  disturbance 
from  this  equilibrium  state  will  be  attenuated.  Since  the  angles  through  which  stable  parachutes  oscillate  depend  on 
the  amplitude  of  the  forces  which  disturb  them  from  equilibrium,  during  descent  they  cannot  be  stated  with  any 
exactitude.  However,  the  oscillatory  motion  which  ensues  is  often  heavily  damped,  for  reasons  to  be  explained  in 
Section  4.3.2.  Thus  the  observed  oscillatory  motion  can  be  quite  minimal  and  in  the  literature  it  is  customary, 
though  inaccurate,  to  specify  these  small  pitching  angles  through  which  srable,  as  well  as  unstable,  parachutes 
oscillate. 

It  is  clear  from  Fig.2.15  that  die  steady-state  aerodynamic  characteristics  of  cross  parachutes  are  functions  of  arm 
ratio,  defined  in  Fig.  Z10  as  the  constructed  length  to  width  ratio  of  one  of  the  two  canopy  arms.  Normal  force 
coefficient  variations  with  angle  of  attack  for  a  variety  of  arm  ratios  are  shown  in  Fig.2.16.  Equations  22  and  2.29 
indicate  that  a  necessary  condition  for  a  conventional  parachute  to  be  stable  in  pitch  is  that  ACJdo.  must  be  positive  at 
a=0\  For  the  cross  parachute  canopies  shown  in  Figs.  2.16,  this  condition  is  satisfied  by  appropriate  combinations 
of  arm  ratios  and  fabric  porosities,  for  example  the  curves  drawn  in  Fig.  2.16a  are  for  im porous  canopies,  showing 
that  a  condition  for  static  stability  in  pitch  is  that  the  arm  ratio  of  imporous  cross-shaped  canopies  should  exceed 
3.0:1. 
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Of  course,  cross  parachutes  are  not  the  only  canopies  which  exhibit  static  stability  in  pitch.  As  Figs.  2.3  and 
2.15  have  already  indicated,  a  tendency  to  parachute  stability  is  a  consequence  of  canopies  increasing  in  porosity. 
However,  as  Fig.  2.5  indicates,  thij  is  accompanied  by  a  corresponding  reduction  in  drag  coefficient  and  as  explained 
in  Section  2.3.3,  leads  to  deteriorating  descent  characteristics.  A  considerable  experimental  test  programme  on 
canopy  stability  was  undertaken  in  1962  by  Heinrich  and  Haak^*.  A  number  of  its  results  arc  given  in  Chapter  4 
of  the  1963  Parachute  Design  Guide*7. 


stability  in  pitch  of  certain  solid  textile  canopies  and  that  of  a  number  of  slotted  parachute  canopies.  In  the  reference 
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quoted  he  indicated  that  there  was  some  Reynolds  number  dependency  in  his  results.  He  also  showed  that  substantial 
interference  to  the  canopy’s  aerodynamic  characteristics  is  often  caused  by  the  presence  of  the  payload  forebody. 

Aerodynamic  characteristics  of  unstable  conventional  parachute  canopies  can  be  considerably  changed  by  removal 
of  a  panel  or  a  portion  of  a  panel.  The  consequent  loss  in  axial  symmetry  results  in  the  parachute  acquiring  a 
horizontal  component  of  velocity  or  ‘drive’  during  descent  and  thus  becoming  a  gliding  parachute.  Jorgensen  and 
Cockrell*1*  have  demonstrated  that  the  resultant  relative  airflow  is  then  at  a  high  and  staticaliy-stable  angle  of  attack, 
thus  the  parachute  acquires  a  satisfactory  descent  performance. 

2.3.3  Steady  Descent 

In  a  steady  equilibrium  descent  it  has  been  shown  in  Section  2.3.1  that,  provided  the  payload  drag  is  negligible 
compared  with  that  of  the  parachute  canopy,  the  attitude  of  the  axis  of  symmetry  Ox  for  a  conventional  parachute  and 
of  the  axis  Oz  for  a  gliding  parachute  is  vertical.  Under  these  conditions,  for  conventional  parachutes  equation  2.12 
reduces  to: 


mg  =  T  (2.33) 

If  the  parachute  is  stable  then  during  descent  its  angle  of  attack  is  zero  and  from  equation  2.7  C,  =  Q,.  Hence, 
from  equation  2.6,  for  a  stable  cwiventional  parachute: 

mg  =  ’/jPVjS.Cb  (2.34) 

From  this  equation  it  is  evident  that  the  descent  velocity  V0  is  inversely  proportional  to  the  square  root  of  the 
drag  coefficient  Q,.  It  is  also  proportional  to  the  square  root  of  the  parachute  and  payload  weight  and  as  the  former  is 
generally  negligible  compared  with  the  latter  the  descent  velocity  can  be  considered  to  be  a  function  of  the  weight  of 
the  payload.  It  decreases  as  the  ground  is  approached  and  the  air  density,  p  increases  and  it  is  inversely  proportional  to 
the  square  root  of  the  canopy  surface  area,  S.,  By  measuring  the  descent  velocity  VD  and  knowing  the  other 
parameters  the  drag  area  of  the  parachute,  CoS.,  can  thus  be  calculated. 

It  is  normal  practice  to  apply  equation  234  to  unstable  as  well  as  to  stable  conventional  parachutes.  This  is  then 
an  average  determination  made  for  a  body  whose  angle  of  attack  during  descent  will  vary.  It  may  therefore  differ  in 
magnitude  from  wind  tunnel  derivations  made  at  fixed  angles  of  attack.  Although  with  stable  parachutes  this 
difference  is  not  marked  with  unstable  canopies  k  could  be  significant 
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Fig.  2.18  STEADY  DESCENT  OF  A  GLIDING  PARACHUTE 
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For  gliding  parachutes,  from  equation  2.16: 

mg  =  L  cos  ya  +  D  sin  yt  (2.35) 

Thus,  from  the  definitions  of  Q,  and  Qj-in  equations  2.5  and  26  and  the  expression  for  cos  %  in  equation  2.19b: 
mg  *  l/2pViS.[Ct/  coa  yj 


(2.36) 
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and  thus 


V»  -  (2.37) 

l  pS.CL  J 

As  Fig.2.18  shows,  the  horizontal  and  vertical  components  of  the  parachute's  resultant  velocity,  V*  are  given 
respectively  by  u  and  w,  where 

a  =  V ,  cos  (2.38a) 

and  w  =  V,  stay,  .  (2.38b) 

From  equation  2.37,  the  velocity  down  the  glide  path,  V*  increases  with  increasing  altitude  and  increasing  wing 
loading,  mg/S..  But  for  a  given  height  loss,  the  horizontal  distance  travelled  is  a  function  of  the  angle  of  descent  % 
which,  equation  2.21  shows,  is  solely  a  fuction  of  the  gliding  parachute's  aerodynamic  characteristic,  CJCo- 
2.3.4  Froude  Number,  F 
The  Froude  number  F  for  a  parachute  is  defined  as: 

F  =  Vl  /  D.g  (239) 

Since  it  relates  the  descent  speed  in  a  given  gravitational  field  to  the  required  size  of  the  parachute  canopy,  this 
similarity  parameter  is  a  very  useful  performance  number.  It  is  not  wholly  an  aerodynamic  characteristic  of  the 
parachute,  since  it  is  a  combination  of  the  system's  drag  coefficient  with  its  mass  ratio,  R,: 

F  «  1/(R„C0)  (2.40) 


where  the  mass  ratio  R.  is  a  measure  of  the  ratio  of  air  mass  enclosed  in  the  fully-inflated  parachute  canopy  to 
the  payload  mass  m,,  thus: 


R»  =  pV/m,  . 


(2.41) 


In  equation  2.41  the  symbol  V  denotes  the  representative  displaced  volume  of  an  immersed  body.  For  a  parachute 
canopy  the  representative  displaced  volume  is  considered  to  be  that  of  a  hemisphere  which  has  a  diameter  equal  to  the 
canopy  nominal  diameter  Dr  The  concept  of  the  representative  displaced  volume  is  further  discussed  in  Section  43.1. 


This  relationship  between  the  Froude  number  and  the  mass  ratio  for  parachutes  is  illustrated  in  Fig.  2.19,  devised 
by  Doherr  as  a  means  of  illustrating  the  descent  characteristics  for  four  different  parachute  classes.  It  illustrates 
equation  2.40  which  shows  that  for  a  given  drag  coefficient  the  Froude  number  F  is  inversely  proponional  to  the  mass 
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the  latter  is  determined  from  the  product  of  the  mouth  area  and  the  rate  of  canopy  descent,  VD.  The  volumetric 
outflow  can  be  estimated  float  (be  air  velocity  U  discharging  through  an  area  which,  because  of  cut  outs,  is  less  than 
the  nominal  area  of  the  canopy ,  together  with  the  discharge  determined  from  the  product  of  the  descent  velocity  and  the 
cut-out  area.  The  unknown  velocity  U  at  a  pressure  difference  dp  equal  to  l/2pV£  can  be  calculated  from  equation 
2.44,  using  previously  determined  values  of  Iq  and  kg  for  the  canopy  fabric. 

2.5  STOKES-FLOW  PARACHUTE  SYSTEMS 

When  the  application  is  one  in  which  very  low  descent  velocities  or  very  high  altitude  performance  is  required,  the 
resulting  low  Reynolds  number  leads  to  in  a  very  different  parachute  concept,  that  of  a  Stokes-flow  pcrachiue  system, 
which  Mihora"3  has  described. 

If  the  Reynolds  number  is  so  small  that  inertia  forces  are  negligible  compared  with  those  caused  by  viscosity  the 
parachute's  drag  mechanism  and  therefore  its  basic  shape  mil  be  changed.  To  achieve  a  low  descent  velocity  the 
canopy  must  be  of  an  extremely  lightweight  construction  of  open-mesh  fabric  or  else  formed  flora  very  small  diameter 
widely-spaced  filaments.  Its  nominal  drag  coefficient,  based  on  its  total  projected  area,  could  well  exceed  1.0.  While 
its  overall  Reynolds  number  might  exceed  10  000  when  based  on  the  canopy  diameter,  if  it  were  based  on  the 
filaments’  diameter  it  could  well  be  less  than  0.5,  a  necessary  criterion  for  Stokes  flow.  Because  of  the  large  voids 
present  in  such  a  canopy  it  can  be  packed  into  a  very  mall  container. 

Niederer**  describes  studies  for  a  Stokea-flow  decelerator  system,  designed  to  provide  subsonic  rates  of  descent  at 
altitudes  up  to  90  km  (56  miles).  Experiments  were  conducted  in  a  high  althude  test  chamber  in  which  either  full  scale 
Reynolds  number  or  full  scale  Knudsen  number  could  be  reproduced  bat  not  both  at  the  same  time.  The  Knudsen 
number  is  the  ratio  X/L,  where  X  is  the  mean  flee  path  for  the  molecules  of  the  gas  in  which  parachute  filaments, 
having  a  typical  linear  dimension  L,  are  immersed.  One  of  Niedeter’s  sample  designs  is  for  a  Stokes-flow  parachute 
which  weighs  0.44  Newtons  (0.1  Ibf)  and  which  supports  a  4.0  Newton  (0.9  Ibf)  payload.  At  an  altitude  of  45.8 1cm 
(28.5  miles)  its  me  of  descent  would  be  24.4  m/sec.  (80  flftec.). 

Although  suited  to  high  altitude  applications  because  there  the  kinematic  viscosity  is  high  and  a  low  Reynolds 
number  is  thus  more  readily  achieved.  Stokes-flow  parachutes  can  be  used  successfully  at  much  lower  altitudes.  An 
example  in  nature  of  the  Stokes-flow  parachute,  found  at  sea  level,  is  the  dandelion  seed. 
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APPENDIX  TO  SECTION  2 

SOME  COMMONLY  -ADOPTED  SHAPES  FOR  CONVENTIONAL  PARACHUTE  CANOPIES 
A  Summary  of  Aerodynamic  Characteristics,  drawn  from  references  1.7, 1.9  and  2.7. 

2A.1  SOLID  TEXTILE  CANOPIES 


Description 

FLAT  CIRCULAR  CANOPY 
cyanucnoM  XMt»A?ic 


Usaee 


A  low  cost,  light  weight 
canopy  which  is  suitable  for 
deployment  at  low  dynamic 
pressures.  It  is  most  widely 
used  by  personnel  and  for 
cargo  recovery. 


Aerodynamic 

Characteristics 


Exhibits  instability  in 
pitch,  possessing  a  conic 
oscillation  of  some  20-30' 
semi-apex  angle.  Clusters 
of  flat  circular  canopies  do 
not  exhibit  this  characteristic, 
since  each  canopy  in  wc 
cluster  is  no  longer  at  O' 
angle  of  attack. 
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Description 

EXTENDED  SKIRT  CANOPY 


Usage 


Suitable  for  ibe  airdrop  of 
material  and  for  the  recovery 
of  drones. 


Frequently  used  as  a  pilot 
parachute,  for  extraction 
purposes.  Has  been  used 
with  both  very  low  and 
very  high  dynamic  pressures, 
the  Mach  number  ranging 
from  0  to  3.0 


2A 2  SLOTTED  TEXTILE  CANOPIES 


Description 
CROSS  CANOPY 


Usage 


coMsnKicnoN  scmmatk 


WLATtOfhOriLS 


An  easily-constructed  low 
cost  parachute  which 
is  frequently  used  for  the 
deceleration  of  ground 
vehicles  and  aircraft.  Its 
applications  also  include 
weapon  deceleration  and 
stabilisation. 


I 


Aerodynamic 

Characteristics 


This  canopy  has  a  higher 
drug  than  a  flat  circular 
canopy.  It  also  has  a 
longer  opening  time,  but  it 
develops  correspondingly 
lower  opening  forces. 


Has  a  high  geometric 
porosity  (15*  to  30%), 
which  gives  it  good 
stability  in  pitch.  It  is  also 
very  reliable  in  inflation. 


Aerodynamic 

Characteristics 


Good  drag  characteristics 
and  provided  that  its 
arm  ratio  and  porosity 
are  properly  selected  it 
has  an  excellent 
Stability  in  pitch. 
However,  it  displays 
a  tendency  to  rotate  about 
its  axis  of  symmetry. 
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Description 
RIBBON  CANOPY 


RING  SLOT  CANOPY 


Usage 


Aerodynamic 

Characteristics 


A  very  strong  canopy  which 
is  tberefbre  used  for  aircraft 
end  vehicle  deceleration  at 
relatively  high  dynamic 
pressures  'up  to  270  kN/m1 
(5  700  lbf/ft^  1 


Has  a  lower  drag 
efficiency  thin  a  solid 
textile  canopy  but  it 
also  has  an  excellent 
stability  in  pitch. 


This  is  like  a  ribbon  canopy  but 
with  much  wider  hauex.  Since 
these  lead  to  it  being  relatively 
che^»  to  manufacture  it  is 
fieonenthr  used  for  high  dynamic 
pressure  cargo  delivery  and 
aircraft  deceleration. 


Its  aerodynamic 
characteristics  ate 
comparable  with  those 
of  the  ribbon  canopy. 


2J  GLIDING  PARACHUTE  CANOPIES 
Hie  Ram-Air  Canopy,  drawn  from  references  15, 2.3  and  9.1. 


Description 


Usage 


Used  when  ‘drive*  is  requited. 

Of  the  high-glide  parachutes 
described  in  ref.  15  the  most  common 
are  canopies  which  are  ram-air  inflated. 
They  have  s  variety  of  commercial 
names. 

They  are  characterised  by  aerofoil 
cross-sectional  shapes  and  pbmforms 
whose  aspect  ratios  are  as  large  as 
P«**ible  on  aerodynamic  grounds. 
However,  in  order  that  good  control  is 
maintained  they  are  limited  to  a 
m**in>wn  of  about  3:1. 


Aerodynamic 

Characteristics 


The  angle  of  descent  y4  is 
inversely  proportional  to  the 
ratio  l/D  (equation  2.19a). 
Typically,  (L/D)^,  is  about 
3:1,  but  sec  Section  9.5.1 


Rear  Suspension  Lines  have  not  been  shown 
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3.TRAJECT0RY  DYNAMICS 

Tbe  object  of  Trajectory  Dynamic*  is  to  predict  the  flight  path  for  the  system  which  comprises  the  combination  of 
a  parachute  and  payload.  In  principle,  provided  the  parachute's  significant  aerodynamic  characteristics,  in  particular  the 
canopy  drag  area  CoS.,  defined  in  Sec  lion  2.3.3,  are  known  as  a  function  of  time,  then  the  appropriate  equations  of 
motion  for  the  parachute-payload  system  can  be  written  and  approximate  numerical  solutions  obtained.  In  the  cr.urse 
of  the  time  over  which  these  equations  are  valid  tbe  parachute  might  be  undeployed,  undergoing  deployment,  inflating 
or  be  fully  inflated. 

For  all  but  the  simplest  models  of  this  system,  tbe  equations  of  motion  are  not  straightforward.  Consequently, 
their  solutions,  too,  are  complex. 

3.1  THE  TRAJECTORY  SYSTEM.  AXES  SYSTEMS.  DEGREES  OF  FREEDOM 

The  nature  of  the  trajectory  system  which  it  is  necessary  to  consider,  the  axis  system  to  be  adopted  and  the  number 
of  degrees  of  freedom  which  are  required  in  tbe  consequent  analysis  depend  on  both  the  input  data  which  are  available 
and  the  complexity  of  the  solution  which  is  required.  The  number  of  equations  required  to  describe  the  motion  of  the 
system  depend  on  the  number  of  degrees  of  freedom  which  the  system  possesses.  A  single  body  which  can  move 
freely  in  a  plane  possesses  three  degrees  of  freedom.  It  requires  three  independent  variables  to  define  its  position 
relative  to  fixed  axes,  two  co-ordinates  to  locate  a  chosen  point  in  tbe  body  and  a  further  co-ordinate  to  orientate  the 
body.  Tbe  total  number  of  co-ordinates  required  to  specify  the  configuration  of  n  unconnected  bodies  is  3n.  but  if 
these  bodies  are  connected  together  by  various  mechanisms  then  degrees  of  freedom  «e  lost 

In  order  to  establish  tbe  equations  of  motion  it  is  important  to  choose  carefully  the  frames  of  reference  which  are 
to  be  adopted.  In  Newton's  laws,  on  which  rigid-body  dynamics  is  based,  all  motion  is  ultimately  considered  relative 
to  a  stationary  reference  frame.  These  laws  state  that  the  external  forces  and  moments  which  act  on  a  system,  together 
with  the  system's  inertia  forces  and  moments,  are  in  a  state  of  equilibrium.  Inertia  forces  and  moments  are  the 
reversed  rate  of  change  of  the  system's  linear  and  angular  momenta. 

For  most  parachute  applications,  the  earth  can  be  considered  to  provide  an  absolute  frame  of  reference  and  axes 
which  are  fixed  relative  to  the  earth  are  termed  earth-fixed  axes.  When  the  frame  of  reference  is  fixed  relative  to  a 
moving  body,  as  could  occur  in  a  system  consisting  of  a  parachute  canopy  rigidly  connected  to  its  payload,  they  are 
termed  body  axes.  In  general,  on  both  parachute  and  payload  the  aerodynamic  and  the  inertial  forces  and  moments 
which  act  oo  the  system  can  be  determined  relative  to  these  body  axes,  since  they  are  functions  of  the  body  resultant 
velocity  and  its  attitude.  However,  not  only  are  tbe  gravitational  forces  and  moments  determined  relative  to  earth-fixed 
axes  but  the  trajectory  of  the  descending  system  is  ultimately  required  relative  to  these  axes.  Thus,  in  obtaining 
solutions  to  trajectory  dynamics  problems,  it  is  usually  necessary  to  adopt  both  earth-fixed  axes  and  body  axes, 
together  with  the  geometric  transformation  relationships  between  these  two  axes  sets.  Since  the  payload  can  move 
relative  to  the  parachute  it  may  be  necessary  to  adopt  more  than  one  set  of  body  axes,  establishing  some  idealised 
relationship  to  describe  the  mdde  of  coupling  between  tbe  parachute  and  iu  payload. 

There  is  no  universally-agreed  method  of  modelling  the  parachute-payload  system.  It  is  important  to  model  the 
system's  mechanics  in  no  more  complex  a  manner  than  is  appropriate  to  develop  the  required  solution  within  the 
constraints  imposed  by  tbe  available  data.  The  approach  adopted  here  owes  much  to  Purvis11.  Two-degree  of 
freedom  equations  of  motion  are  first  described  for  die  motion  in  two  dimensions  of  a  point  mass  representative  of 
both  the  parachute  and  the  payload,  subjected  to  both  aerodynamic  and  gravitational  forces.  Next,  a  three-degree  of 
freedom  set  of  equations  is  developed  for  two-dimensional  motion  of  a  payload  which  is  simply-connected  to  a 
parachute  of  negligible  mass.  Using  these  equations,  the  consequences  are  considered  of  including  the  parachute 
canopy  mass  in  this  trajectory  dynamics  model.  Tbe  third  case  described  is  a  three-degree  of  freedom  model  for  a 
parachute  canopy  assumed  to  be  rigidly-connected  to  its  payload.  Finally,  reference  is  made  to  publications  which 
describe  trajectory  models  with  six  or  more  degrees  of  freedom  and  to  the  use  of  finite-element  analysis  in  order  to 
model  parachute  canopies  during  their  deployment  and  inflation  phases  as  if  they  were  a  series  of  elastically-connected 
mass  nodes. 

3.2  TWO-DEGREE  OF  FREEDOM  MODEL 

Consider  the  motion  of  a  point  mass  in  the  x-z  plane,  shown  in  Fig.3.1.  As  this  mass  is  considered  to  act  ai  a 
point,  the  system  which  it  represents  cannot  possess  any  moments  of  inertia;  hence  no  moments  arc  exerted  on  it.  In 
the  plane  Oxz,  therefore,  the  system  possesses  only  two  degrees  of  freedom,  translational  motion  in  tbe  direction  O-x 
airi  translational  motion  in  the  direction  0-z. 

On  symmetry  grounds,  the  point  mass  cannot  sustain  any  component  of  aerodynamic  force  normal  to  the  line  of 
flight  Thus,  the  resultant  component  forces  which  act  oo  it  in  tbe  O-x  and  O-z  directions,  F,  and  F„  are 
respectively  given  by: 


F,  s  -{D,+D,)coi  y4 


(3.1) 
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and  F,  =  mg  •  (Dr+D,)sin  y,  (3.2) 

where  m  denotes  the  mass  of  the  system  inclusive  of  any  added  mass,  D,  is  the  drag  developed  by  the  parachute, 
D,  is  the  drag  developed  by  the  payload  and  yt  is  the  angle  of  descent  Added  mass  is  discussed  in  Section  4. 


Measuring  the  displacements  x  and  z  relative  to  the  earth-fixed  axes  O-x  and  O-z,  since 


» 

X  s 

v  cosy,, 

(33) 

and 

* 

X  « 

V  stay. 

(3.4) 

then  the  trajectory  equations  are: 

•• 

mx  = 

-  (l/2)pVx  (CBS,+  CAS() 

(3.5) 

and 

mz  = 

mg  -  (1/2 )pVi  (CDS.+  CAS,) 

(3.6) 

where  CqS.  is  the  parachute  canopy  drag  area,  S(  denotes  the  payload  reference  area  and  CA  the  payload  axial 
force  coefficient. 

From  known  initial  conditions  of  the  parachute  and  payload,  equations  33  and  3.4  give  initial  values  of  x  and  z. 
Then,  for  known  values  of  the  payload  reference  area  and  axial  force  coefficient,  the  system  weight  mg,  the  canopy 
drag  area  as  a  function  of  time,  together  with  the  parachute  system's  initial  altitude  and  velocity,  equations  3.5  and  3.6 
are  soluble.  They  yield  as  functions  of  time  the  system's  horizontal  and  vertical  co-ordinates,  its  velocity  and  its 
flight  path. 

Initial  values  of  x  and  z  can  then  be  determined  from  the  equations  3.5  and  3.6.  These  relationships  are  first-order 
differential  equations  in  x  and  z.  Over  short,  finite  time  increments,  At  their  solutions  for  x  and  z  can  be  used  to 
update  x  and  z.  To  achieve  this  end,  a  number  of  appropriate  numerical  integration  schemes  exist,  of  which  the 
following  is  illustrative: 


(*)i  =  (i)i  ♦  (x)tdt  (3.7) 

and  (*)i  *  (*)i  ♦  (a),  At  .  (3.8) 

3.3  THREE-DEGREE  OF  FREEDOM  MODEL 

When  angular  as  well  as  linear  motion  in  the  Oxz  plane  is  required  for  the  parachute  and  its  payload,  then  a  three- 
degree  of  freedom  model  is  adopted.  Purvis  notes  that  at  this  stage  the  method  adopted  to  represent  the  parachute  and 
the  payload  must  be  considered  carefully  and  lists  four  different  approaches: 


WSWj  *»»w 
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ffl.  *  massless  parachute  is  joined  by  a  massless  rigid  link  to  its  payload,  the  link  being  pin-joined  at  its 
attachment  points; 

(iQ.  the  parachute  and  its  payload  constitute  a  single  rigid  body; 

(Hi).  the  parachute  and  its  payload  each  possess  mass  and  each  constitute  a  single  rigid  system,  the  systems 
being  joined  by  a  massless  rigid  link,  and 

(iv).  the  parachute  is  represented  as  an  elastic  system,  the  payload  as  a  single  rigid  body. 

3.3.1  The  Massless  Parachute  Joined  to  the  Payload 

For  the  point  mass  system  described  in  Section  3.2  the  resultant  aerodynamic  force  must,  on  symmetry  grounds, 
act  along  the  flight  path.  But  if  the  parachute  is  now  represented  as  a  solid  body,  it  need  not  possess  axial  symmetry 
and  then,  in  general,  aerodynamic  force  components  along  and  at  right  angles  to  the  flight  path  will  act  on  it. 
However,  if  m  is  made  equal  to  zero,  Fig  3.2  shows  that  the  parachute  is  modelled  as  possessing  no  mass,  then  the 
resultant  force  T,  which  it  develops  on  the  payload,  must  be  equal  and  opposite  to  the  parachute  drag  Dr  acting  along 
the  flight  path. 


To  develop  the  translational  equations  of  motion,  earth-bound  axes  are  used  as  in  Section  32,  but  for  the 
rotation.*  equation,  body  axes  will  be  adopted. 


In  Fig.33  u  body  rotates  with  angular  velocity  q  about  an  origin  O’.  Relative  to  the  earth-fixed  axes  Ox  and  Oz, 
the  velocity  components  at  the  origin  of  the  body  are  u  and  w  respectively.  Then,  relative  to  those  axes,  the  velocity 
components  at  a  general  point  P  within  the  body  are: 


I 
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a  +  qt. 


(3.9) 

(3.10) 


*p 

*p  *  w  •  q  *p - 

Relative  to  these  earth-fixed  axes,  the  component  rates  of  change  of  linear  momentum  for  the  body  are  therefore: 


=  Zm,u  +  qZmrzr  =  Zm,u 

(3.11) 

and 

=  Zm.,w  •  qXm,xf  =  £mfw  . 

(3.12) 

Whereas,  relative  to  the  body  axes  origin  O',  the  rate  of  change  of  angular  momentum  in  two  dimensions  can  be 
shown  (e.g.  Duncan”)  to  equal  1^4  where  Ia  is  tire  system  moment  of  inertia  about  the  body  axis  O'y'. 


From  Fig.3.4,  the  equations  of  motion  for  the  payload  in  three  degrees  of  freedom  are  therefore: 


mx 

s  •  FftSln  6  -  FAcos  0  *  VjpVxCoS. 

(3.13) 

mz 

=  •  FkCOS  0  •  FAsin  0  -  VjpVzCnS.  ♦  mg 

(3.14) 

and 

=  M0  -  ,/JpV,CAS,(L-L0)  sin  a  . 

(3.15) 

In  these  equations  the  drag  of  the  parachute  D,  has  been  expressed  as  1/jpVJCDS0,  while  m  and  In  respectively 
denote  the  total  mass,  inclusive  of  the  added  mass  and  the  total  moment  of  inertia,  inclusive  of  the  added  moment  of 
inertia,  for  the  system  under  acceleration.  Added  masses  and  moments  of  inertia  are  explained  later,  in  Chapter  4. 

3.3.2  The  Parachute  and  its  Payload  Modelled  as  a  Single  Rigid  Body 

The  classical  approach  to  the  derivation  of  the  equations  of  motion  for  the  trajectory  of  a  rigid  body  is  to  use  body 
axes  for  all  the  required  equations.  It  is  dearly  desirable  to  use  body  axes  for  the  rotational  equations  of  motion:  if 
they  are  also  used  for  translational  motion  the  resulting  equations  can  readily  be  linearised  to  determine  the  system 
response  to  small  disturbances  Thus  Duncan”  and  Elkin”  both  use  body  axes  in  order  to  develop  the  equations 
of  motion  for  a  rigid  body  such  as  an  aircraft,  moving  through  space. 

In  developing  the  dynamical  equations  of  motion  for  descending  parachutes,  there  are  two  problems  which  do  not 
occur  with  aircraft.  The  first  is  the  large  angles  though  which  parachutes  can  oscillate  during  their  descent  These 
might  well  limit  the  usefulness  of  any  linearisation  techniques  which  are  developed  as  a  part  of  the  solution  procedure, 
t 
( 

t 


..  F-  v 
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The  second  problem  is  the  necessity  for  the  introduction  of  added  mass  terms  in  the  parachute's  equations  of 
motion.  The  fluid  through  which  the  parachute  and  store  descend  is  real  rather  than  ideal  and  it  will  be  shown  in 
Section  4  that  it  is  necessary  to  add  experimentally-obtained  values  of  certain  added  mass  components  into  the 
equations  of  motion.  These  will  add  to  their  complexity. 

When  adopting  the  more  classical  approach  of  using  body  axes  for  both  the  translational  and  rotational  equations 
of  motion  the  fundamental  problem  is  that,  relative  to  fixed  earth-bound  axes,  the  body  axes  route.  However,  in  both 
establishing  the  equations  of  motion  and  is  the  presentation  of  their  solution,  the  motion  of  the  system  must  be 
referred  to  an  earth-fixed  axes  system. 

For  parachute  trajectory  analysis  Purvis14  has  recommended  the  use  of  body  axes  for  the  associated  routional 
equations  but  earth-fixed  axes  for  the  translational  equations,  together  with  the  necessary  axis  conversion  matrix  set 
The  translational  motion  equations  expressed  relative  to  earth-fixed  axes,  as  in  equations  3.13  and  3.14,  do  not  include 
terms  which  include  the  product  of  linear  and  angular  velocities.  Equation  3.16  which  follows,  expresses  the 
parachute's  equation  of  motion  relative  to  the  body  axis  Ox  and  in  so  doing  it  includes  the  terms  (m+c^Xrv-qw).  The 
presence  of  terms  of  this  form  adds  to  the  equation's  complexity  and  also,  when  the  parachute's  angular  velocity 
components  p.q  and  rare  large,  to  the  length  of  the  process  for  numerical  solution  of  tire  equation.  Using  earth-fixed 
axes  rather  thin  body  axes  for  the  system's  translational  motion  can  therefore  simplify  both  the  presentation  and 
solution  of  these  equations,  though  it  will  complicate  the  task  of  expressing  the  relationships  for  the  system's  attitude 
angles. 

In  the  presentation  of  the  equations  of  motion  inclusive  of  added  mass  terms,  in  accordance  with  Section  43.2, 
only  the  two  added  mass  coefficients  ku  and  k„  have  been  retained  and  these  have  been  assumed  to  be  known 
constants,  determined  experimentally  for  the  parachute  system  under  consideration. 
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Fig.3.5  EQUATIONS  OF  MOTION  FOR  A 
PARACHUTE-PAYLOAD  MODELLED  AS 
A  SINGLE  BODY 


Fbr  a  conventional  parachute  subjected  to  external  aerodynamic  and  gravitational  force  components  X,Y  and  Z  and 
external  aerodynamic  and  gravitational  moment  components  L,  M  and  N,  as  shown  in  Fig.3.5,  relative  t n  the  body 
axes  Ox,  Oy  and  Oz  the  appropriate  equations  of  motion  have  been  shown  by  Cockrell  and  Doherr34  to  be: 


X 

at 

(m  ♦  ct||)u  *{m  +  oM)(rv  •  qwMnx,(q*  ♦  r1) 

(3.16) 

Y 

* 

(m  4-  aM)(v  •  pw)  Mat  a„)ru  +  mx,(r  +  pq) 

(3.17) 

Z 

s 

(m  ♦  OmXw  ♦  pv)  -  (m  +  a„)qn  -  mx,(q  -  pr) 

(3.18) 

L 

M 

I..P 

(3.?9) 

M 

sc 

•sA  ’  «*»(*  •  qn  ♦  pv)  ♦  (I,,  -  I„)?r 

(330) 
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N  =  Ijjf  +  mx,(v  +  ni  •  pw)  -  (!„  •  I]])pq  .  (3.21) 

In  this  family  of  equation*  the  origin  O  has  been  located  at  the  canopy  centroid  and  the  store,  at  a  distance  x,  from 
O,  has  been  assumed  to  be  rigidly  connected  to  the  canopy.  The  symbols  I„  and  !„  have  been  used  to  denote 
moments  of  inertia  of  the  entire  system  about  the  axes  Ox  and  Oy  (or  Oz)  respectively  and  o„  and  eta 10  denote  the 
added  mass  components  of  the  canopy  in  the  directions  of  the  axes  Ox  and  Oy  (or  Oz)  respectively. 

(n  accordance  with  equation  4.3,  the  added  mass  components  a, ,  and  a,]  are  given  by: 

o.i  =  PfVk„  (3.22) 

and  ctjj  =  pfVkn .  (3.23) 


The  symbols  p(and  V  respectively  denote  the  density  of  the  fluid  in  which  the  parachute  is  immersed  and  the 
representative  displaced  volume  of  the  parachute  canopy.  In  comparison  with  other  terms,  the  added  masses  and 
added  moments  of  inertia  of  the  store  have  been  neglected.  The  symbols  u,v  and  w,  also  the  symbols  p.q  and  r  refer  to 
linear  and  angular  velocity  components  respectively,  along  and  about  the  axes  Ox,  Oy  and  Oz,  as  shown  in  Fig.3.5. 

In  this  presentation,  the  reasoning  of  Sedov11  and  others  has  been  followed  and  the  steady  translational  inertial 
moments,  such  as  (a,!  -  ctjj)uw,  have  been  neglected. 


Axis  of 
Symmetry 


T 

positive  direction 
of  moment  M 


Aerodynamic  forces 
developed  on  the 
payload  have 
been  neglected 


O  Origin  of  ro-ordinate 
system 

G  Parachute  system  centroid 
cp  System  centre  of  pressure 


tR 

a 

6 


N 


Resultant  velocity  at 
origin 

Angle  of  attack 
Inclination  angle  of 
parachute  system 
Resultant  aerodynamic 
fo;ce 

Normal  aerodynamic 
force  component 
Tangential  aerodynamic 
force  component 


Fig.3.6  AERODYNAMIC  AND  GRAVITATIONAL 
FORCES  AND  MOMENTS 


Negiecting  any  aerodynamic  forces  developed  on  the  payload,  for  a  parachute  which  is  oscillating  in  the  Oxz  plane, 
Fig.3.6  gives  for  the  external  forces  and  moment  about  O  acting  on  it: 


X  3 

mg  cos  0  •  T 

(3.24) 

z  = 

mg  sin  0  •  N 

(3.25) 

M  = 

N  x,  +  (mg  sin  0)  x. 

(3-24) 

and  expressions  like  these,  or  ihe  corresponding  ones  in  three  dimensions  which  are  given  by  both  Duncan”  and 
Elkin”,  should  be  inserted  into  equations  3.16  to  3.21,  the  equations  of  motion. 

When  solutions  are  required  in  all  three  planes,  Euler  angle  transformations  appropriate  both  to  the  gravitational 
forces  and  moments  acting  on  the  system  and  also  to  the  trajectory  solution  which  is  being  sought,  are  used  to  relate 
the  body  axes  to  the  earth-fixed  set  of  axis. 

The  equations  of  motion  given  in  equations  3.16  to  3.21  were  developed  from  a  model,  originally  published  by 
Tory  and  Ayres1*,  of  a  parachute  which  was  rigidly-connected  to  its  payload.  A  somewhat  similar  five-degree  of 


! 
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freedom  model  for  this  system,  but  with  rather  different  assumptions  about  the  magnitudes  of  the  unsteady 
aerodynamic  forces  and  moments,  has  been  developed  and  published  by  White  and  Wolf”. 

3.3.3  The  Parachute  and  its  Payload  Modelled  as  Two  Rigid  and  Linked  Systems 

If  the  mass  of  the  parachute  canopy  cannot  be  considered  to  be  negligible,  so  that  the  parachute  and  its  payload 
comprise  two  separate  but  linked  rigid  systems,  then  appropriate  modifications  to  equations  3.13  to  3.1S  for  the 
payload  can  be  written  and  solved  in  conjunction  with  the  equations  given  in  Fig.3.2  for  the  parachute. 

For  a  six-degree  of  freedom  coupled  payload  and  par  'chute  model,  Cutchins,  Purvis  and  Bunton110  have 
developed  the  concept  discussed  in  Section  3.3.1  of  using  eartn-fixed  axes  for  the  translational  equations  of  motion  and 
body  axes  for  the  rotational  equations  of  motion. 

Traditionally  however,  this  problem  has  been  tackled  by  the  method  developed  in  Section  3.3.2,  using  body  axes 
for  both  the  rotational  and  translational  motions.  For  example,  using  two  different  sets  of  body  axes,  Schatzle  and 
Curry”  have  developed  a  nine-degree  of  freedom  model  consisting  of  a  system  of  equations  for  a  forebody  (or 
payload)  coupled  to  a  parachute.  They  considered  the  aerodynamic  forces  and  moments  developed  on  each  body, 
together  with  the  weight  of  the  forebody.  Similarly,  allowing  for  the  weight  of  the  parachute  as  well  as  that  of  the 
forebody,  Doheir110  developed  a  body  axis  model  and  among  others.  Wolf”  has  published  a  model  for  a  coupled 
payload  and  parachute.  All  these  models  start  from  the  same  premises  but  differ  slightly  in  the  way  that  they  treat  the 
unsteady  aerodynamic  forces,  in  the  presentation  of  their  equations  of  motion  and  in  the  subsequent  linearisation 
techniques  which  they  propose. 

3.3.4  The  Parachute  as  an  Elastic  System,  Linked  to  a  Rigid  Body  Payload 

Sundberg”1  has  explained  the  application  to  trajectory  dynamics  of  finite-element  methods,  which  enable  both 
the  canopy  and  its  suspension  lines  to  be  modelled  as  flexible,  distributed  mass  structures,  coupled  to  a  rigid  payload. 
This  is  a  particularly  appropriate  model  for  die  deployment  and  inflation  phases  of  the  parachute.  An  extension  of 
this  earlier  work  has  been  made  by  Purvis1'*. 
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4.  UNSTEADY  AERODYNAMICS 

During  the  deployment  and  inflation  stages  and  even  during  descent,  much  of  a  parachute's  motion  is  unsteady. 
Whetha  the  axis  system  which  is  adopted  is  earth-bound  or  i*  fixed  within  the  canopy-store  system,  linear  an*4  angular 
components  of  acceleration,  along  and  about  these  axes,  will  occur.  In  this  chapter  the  aerodynamic  consequences  of 
this  unsteady  motion  are  discussed. 

4.1  INTRODUCTION  TO  THE  UNSTEADY  FLOW  PROBLEM 

In  describing  the  unsteady  motion  of  a  body  immersed  in  a  fluid  care  must  be  taken  in  defining  the  constituent 
pans  of  the  system  under  consideration. 

First,  consider  it  to  consist  of  a  sphere  of  mass  m  immersed  in  a  fluid.  Suppose  this  sphere  to  be  driven  through 
the  fluid  by  a  thrust  T  so  that  at  time  t  the  sphere  is  moving  with  an  instantaneous  linear  velocity  V(t)  end  an 
instantaneous  linear  acceleration  V(t).  On  the  sphere  an  aerodynamic  force  D'(t)  is  developed.  This  force  is  larger  than 
the  drag  force  D,  which  would  be  developed  if  it  were  to  move  steadily  through  the  fluid  at  velocity  V.  The  difference, 
D'-D,  certainly  depends  the  instantaneous  acceleration, V(t).  Ax  is  discussed  in  Section  4.2.2,  it  may  also  depend  on 
die  nature  of  die  accelerated  motion,  though  this  may  not  be  an  important  dependency. 

The  reason  for  this  increased  drag  is  that  the  fluid  which  surrounds  the  moving  sphere  will  also  acquire 
momentum.  Thus,  as  the  sphere  accelerates,  not  only  is  there  a  rate  of  change  in  the  momentum  of  the  sphere,  mV, 
but  in  that  of  die  surrounding  fluid.  If  this  fluid  were  ideal,  that  is  if  it  were  incompressible  and  inotatior.il.  then  the 
rate  of  change  of  fluid  momentum,  M,,  is  determinable  by  methods  to  be  explained  in  Section  7.  It  can  be  expressed 
by: 


Mr  =  D’-D  =  aV  (4.1) 

where  a  is  described  as  an  added  mass  component  for  the  system.  For  a  sphere  whose  volume  is  denoted  by  V, 
linearly  accelerating  through  an  ideal  fluid  of  density  p„  the  added  mass  component  a,  can  be  shown  to  be  equal  to 
0.5p,V. 

If  the  sphere  were  to  accelerate  through  a  real  fluid  at  distinct  from  one  which  is  ideal,  a  similar  added  mass 
component,  a„  but  of  a  different  numerical  value  from  a,,  would  be  developed.  A  real  fluid  is  one  whose  viscosity  is 
the  cause  for  it  to  possess  vorticity,  hence  to  be  rotational.  Writing  the  rate  of  change  of  momentum  for  the  body 
immersed  in  the  real  fluid  as  M*  and  assuming  a,  to  be  independent  of  time  gives: 

M»  =  T-D’ 

where  from  equation  4.1,  since  D-D  «  oV; 

then  =  T  -  D  -  oV 

and  thus  (m  +  ot)V  =  T  -  D  .  (4.2) 

where  D  is  die  dreg  force  which  the  sphere  develops  in  steady  motion.  Unlike  equation  4.1,  equation  4.2  treats  the 
immersed  body  of  mass  m  together  with  the  fluid  in  which  it  is  immersed  as  two  constituent  parts  of  a  single 
dynamical  system.  The  approach  has  a  120-year  long  history,  extending  from  both  Thomson  and  Tail*-1  and 
Kirchhoff*2.  Lamb4"  remarked  that,  “it  avoids  the  troublesome  calculation  of  the  effect  of  fluid  pressures  on  the 
suriaces  of  solids”,  which  would  be  a  necessary  procedure  if  she  aerodynamic  force  D‘  were  to  be  determined  directly  for 
the  unsteady  motion  of  a  sphere  through  an  ideal  fluid.  Both  the  importance  and  the  method  of  determining  added 
mass  components  relevant  to  tbe  motion  of  a  parachute  canopy  through  ideal  fluids  remains  to  be  discussed  in  Section 
7.1.  No  analytical  methods  exist  by  which  aerodynamic  forces  can  be  determined  for  unsteadily-moving  bodies 
immersed  in  real  fluids:  the  ideal  fluid  concept  gives  an  approximate  model  which,  for  a  Muff  body  like  a  conventional 
parachute,  is  of  uncertain  value. 

In  exactly  the  same  way  as  that  outlined  in  equation  4.1,  when  any  immersed  body  moves  unsteadily  in  any 
direction  through  a  fluid  the  aerodynamic  forces  and  moments  are  developed  on  it  aiger  in  magnitude  from  their  steady 
state  values.  In  a  given  problem,  whether  or  not  this  difference  is  of  any  enginesing  significance  depends  on  the 
relative  magnitudes  of  the  rue  of  change  of  momentum  for  the  immersed  body,  in  U*M  of  the  fluid  in  which  it  is 
immersed,  M*.  This  ratio  is  a  function  of: 

(0.  the  immersed  body  density  compared  with  that  of  the  fluid  which  it  displaces; 

(ii) .  the  immersed  body  shape  inclusive  of  its  porosity,  if  any;  and 

(iii) .  the  direction  in  which  these  rues  of  change  of  momenta  occur. 
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For  parachutes,  the  ratio  is  certainly  of  engineering  significance  when  considering  parachute  inflation  and  also 
when  modelling  the  oscillatory  motion  of  parachute  canopies  in  their  descent  through  the  air. 

4.2  THE  ADDED  MASS  CONCEPT  IN  PARACHUTE  UNSTEADY  MOTION 
4.2.1  General  Considerations 

There  are  two  important  applications  of  the  added  mass  concept  to  unsteady  parachute  motion.  The  first  is  to 
canopy  inflation.  Knacke14  states  that  Scheutel^/epotted  in  1946,  was  the  first  to  investigate  added  mass  effects  on 
parachute  opening  processes  and  these  will  be  discussed  here  in  Sections  43.2  and  S3.  The  second  application  is  to 
the  parachute’s  dynamic  stability  characteristics  in  pitch  during  its  fully-inflated  descent  motion.  Whereas  static 
stability  has  been  shown  in  Section  2.33  to  be  concerned  with  the  direction  of  the  moments  developed  on  a  parachute 
system  after  it  has  been  disturbed  from  the  equilibrium  state,  dynamic  stability  characteristics  are  related  to  the 
frequency  and  to  the  damping  of  the  oscillatory  unsteady  motions  which  then  ensue.  To  analyse  these  motions 
satisfactorily  it  is  necessary  to  write  and  to  solve  systems  of  equations  similar  in  form  to  those  expressed  in  Section 
33.  For  the  adequate  analysis  of  both  the  canopy  inflation  phase  and  the  oscillatory  motion  of  the  parachute  system 
what  the  canopy  is  fully  inflated,  it  is  first  necessary  to  determine  the  added  mass  components  appropriate  to  a  given 
canopy  shape.  If  analytical  methods  are  to  be  adopted,  it  is  necessary  to  assume  that  the  fluid  field  in  which  the 
parachute  canopy  is  immersed  is  ideal 

With  the  application  to  parachute  canopies  in  mind  and  following  Lamb  and  Milne-Thomson*-5  among  other 
authors,  Ibrahim44  in  1965  evaluated  some  of  the  added  mass  components  for  cup-shaped  parachute  canopies  in  ideal 
fluid  flow.  Later,  Klimas4-1  demonstrated  that  these  analytical  values  are  strongly  dependent  on  the  porosity  of  the 
canopy.  By  methods  such  as  they  adopted  it  can  be  shown  that  the  canopy  added  mass  component  a  is  a  function  of 
iu  shape  and  attitude.  It  is  also  proponioral  to  the  body  size  and  the  fluid  density,  p,  thus: 

at,  *  k,pfV,  (43) 

where,  in  ideal  fluid  flow,  the  added  mass  coefficient  k,  is  a  constant  for  a  given  immersed  body  shape  and 
attitude.  The  symbol  V  is  used  fer  the  representative  displaced  volume.  Where  a  straightforward  definition  is 
possible,  e.g.  for  an  immersed  sphere,  the  representative  displaced  volume  is  the  actual  displaced  volume.  For  a 
parachute  canopy  it  is  conventionally  considered  to  be  that  of  a  hemisphere  (though  some  authorities  define  the  added 
mass  coefficients  for  a  canopy  in  terns  of  a  sphere  rathe  than  a  hemisphere)  having  a  diameter  D.  equal  to  the 
nominal  diameter  of  the  canopy,  Lr 
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Some  appropriate  analytical  values  of  added  mass  coefficients  are  given  in  Table  4.1.  In  this  Table  the  added 
moment  of  inertia  coefficient  kg,  measured  about  an  axis  at  right  angles  to  the  plane  containing  k„  and  k».  is  defined 
in  terms  of  the  added  moment  of  inertia  component  as 

*  *%At  (4.5) 

where  I,  is  the  moment  of  inertia  of  the  fluid  in  the  representative  volume.  For  a  parachute  canopy: 

t,  *  <pp:wi6  .  «.6) 

the  density  of  the  fluid  in  which  the  canopy  is  immersed  being  denoted  by  Pc 

A  number  of  experiments  have  been  performed  to  determine  the  added  mass  components  for  bluff  bodies  moving 
unsteadily  through  funds  since  at  ksact  as  long  ago  as  1836,  when  Bessel”  tested  the  periodic  motion  of  a  spherical 
pendulum,  both  in  air  and  under  water.  As  shown  in  Table  4&  many  of  these  results  have  proved  to  be  inconclusive. 
Often  the  added  mass  components  which  have  been  evaluated  differed  only  marginally  from  those  obtained  analytically 
for  unsteady  motion  through  an  ideal  fluid.  Howeter,  other  tests  have  shown  some  marked  differences. 
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Table  4 2  ADDED  MASS  COEFFICIENTS  DETERMINED  EXPERIMENTALLY 
FOR  BLUFF  BODIES  IN  REAL  FLUIDS 

4.2,2  Determining  the  Added  Mass  Components 

In  their  determination  of  the  added  mass  components  from  suitably-designed  experiments  most  workers  adopted 
methods  similar  to  those  which  Iversen  and  Balew*4'  deserved  in  1951.  From  equation  42; 

(m  *  o,)V  *  T  -  D 

and  writing  the  steady-state  drag  D  as  equal  to  '/jpV^S.Co,  they  considered  the  appropriate  component  of 
aerodynamic  force,  F(t)  to  consist  of  two  parts:  oae  which  could  be  expressed  in  terms  of  the  instantaneous  velocity, 
V(t)  and  the  other  which  could  be  expressed  in  terms  of  the  instantaneous  acceleration,  V(t).  Thus  in  unsteady  flow 
the  tool  aerodynamic  force  F(t)  at  time  t  was  written  in  terms  of  coefficients  a  and  b  as: 

F(t)  -  aV*(t)  ♦  bV(t)  .  (4.7) 
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It  remains  to  devise  an  appropriate  experimental  programme  from  which  10  determine  a  and  b. 

One  way  of  doing  this  is  would  be  to  measure  the  aerodynamic  farce,  F,  at  a  steady  velocity  V  fust,  so  obtaining 
the  component  a.  The  component  b  would  then  be  given  from  the  difference  in  magnitude  between  the  aerodynamic 
forces  developed  on  the  immersed  body  in  corresponding  unsteady  and  steady  motion,  F(t)-Fy  Uring  equation  4.3,  the 
appropriate  added  mass  coefficient  could  then  be  calculated  from  the  component  b.  This  method  assumes  that  the 
velocity-dependent  component  of  aerodynamic  force  in  unsteady  motion  is  unchanged  from  its  corresponding  value  in 
steady  motion,  an  assumption  which  is  difficult  to  justify  on  physical  grounds. 

In  a  method  applicable  to  periodic  motion,  which  owes  much  to  both  Smpkaya  &  Isaacson41’  and  also  to  Bcarman 
et  al*31,  avenge  values  of  the  components  a  and  b  are  calculated  by  first  multiplying  equation  4.7  either  by  V(t)  or  by 
V(t),  then  taking  mean  values  of  the  resulting  relationship  over  a  period  T ,  so  that 


»UTF(t)V(t)  dt 

«  VI  aV»(t)  dt  *  VI  bV(t)V(t)  dt 

(4.8) 

and  VT/,TF(t)V(t)  dt 

*  VI  aV*(t)V(t)  dt  +  %!I  bVHt)  dt  . 

(4.9) 

Avenge  values  of  a  and  b  can  then  be  determined,  since 

Si  V(t)V(t)  tit 

a  SI  V\t)V(t)  dt  *  0  . 

(4.10) 

A  simpler  procedure  would  be  to  apply  the  method  described  in  Section  2.4  by  first  writing  equation  4.7  in  the 
form: 


{F(t)/V‘(t)>  =  a  +  b(V(t)/VJ(t)}  .  (4.11) 

This  is  » linear  relationship  between  the  two  variables  {F(t)Af2(t))  and  ( V(tyV*(t)}.  Having  plotted  experimental 
data  in  this  fotm,  regression  analysts  can  be  used  to  yield  values  for  tbe  two  components  a  and  b. 

An  advantage  of  this  latter  technique  over  the  one  described  previously  is  that  it  is  immediately  possible  to  see 
from  the  graph  of  equation  4.1 1  if  the  two  components  a  and  b  are  sensibly  constant  over  an  entire  period,  or  if  they 
vary.  For  example,  they  could  possess  different  values  when  the  acceleration  of  the  system  is  positive  from  those 
when  it  is  negative.  An  illustration  by  Harwood0*,  in  the  course  of  his  determination  of  Table  4 3,  is  shown  in  Fig. 


4.1.  This  graph  indicates  that  although  equation  4.7  is  an  approximation  it  can  serve  as  an  acceptable  engineering 
relationship  for  the  unsteady  aerodynamic  force  F(c)  developed  ou  parachute  canopies. 

Experimental  conditions  under  which  the  various  testa  reported  in  Table  42  were  performed  varied  greatly.  In 
some,  the  unsteady  motion  was  unidirectional.  In  others  it  was  either  wholly  oscillatory  or  a  combination  of  these 
two  modes.  Methods  of  dctenninaiion  of  the  added  mass  coefficient  varied  too;  in  some  experiments  average  values  of 
aerodynamic  force  were  determined,  in  others  instantaneous  values  were  obtained.  In  most  of  these  tests,  the 
experimental  uncertainty  with  which  the  added  mass  coefficients  were  obtained  was  high. 

Among  others,  Hamilton  and  Lindeli02  have  suggested  that  the  manner  in  which  unsteady  motion  is  imparted  to 
an  immersed  body  is  significant  and  there  is  a  significant  weight  of  experimental  evidence  to  support  their  view. 
They  suggested  that  only  if  there  are  long  periods  of  constant  velocity  or  constant  acceleration  can  the  instantaneous 


aerodynamic,  force  F(0  be  well  aprened  solely  in  terms  of  the  instantaneous  velocity,  V(t)  and  tki  incise  taneous 
accetofadon,V(t),  in  the  manner  proposed  in  equation  4.7.  Hamilton  and  LiadsU  stated  that  ‘the  practice  of  expressing 
Add  force  as  a  function  of  the  instantaneous  velocity  and  the  acceleration  of  the  body  (or  the  fluk?)  has  carried  over  to 
engineering  problems  in  which  convective  accelerations,  flow  separation  and  wakes  are  important  (nonlinear)  cases. 
The  force  is  dctem.ioed  experimentally  and  iu  relationship  to  velocity,  acceleration  snd  other  parameters  is  expressed 
in  a  variety  of  ways.  The  equations  may,  or  may  not,  include  an  added  mass  terra  and,  if  they  do,  it  may  be  a  constant 
or  a  function  of  various  perimeters.  The  choice  depends  oo  the  kind  of  motion  and,  to  some  extent,  on  the  view  of 
the  author1.  Instead  they  proposed  that  a  more  general  form  of  equation  4.7  should  be  written: 

F(t)  *  aVS't)  ♦  bV(t)  ♦  Hew  history  term  (4.11a) 

the  form  of  the  latter  term  depending  on  the  manner  it,  which  the  unsteady  motion  is  imparted  to  the  immersed 
body. 

43  SIGNIFICANCE  OF  ADDED  MASS  COEFFICIENTS  TO  PARACHUTE  UNSTEADY 
MOTION  PREDICTION 

4.3.1  Their  Historical  Importance 

For  fluid  dynamical  problems  concerning  parachute  canopies,  whose  drag-producing  capabilities  are  among  their 
most  significant  aerodynamic  characteristics,  to  base  analysis  oo  the  assumption  that  the  fluid  through  which  the 
canopy  descends  is  ideal  would  appear  to  be  unrealistic,  for  in  such  a  fluid  no  separated  wake  could  be  developed  and 
steady  motion  descent  would  necessarily  be  drag  flee.  It  is  thus  ctear  that  results  obtained  from  such  a  mathematical 
model  must  be  validated  through  appeal  to  appropriate  experimental  programmes. 

Considering  the  canopy  deployment  and  inflation  application  and  based  on  Hot,  visualisation  studies  for  inflating 
parachute  canopies,  LinganP”  and  others  have  suggested  that  since  cancpy  deployment  is  a  very  rapid  process  it 
may  well  be  realistic  to  consider  the  surrounding  flow  field  as  though  it  were  iirocarional  and  thus,  if  It  were  also 
effectively  incompressible,  as  ideal.  0'Hara°\  for  eumpls,  assumed  that  the  added  mass  coefficient  for  an  inflating 
parachute  canopy  was  equal  to  that  for  a  flat  disc  immersed  in  an  ideal  fluid,  with  a  diameter  equal  to  that  of  the 
canopy  and  thus  during  inflation  its  added  mass  would  increase  as  the  cube  of  its  diameter. 

In  his  1944  p^ter  on  parachute  descent  behaviour,  Henh1*  discussed  the  significance  of  added  mass  components  on 
parachute  dynamic  stability.  He  called  these  components  ‘co-acceierated  air  masses'  and  sued  that  they  wete  of  the 
order  of  size  of  the  combined  mass  of  the  canopy  and  itt  load.  For  their  determination  in  two  degrees  of  freedom  he 
introduced  an  ellipsoid  having  the  same  major  diameter,  volume  and  location  as  the  parachute  canopy,  using  the 
analytical  values  for  this  ellipsoid  of  the  axial,  transverse  and  rotational  added  mass  coefficients  which  Lamb43  bad 
previously  determined  and  which  have  been  given  in  Table  4.1.  He  concluded  that  added  masses  have  a  significant 
effect  oo  the  deteminaffoa  of  parachute  dynamic  stability.  Larr,  Lester434  reformulated  Henn’s  equations.  Both  Henn 
and  be  had  shown  that  added  mass  coefficients  were  of  importance  in  determining  parachute  dynamic  stability 
characteristics.  Lester  commented  on  the  unsatisfactory  procedure  of  using  analytical  values,  derived  from  the 
behaviour  of  ideal  fluids  around  bluff  parachute  canopies,  for  these  coefficients. 

White  and  Wolfs’-7  1968  paper  on  parachute  dynamic  stability  and  Wolfs’*  later  1971  contribution,  while 
recognising  that  the  added  mass  coefficients  were  tensors,  nevertheless  over-simplified  the  problem  of  representing 
parachute  unsteady  motion.  More  recently  Eaton4*  has  reformulated  the  analytical  problem,  discussing  the  relative 
significance  of  the  analytical  values  far  the  added  mass  coefficients  which  various  authors  have  obtained.  However,  he 
presented  its  solution  without  recourse  to  experimentally -determined  added  mass  coefficients. 

In  1965  Ibrahim447  published  a  paper  describing  an  experimental  method  of  determining  added  moments  of  inertia 
for  parachute  canopies.  Apart  from  this  earlier  work,  systematic  experiments  to  determine  added  mass  coefficients  for 
parachute  canopies  were  not  reported  until  Yavuz4*  first  published  his  work  in  1982. 

4.3.2  Their  Contemporary  Importance 

During  canopy  inflation  there  is  no  doubt  that  the  effects  of  added  mass  are  of  significance.  If  they  were  neglected 
and  instead  a  series  of  steady  flow  solutions  were  obtained  for  canopies  with  increasing  degrees  of  inflation  in 
potential  flow,  it  is  highly  unlikely  that  good  qtproximations  would  be  obtained  for  pressure  distributions  within  the 
inflating  canopies.  In  the  continuity  relationship  given  as  equation  5.10  in  Section  5.3.2  Heinrich  introduced  the 
added  mats  of  the  inflating  canopy.  However,  as  the  aerodynamics  of  canopy  inflation  is  still  very  much  in  its 
infancy,  even  very  approximate  estimates  for  the  pressure  distribution  round  the  inflating  canopy  would  be  of  real 
value.  Insufficient  material  bas  been  published  to  warrant  any  further  discussion  here  of  canopy  added  manes  during 
the  inflation  process,  in  consequence  only  the  dynamic  stability  pitching  characteristics  of  fully-inflated  canopies  will 

For  a  body  which  moves  unsteadily  through  a  fluid  the  added  mass  coefficients  form  a  second-order  tensor  with 
twenty-one  independent  components,  comprising  six  in  which  both  the  force  and  the  acceleration  components  me 
tnasbaory,  six  in  which  both  moments  and  angular  acceleration  components  are  rotary  and  nine  which  describe  mixed 
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uanslatory  and  rotary  behaviour.  Following  the  form  of  the  equations  derived  in  Section  3.4  for  a  parachute  system 
with  six  degrees  of  freedom  and  after  making  a  number  of  simplifying  assumptions,  Cockrell  and  Doherr1'*  argued  that 
for  a  conventional  parachute  canopy  considered  to  be  rigidly  connected  to  its  payload  and  for  which,  as  shown  in 
Fig.2. 1,  Ox  is  the  axis  of  symmetry  with  planes  of  symmetry  Oxy  and  Oxz,  only  four  indej.  -ndent  and  significant 
added  mass  coefficients  need  be  considered.  In  the  tensor  nomenclature  which  they  adopted  the  first  subscript  denoted 
the  direction  in  which  the  unsteady  force  was  measured  and  the  second  (hat  of  the  acceleration  causing  the  added  mass 
component  under  consideration.  The  numbers  1,2  and  3  implied  linear  motion  in  the  directions  of  the  axes  Ox,  Oy 
and  Oz  respectively  while  4,5  and  6  described  angular  motion  respectively  about  the  axes  Ox,  Oy  and  Oz.  Thus  k„  is 
the  added  mass  coefficient  (referred  to  as  the  axial  added  mass  coeffickwu)  which  is  determined  when  an  x-dirccted  force 
is  measured  on  a  canopy  undergoing  an  x-directed  linear  acceleration,  Le  along  the  axis  of  symmetry. 

There  four  significant  added  mass  coefficients  are: 


kn 

= 

a/oV ; 

(4.12) 

(is). 

k  „ 

= 

kji  ;. 

(4.13) 

There  latter  are,  respectively,  the  added  mass  coefficients  determined  when  a  y-directed  force  is  measured  on  a 
canopy  undergoing  a  y-directed  linear  acceleration  and  when  a  z-directed  farce  is  measured  on  a  canopy  undergoing  a  z- 
directed  linear  acceleration.  They  are  non-dimensionalised  from  their  appropriate  added  mass  components  (Xg  and  a,3 
in  exactly  the  same  way  as  was  k„  in  equation  4.3; 

(iii).  k„  =  k„  (4.14) 

These  are  the  added  moment  of  inertia  confidents  determined  when  a  moment  about  the  y-axis  is  measured  on  a 
canopy  undergoing  mi  angular  acceleration  about  the  y-axis  and  when  a  moment  about  the  z-axis  is  measured  on  a 
canopy  undergoing  an  angular  acceleration  about  the  z-axis,  respectively.  As  has  been  shown  in  equations  4.4, 4.5 
and  4.6,  they  are  non-dimeasionalired  from  their  respective  added  moment  of  inertia  components,  otj,  and  a*  thus: 

kM  =  all/(xplD’/192)  (4.15) 

and  k«  a46/(7tpID’/192)  .  (4.16) 


(iv).  The  remaining  significant  added  mass  coefficient  is: 

k*  =  kM,  (4.17) 

non-dimensionaUsed  from  a**  Ojj.  This  is  the  added  mass  coefficient  which  is  determined  when  a  y- 
directed  force  is  measured  on  a  canopy  undergoing  an  angular  acceleration  about  the  O-z  axis,  or  when  a  z-directed  force 
is  measured  on  a  canopy  undergoing  aa  angular  acceleration  about  the  O-y  axis.  Yavuz*-3*  has  shown  that  if  the 
origin  of  the  co-ordinate  system  is  located  close  to  the  canopy  centre  of  pressure  then  this  latter  coefficient  is  of 
negligible  magnitude. 

Hence  the  problem  of  determining  added  mass  coefficients  has  become  one  of  determining  experimentally  only 
three  added  mass  coefficients,  ku,ks=kuandks:sktt.  These  three  can  be  further  reduced  to  two  coefficients.  If  the 
origin  of  the  co-ordinate  system  is  locatea  dose  to  the  canopy  centre  of  pressure  then  the  apparent  moment  of  inertias 
of  the  canopy,  Os,  about  the  axis  O-y  and  a*  about  the  axis  O-z  will  be  totally  dominated  by  the  moments  of  inertia 
of  the  payload  and  hence  they  can  be  neglected. 

The  two  remaining  added  mass  coefficients  are  kn  and  kg  =  k„.  In  the  determination  of  parachute  dynamic 
stability  characteristics  these  cam  be  shown  to  be  significant. 

in  tests  performed  on  parachute  canopies  which  moved  steadily  under  water  while  they  were  forced  to  oscillate  in 
either  their  axial  or  transverse  direction,  Cockrell,  Shea,  Harwood  and  Baxter*-"  obtained  the  average  values  for  the 
added  mass  coeficients  k„  and  k»  which  are  given  in  columns  5  and  6  of  Table  4.3.  From  these  experiments  it  is 
evident  that  the  real  fluid  flow  values  obtained  for  them  added  mass  coefficients  substantially  exceed  the  potential  flow 
evaluations  which  Lamb,  Ibrahim  and  others  earlier  determined. 

Using  Table  4.1  and  idealising  a  parachute  canopy  into  an  ellipsoid  having  a  length/diameter  ratio  of  0.5, 
analytical  values  for  the  added  mass  coefficients  k,,  *  0.70  and  kn  «  0.21  can  be  obtained.  These  ate  based  on  the 
ellipsoid's  displaced  volume,  which  is  equal  in  magn  itude  to  that  of  a  hemisphere  whose  diameter  is  equal  to  length  of 
the  ellipsoid’s  major  axis.  Hence,  for  this  model  for  the  parachute  canopy,  analytical  values  based  on  the  canopy 
projected  diameter  are  abo  0.70  and  02 1  respectively.  Assuming  the  projected  diameter  to  equal  0.7  of  the  nominal 
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projected  diameter  are  also  0.70  and  0.21  respectively.  Assuming  the  projected  diameter  to  equal  0.7  of  the  nominal 
diameter  the  analytical  values,  in  terms  of  the  representative  displaced  volume  based  on  nominal  canopy  diameter  D„ 
would  equal: 


ku(analytical)  =  0.24;  k, , (analytical)  =  0.07.  (4.18) 

Corresponding  values  which  Doherr  and  Salieris00  used,  written  in  terms  of  the  nominal  canopy  diameter  D.  and 
the  present  axis  convention  were: 


ku(analytical)  =  034:  k„(analytlcal)  =  0.17.  (4.19) 


CANOPY 

TYPE 

CT 

Steady 

State 

■cr 

Un- 

-steady 

Motion 

Steady 

State 

/radian 

Ct  -a  0 
Unsteady 
Motion 
/radian 

Added  Mass  Coefficients 

*11 

Axial 

*33 

Transverse 

Round  Canopies: 

Without  Drive  Slots 

0.63 

0.64 

-0.40 

- 

1.1 

02 

With  Drive  Slots 

0.61 

0JS7 

-0.26 

- 

1.0 

Cross  Canopies: 

3:1  Arm  Ratio 

Imporous 

0.64 

0.68 

-0.19 

-0.13 

14 

03 

Very  PorousJU23" 

0.66 

0.67 

+0.46 

+0.41 

0.8 

0.1 

_  4:1  Arm  Ratio 

Imporous 

0.77 

0.76 

+0.52 

+0.43 

1.6 

0.1 

cubic  ft/sq.ft/sec.  measured  at  10  inches  of  water  pressure 
Table  43  EXPERIMENTAL  CHARACTERISTICS  FOR  PARACHUTE 
CANOPIES  IN  STEADY  AND  IN  UNSTEADY  MOTION 
(after  Cockrell,  Shen,  Harwood  and  Baxter'-”  and  later  revised  by  Harwood'-”) 

Ibrahim",  modelling  the  canopy  as  if  it  were  a  spherical  cup,  obtained  a  value  of  k„  (based  on  projected  diameter) 
of  about  2.6.  In  terms  of  nominal  diameter  D.  this  becomes  a  value  of  about  three  times  that  in  equation  4.19;  of  the 
order  of  the  experimental  values  shown  in  Table  43.  They  vary  with  the  volume  of  air  enclosed  by  the  canopies  and 
they  decrease  appropriately  with  increasing  canopy  porosity.  The  uncertainties  in  the  measurements  of  the  transverse 
added  mass  coefficients  kg  are  high  but  they  are  seen  to  be  of  the  order  of  one  fifth  to  one  tenth  the  corresponding 
values  of  k„. 


Fig.  4.2  EFFECT  OF  ADDED  MASS 
COEFFICIENTS  ON  CANOPY 
DYNAMIC  STABILITY 
after  Doberr  &  Saliaris,  rtf.4J8 
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The  vdodty-dependent  tangential  force  coefficients  measured  in  these  tests  and  presented  in  cotumn  2  of  Table  43 
arc  approximately  equal  the  corresponding  tanperuisl  force  coefficient:  in  steady  motion  given  in  column  1.  Values  of 
(dC* /da  determined  in  unsteady  mod  .and  presented  in  column  4  of  the  table  do  not  differ  appreciably  from 


ssitSL3’2  “  *"*>  ■"*iwi"«  <»»■* 

«s dwwn m Section 13.2.  is  tbe  same condition as  that  for Sk sSjhTnS,’ mustbe!ar8eand Positive.  TOs, 

'^ue  ^  ^  of  pitching  oscillations  ov^tL  the/SdT^  of  "»tion  «  to 

ff^  °"tbe  o^olatwo  damping  rate.  A  typical  roots-locus  diagran  from  .efa^^'  Tth  Uuk  corresponding 
In  the  curve  a  shown  in  this  figure  both  the  axial  and  the 4-29  ,sshown  foRg.4.1 
By  considering  the  added  mas,  coefficients  k„  SSL  m fddedtmass  coefficients  have  been  neglected, 
characiehstic  there  is  a  greater  fierpiency  of  oscillation  for  a  gfren  valued  dC^T'  bhhas.been  opined.  With  this 
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5 .PARACHUTE  DEPLOYMENT  AND  INFLATION 

Since  parachutes  must  be  designed  to  be  sufficiently  strong  to  withstand  the  opening  loads  which  occur  as  part  of 
the  inflation  process,  when  considering  parachute  aerodynamics  ah  appreciation  of  the  parameters  which  influence  these 
loads  and  the  manner  in  which  they  relate  to  these  independent  parameters  is  of  fundamental  concern.  Because  inflation 
is  an  unsteady  aerodynamic  phenomenon  in  which  there  are  necessary  changes  in  the  canopy  shape  during  the  process, 
it  is  no  simple  matter  to  develop  an  adequate  inflation  model.  Aerodynamic  models  which  are  available  range  from 
rcadily-ustble  empirical  methods  to  models  which  by  comparison  are  so  complex  as  to  be  almost  useless  to  the 
parachute  designer. 

5.1  INTERACTION  FORCE  BETWEEN  CANOPY  AND  PAYLOAD 

During  the  canopy  deployment  phase,  the  time-varying  interactive  force  developed  between  the  canopy  and  its 
payload  significantly  exceeds  its  steady-state  value  and  it  is  therefore  of  importance.  This  interactive  force  can  be 
measured  and  interpreted  either  when  the  parachute-payload  system  is  in  free  flight  or  when  it  is  in  a  wind  tunnel. 
However,  it  is  very  difficult  to  establish  a  mathematical  model  which  contains  all  the  relationships  necessary  for  it  to 
be  adequate  for  design  purposes.  From  such  a  model  the  maximum  loading  during  deployment  and  inflation  is 
certainly  required.  Some  information  is  desirable  so  that  the  trajectory  of  th;  parachute  and  payload  during  the  inflation 
process  can  also  be  predicted.  Rom  more  sophisticated  models  it  might  be  possible  to  predict  with  confidence  the 
stress  distribution  over  the  canopy.  It  certainly  would  be  desirable  if,  during  the  inflation  process,  the  model  were 
capable  of  predicting  the  canopy  shape  and  size. 

Methods  of  analysis  have  been  developed  in  three  different  directions: 

(i) .  the  determination  of  physically-based  relationship-'  the  closure  of  which  depends  on  the  development  of 

appropriate  empirical  expressions; 

(ii) .  the  reproduction  of  dynamic  similitude  in  either  fufl-xale  or  model  test  situations  through  the  establishment 

of  appropriate  functional  relationships  for  the  required  unknown  parameters;  and 

(iii) .  the  construction  of  complex  mathematical  mode’a  which  depend  on  knowledge  of  the  pressure  distribution 

around  the  inflating  canopy  in  order  ur  establish  the  force  required  to  drive  the  inflation  process. 

5.1.1  Expression  for  the  Interactive  Forces,  F, 

In  order  to  develop  expressions  for  the  time-varying  force  which  is  developed  between  the  canopy  and  its  payload, 
consider  in  the  direction  of  the  flight  path  the  separate  equations  of  motion  for  the  parachute  canopy  and  the  payload,  as 
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are  shown  in  Fig  5.1. 

For  the  parachute  canopy,  neglecting  its  component  of  weight  down  the  flight  path  in  comparison  with  the  canopy 
drag  D,  and  the  interactive  force  F,: 

F,  -  D,4.m,V*A,V  (5.1) 

As  Lingard**  has  indicated,  the  total  interactive  force  F,  depends  on  both  the  aerodynamic  force  D,  and  the 
inertial  force  -<m,V  +  m,V).  When  establishing  the  inertial  force  during  the  deployment  and  inflation  phases. 
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developing  an  expression  for  the  canopy  mass,  inclusive  of  its  added  mass  (the  latter  defined  in  Section  4.1)  as  a 
function  of  time  presents  real  difficulties.  The  aerodynamic  force  could  be  measured  directly  in  an  appropriate  wind 
tunnel  test,  as  could  the  inertial  force,  e.g.  Lingard54,  but  if  the  canopy  is  rigidly  mounted  in  the  wind  tunnel  the 
measured  interactive  force  would  not  include  the  inertial  term  in  equation  S.l. 

In  figure  5.2a  conditions  in  such  a  wind  tunnel  test  are  shown.  Under  these  circumstances,  not  only  is  the  inertia 
force  excluded  but  as  the  canopy  inflates  the  velocity  of  air  relative  to  it  does  not  decrease  but  will  remain  of  constant 
magnitude.  The  measurement  of  the  interactive  force  is  described  as  having  been  made  under  an  irftnite  mass 
condition,  since  if  the  canopy  mass  were  infinite  then  under  the  action  of  a  finite  force  there  could  be  no  deceleration  of 
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Fig  5.2  Forces  Developed  Between  the  Canopy  and  the  Payload 

During  Inflation 
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the  relative  airflow. 

By  contrast,  in  figure  5.2b  canopy  inflation  is  shown  as  taking  place  under  free  flight  conditions.  Under  such 
circumstances  deceleration  of  the  relative  airflow  occurs  and  there  is  a  consequential  reduction  in  the  interactive  forces. 
Appropriate  measurements  are  now  described  as  having  been  made  under finite  mass  conditions,  when  the  peak  opening 
force  F„  occurs  considerably  earlier  than  at  the  instant  when  the  canopy  becomes  fully  inflated. 

The  processes  of  canopy  deployment  and  inflation  can  be  divided  into  a  number  of  distinct  stages.  The  first  of 
these  is  before  deployment,  when  no  aerodynamic  force  is  developed  on  the  canopy.  The  trajectory  of  the  canopy  and 
the  payload  system  is  then  solely  determined  by  the  system's  initial  conditions  and  by  the  payload's  aerodynamics  and 
weight.  Next  follows  the  initial  deployment  stage,  which  ends  when  the  canopy  rigging  lines  are  fully  stretched.  As 
this  occurs  the  interactive  force  peaks  to  a  local  maximum,  called  the  lines -taut  snatch  force.  Then  follows  the 
inflation  stage:  this  is  from  the  instant  at  which  the  parachute  canopy  begins  to  open  until  it  first  reaches  its  normal 
fully-inflated  projected  area.  At  the  commencement  of  this  third  stage  the  canopy  skirt  forms  a  mouth  which  begins  to 
open  and  inhale  the  air.  Within  the  canopy  this  inhaled  air  forms  a  bail  which  moves  down  the  length  of  the  partly- 
inflated  canopy  to  the  vent,  with  which  it  impacts.  Once  this  impact  has  occurred  the  canopy  begins  to  inflate  radially 
and  during  this  process  the  interactive  force  achieves  a  peak  opening  load  As  figure  5.2  indicates,  this  peak  opening 
load  is  the  maximum  interactive  force  which  me  system  experiences. 

In  this  figure  the  interactive  forces  are  shown  as  functions  of  time,  the  inflation  time  t,  being  that  between  the 
peak  of  the  lines-laut  snatch  force  and  the  peak  value  of  the  inflation  force  or,  mote  strictly,  between  the  peak  values  of 
their  respective  coefficients.  Though  the  inflation  time  is  shorter  than  the filling  time  if,  defined  as  the  time  taken 
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production  of  the  snatch  force  to  the  instant  at  which  the  canopy  first  reaches  its  steady-state  diameter  and  is  thus 
considered  to  be  fully  inflated,  in  an  experimental  situation  the  inflation  time  is  more  clearly  defined  than  is  the  filling 
time. 

5.2  CRITICAL  OPENING  AND  CLOSING  SPEEDS.  SQUIDDING 

Consider  an  uniuflated  canopy  in  a  wind  tunnel.  If  the  relative  velocity  of  the  airflow  is  increased  a  speed  will  be 
reached  at  which  the  canopy  will  just  inflate.  The  practical  applications  of  parachutes  ensure  their  usage  at  much 
higher  relative  air  speeds  than  this  minimum  inflation  speed  so  this  latter  is  not  considered  to  be  one  of  the  canopy's 
critical  speeds. 

If  in  the  wind  tunnel  the  canopy  were  now  fully  inflated  at  a  relative  airflow  which  is  above  this  minimum 
inflation  speed  and  if  this  relative  velocity  were  now  gradually  increased,  conditions  would  eventually  occur  at  which 
the  canopy  collapses  to  a  form  in  which  its  maximum  diameter  was  only  between  one-quarter  and  one-third  of  its  fully- 
inflated  diameter.  Sirvte  the  shape  of  the  canopy  then  resembles  that  of  a  sqid,  this  phenomenon  of  canopy  partial 
closure,  occurring  at  a  critical  closing  speed ,  is  termed  squidding.  The  occurrence  of  this  squid  configuration  is  a 
consequence  of  premature  equilibrium  between  the  radial  pressure  load  and  structural  tension,  thus  the  critical  closing 
speed  varies  from  one  canopy  shape  to  another  and  is  a  function  of  the  canopy  porosity.  Still  further  increases  in  the 
relative  air  speed  do  not  cause  the  squidding  parachute  shape  to  alter  greatly.  When  the  relative  air  speed  is  reduced 
below  this  critical  closing  speed  the  crown  of  the  squidding  parachute  begins  to  inflate  and  at  some  considerably  lower 
relative  air  speed,  called  the  critical  opening  speed,  the  canopy  suddenly  opens  fully. 

Since  it  is  the  maximum  speed  at  which  the  canopy  fully  inflates  the  critical  opening  speed  is  of  importance  and 
the  velocity  at  which  a  canopy  is  required  to  deploy  must  be  less  than  the  critical  opening  speed.  Once  the  canopy  has 
fully  inflated,  the  air  speed  relative  to  it  seldom  increases  any  more  and  thus  the  critical  closing  speed  is  of  much  less 
practical  significance. 

5.3  CANOPY  INFLATION  THEORIES 

In  his  1927  examination  of  canopy  inflation  physics,  M Oiler31  applied  the  principle  of  conservation  of  mass  to  the 
control  volume  defined  by  the  physical  boundary  of  the  parachute  canopy.  In  so  doing,  he  stated  that  the  rate  of 
increase  of  the  canopy  volume  was  equal  to  die  product  of  the  canopy  mouth  area  and  the  canopy  speed,  In  filling  time 
models  which  developed  from  this  work  empirical  expressions  were  established  for  the  variation  of  the  canopy  mouth 
area  with  time  during  inflation.  Through  his  proposal  that  "the  distance  necessary  for  the  complete  inflation  of  a  given 
canopy  is  a  constant  and  is  proportional  to  the  linear  dimensions  of  the  parachute*.  Scheubel*4  idealised  the  inflation 
process  into  one  in  which  die  assumed  shape  of  the  canopy  remains  effectively  constant  as  its  size  increases.  His 
filling  distance  inflation  theory  which  resulted  was  extended  by  O’Hara*15,  who  adopted  a  rather  more  sophisticated 
shape  for  the  inflating  canopy  than  Scheubei  had  proposed.  It  was  then  further  developed  by  others  who  established 
relatively  simple  and  effective  inflation  theories.  Since  Scheubel's  hypothesis  is  generally  valid  and  since  the 
empirical  relationships  upon  which  these  inflation  theories  depend  most  fit  the  circumstances  for  which  they  have  been 
formulated,  these  empirical  inflation  theories  are  reliable  and  have  been  very  widely  adopted.  However,  Roberts  and 
Reddy13  have  commented  that  their  essential  weakness  rests  in  their  acceptance  as  a  necessary  input  of  the  shape  which 
the  inflating  canopy  adopts,  rather  than  this  time-dependent  canopy  shape  being  determined  instant-by-instant  as  a 
significant  output  from  the  inflation  calculation. 

5.3.1  Semi-empirical  Inflation  Models  Based  on  the  Filling  Distance  Concept 

In  a  wind  tunnel  test  at  a  constant  relative  wind  speed  V,  i.e.  under  infinite  mass  conditions,  consider  the 
instantaneous  peak  opening  force  F,  to  be  measured  on  the  inflating  canopy  compared  with  the  corresponding  steady- 
state  force  Fc,  measured  on  the  fully-inflated  canopy.  The  ratio  of  the  peak  to  steady  force  is  called  the  opening  force 
coefficient  C,,  thus 


F,  =  VV’fCoS), 


,/JpV,(CDS).C, 


where  (CoS),  is  the  canopy's  instantaneous  drag  area  and  (CoS),  is  the  drag  area  of  the  fully  inflated  canopy. 

Thus 

C,  =  (C„S),  /  (CBS), .  (5.4) 

Figure  52  shows  that  during  canopy  inflation  in  free  flights .e.  under  finite  mass  conditions  the  peak  opening  force 
is  considerably  reduced  from  its  infinite  mass  value.  The  ratio  of  its  magnitudes,  under  infinite  mass  conditions  to 
finite  mass  conditions,  is  called  the  opening  force  reduction  factor  and  this  is  denoted  by  the  symbol  X,. 

Thus  in  free  flight  the  peak  opening  force  is  written  as 


F. 


VjpV*  (CpSJ.CjXi  (5.5) 


44 


where  V,  is  the  snatch  velocity,  Le.  the  velocity  of  the  parachute  system  at  the  instant  when  the  snatch  force  is 
developed.  Very  approximately,  this  is  the  velocity  at  which  the  canopy  deploys. 

To  use  equation  5.5  it  is  first  necessary  to  determine  both  C,  andX,.  Knacke1* 54  argues  that  for  a  given 
canopy  shape  the  opening  force  reduction  factor  X,  is  a  function  of  the  canopy  loading  W/(CdSJ,  where  W  denotes  the 
payload  weight.  Since  the  canopy  loading  has  dimensions,  its  units  must  be  specified. 

Typical  values  quoted  by  Knacke  for  X,  are  1.0  for  an  aircraft  decelerator  parachute  with  a  canopy  loading  of  14 
kPa;  033  for  a  parachute  retarder  of  ordnance  supplies  with  a  canopy  loading  of  some  200  Pa  and  as  little  as  0.03  for  a 
personnel  parachute  with  a  canopy  loading  of  only  25  Pa. 

S.3.1.1  The  Mass  Ratio  Method  In  the  mass  ratio  method,  the  coefficients  C,  and  X,  are  combined  to  form  an 
instantaneous  shock  factor,  x,  thus: 


*,  =  C„Xt  »  {V*(CaS)>}/  {Vi  (CDS).}.  (5.6) 

The  factor  x,  is  the  ratio  of  the  peak  interactive  force  developed  during  inflation  to  the  interactive  force  when  the 
canopy  is  fully  inflated.  In  Section  5.3.3  the  load  factor  X  will  be  further  described.  Unlike  x„  the  load  factor  is 
defined  as  the  instantaneous  value  of  the  ratio  of  the  aerodynamic  force  developed  on  the  inflating  canopy  to  the  force 
in  die  steady  state.  The  load  Actor  is  thus  a  function  of  time  and  at  its  maximum  value  it  is  equal  in  magnitude  to  the 
opening  force  Ac  tor  x,. 

Following  Schilling15  the  opening  force  factor  is  considered  to  be  an  empirical  function  of  the  mass  ratio  R_.  The 
mass  ratio  is  a  measure  of  the  air  mass  included  in  the  fully-inflated  parachute  canopy  to  the  store  mass  m,  where  m  *> 
W/g. 


Hence 

R.  =  (p(CBSJ,J]/m  . 

(5.7) 

and 

Ck  s  function  (R„) . 

(5.8) 

5-3.1-2  The  C»nnnv  leading  Method  In  the  canopy  loading  method,  the  values  for  X,  ate  given  as  functions  of 
canopy  loading  in  the  manner  explained  in  Section  53.1.  Then,  following  Knacke.  the  opening  force  coefficient  C„ 
can  be  considered  to  be  a  function  of  the  canopy  shape  only,  hence  for  different  canopy  shapes  C,  can  be  tabulated.  For 
example,  far  flat  circular  canopies  Knacke  quotes  C,  as  equal  to  0.7. 

5.3.0  The  Pflanz  Method  Fallowing  Pflanz*4,  in  the  manner  outlined  above  the  opening  force  coefficient  C.  is 
considered  to  be  a  function  of  canopy  shape.  Then,  once  the  increasing  drag  area  of  the  inflating  canopy  has  been 
modelled  as  a  function  of  lime  by  one  of  a  number  of  simple,  definable  relationships,  for  a  rpecific  canopy  shape  the 
opening  force  reduction  factor  X,  is  given  as  an  empirical  function  of  a  relationship  whose  magnitude  can  be 
detennined 

In  any  of  these  semi-empirical  inflation  methods,  once  the  peak  opening  force  has  been  determined  the  filling  time 
can  be  estimated.  Using  Scbeubel's  concept  that  the  filling  distance  s,  is  expressed  in  terms  of  the  canopy's  nominal 
diameter  D.  by: 


»r  =  nD.  (5.9) 

where  th efill  constant  n  can  be  determined  experimentally  for  a  given  type  of  parachute,  the  filling  time  t,is  then 
expressed  in  terms  of  the  velocity  at  deployment  V,  as: 


t,  =  tfy,  =  (nDJ/V, .  (5.10) 

5JJ  More  Sophisticated  Filling  Time  Models 

By  using  the  two  equations  quoted  in  figure  5.1,  a  first  order  linear  differential  equation  can  be  written  for  the 
motion  along  the  flight  path  of  the  parachute  and  payload  system. 

Since  the  interaction  force  F,  is  there  expressed  through: 

F,  -  Dr  =  m,V  ♦  m,V 

and  •  F  i  •  D,  =>  mj  sin  yt  +  ne,V  , 

by  neglecting  both  the  payload  drag  compared  with  that  of  the  canopy  and  the  component  along  the  flight  path  of 
the  store  weight,  then: 


-  D,  «  (■,+*>  ♦  m,V. 


(5.11) 


l 
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During  the  inflation  process,  the  drag  of  the  canopy  D,  could  be  estimated  as  a  function  of  time  by  adopting  steady 
values  of  drag  coefficient  corresponding  to  simulated  canopy  shapes.  In  this  equation  the  mass  of  the  canopy  is 
inclusive  of  its  time-varying  added  mass.  In  order  to  obtain  its  time  rate  of  change,  some  estimation  must  be  made  for 
the  rate  of  change  of  canopy  volume  during  the  inflation  process.  Heinrich17  and  others  have  suggested  that  this  be 
done  through  the  application  of  the  continuity  principle: 

m  vtanRi  •  v#.,'/2JtD\  (5.12) 

where  V  represents  the  rate  of  change  of  fluid  within  the  canopy  control  volume; 

VfattRi  represents  the  inflow  to  the  canopy,  with  velocity  v„  through  the  canopy  mouth,  radius 

and  v^'/jnDj  represents  the  outflow  from  the  canopy,  with  velocity  vM  through  a  porous  hemispherical  canopy 
of  projected  diameter  Dr 

There  remains  the  problem  of  estimating  the  inlet  flow  velocity  vh.  Heinrich,  utilising  the  filling  time  concept, 
assumed  this  to  be  an  empirical  function  of  the  Ailing  time  t,,  determining  this  function  from  the  results  of 
appropriate  wind  tunnel  tests. 

Wolf44  has  argued  that  by  using  such  a  continuity  expression,  requiring  empirical  inputs,  the  prediction  cf 
canopy  inflation  has  become  unnecessarily  restricted.  What  such  filling  time  models  have  neglected  is  the  prospect  of 
a  dynamic  relationship  in  which  the  shape  of  the  inflating  canopy  would  be  determined  instant-by-instant  through 
knowledge  of  the  radial  component  of  fluid  momentum  driving  force,  which  depends  on  the  pressure  difference  across 
the  inflating  canopy,  In  1951  Wemig31  had  introduced  such  ideas  and  these  were  later  developed  both  by  Toni110 
and  by  Roberts111. 

5.3.3  Ludtke's  Parachute  Opening  Force  Analysis 

In  reference  5.12,  by  introducing  the  dimensionless  ballistic  mass  ratio  M,  Ludtke  throws  some  fresh  light  on  the 
interrelationship  between  velocities  and  aerodynamic  faces  during  the  inflation  process.  Neglecting  the  mass  of  the 
parachute  compared  with  that  of  the  payload,  in  equation  5.1: 

F,  =  D,  (5.13) 

Then,  from  equations  5.2  and  5.13,  together  with  the  equation  of  motion  fa  the  payload  in  Fig.  5.1.  assuming  that 
both  the  payload  drag  and  the  angle  of  descent  y,  are  negligible: 

-  D,  *  -  VjpV1  (SC„),  *  m.V  .  (5.14) 

Integrating  equation  5.14  with  respect  to  time  from  t  =  0,  at  the  peak  value  of  the  lines-taut  snatch  face  and  the 
velocity  is  the  snatch  velocity  V„  to  t  =  t,  when  the  canopy  is  fully-inflated  and  its  velocity  is  V: 

/J  (SC*),  dt  =  (2m,  Ip)  li  (-V  /  V1)  tit  (5.15) 

and  since  V  dt  =«  dV: 

=  (2m, /p)  Jv,  (-  dV  /  V1)  .  (5.16) 

Dividing  equation  5.16  through  by  (CoS).  teV, ,  where  to  is  the  inflation  time,  and  by  writing: 

M  «  (2m,/p  (C,S).t,VJ  (5.17) 

where  M  is  called  the  dimensionless  ballistic  mass  ratio.  Equation  5.16  can  then  be  solved  in  terms  of  M  to  give: 

V/Vs  a  1  (5.18) 

(1  ♦  <l/Mt,)Ji  l(SC,V(SC.)J  dt) 

During  inflation  the  expression  &  [(SCe)/(SCD)J  dt}  is  a  known  constant  fa  a  given  type  of  canopy  so,  once 
a  value  has  been  ascribed  to  the  ballistic  mass  ratio,  the  ratio  V/V,  can  be  determined.  Then,  having  calculated  V/V,, 
other  characteristics  of  the  inflation  phase,  such  as  the  instantaneous  shock  facta  x,,  can  be  evaluated.  Ludtke 
argued  that  the  dimensionless  ballistic  mass  ratio  is  the  roost  appropriate  scaling  parameter  with  which  10  consider  the 
aerodynamic  characteristics  of  parachutes  during  the  inflation  process. 
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5.3.4  A  Dynamic  Similitude  Model  for  Parachute  Canopy  Inflation 

Using  dimensional  analysis,  for  a  given  parachute  system  Lingard*-^13  developed  a  semi-empirica'  method 
of  predicting  over  the  entire  operational  envelope,  from  a  limited  number  of  field  trials,  the  relationship  between  the 
total  interactive  force  F„  or  the  inflation  load,  and  the  time  which  elapses. 

Assuming  that  the  drag  D,  of  the  parachute  canopy  at  an  instant  in  time  is  a  function  of: 
canopy  shape  and  size; 

its  instantaneous  velocity  V  and  acceleration  V; 
the  density  p  and  viscosity  4  of  the  fluid  in  which  it  is  inflating; 
then,  by  further  assuming  that  the  acceleration  V  is  a  function  of: 
the  snatch  velocity  V,; 

the  time  t,  measured  from  the  instant  when  the  snatch  force  is  developed; 

the  canopy  drag  D,; 

the  gravitational  acceleration  g; 

the  masses  of  the  payload  m,  and  the  canopy  ny 

the  canopy  rate  of  change  of  mass  riy,  and 

the  instc  uaneous  angle  of  the  parachute  trajectory  relative  to  the  vertical,  or  the 
deployment  angle  0 ; 


then: 


D,  =  function  {shape;  D.  ;  V,  ;  t  ;  p  j  4  ;  g  ;  m,  ;  m,  ;  0  }  .  (5.19) 

For  a  canopy  of  a  given  shape,  by  neglecting  the  effects  of  4,  m,  and  m,  Lingard  showed  that  for  geometrically 
similar  parachute  systems  with  similar  porosity  constants: 

c;  =  D/pDjV1  =  functional,;  Fj  t ;  0  }  (5.20) 

and  X  =  D,/m,g  =  function (M, ;  F ;  x  ;  0  }  (521) 

whe.eCp  is  the  dimensionless  aerodynamic  force  developed  on  the  canopy; 

X  is  the  load  factor,  defined  as  the  ratio  of  the  instantaneous  to  the  steady  aerodynamic  farce  developed  on  the 
parachute  canopy; 

M,  is  the  mass  ratio,  here  defined  as  the  ratio  of  payload  mass  m,  to  a  mass  representative  of  that  included 

within  the  canopy,  pD3 ; 

F  is  the  Froude  number  V^gD.  defined  in  Section  23.4; 

1  is  the  dimensionless  time,  V,t/D. ; 

0  is  the  deployment  angle,  i.e.  the  instantaneous  angle  of  the  parachute  trajectory,  relative  to  the  vertical. 

Smce  he  was  primarily  concerned  with  personnel  parachutes,  in  his  analysis  Lingard  did  not  consider  the  Mach 
number  to  be  a  significant  independent  parameter. 

Unsteady  inflation  force  data  obtained  from  experiments  conducted  on  a  variety  of  canopies  tested  over  a  range  of 
mass  ratios  and  descent  parameters  correlated  well  when  {dotted  in  the  form: 

Cj  =  function  (x)  (522) 

Lingard513  therefore  concluded  that  each  canopy  shape  has  a  unique  dimensionless  inflation  force/time  signature, 
which  can  be  extracted  from  a  limited  number  of  trials  of  a  given  system.  By  employing  this  .nform?tion,  together 
with  the  application  of  Newton's  laws  of  motion,  as  shown  in  Chapter  3,  to  the  system  it  is  possible  to  predict  the 
performance  of  the  parachute  and  the  payload  system  over  its  entire  operational  envelope. 

5.3.5  Xinctic  Models  for  Parachute  Canopy  Inflation 

The  inflation  method,  introduced  in  1971  by  Robots,  represented  the  parachute  canopy  by  a  continuous  elastic 
system.  The  pressure  distribution  required  to  determine  the  force  which  causes  the  necessary  rate  of  change  of  fluid 
momentum  was  calculated  by  assuming  potential  flow  about  an  expanding  and  decelerating  parabolic  shell, 
representative  of  the  canopy.  B 'cause  of  its  geometrical  and  mathematical  sophistication  3uch  an  advanced  model  is 
difficult  to  apply  and  thus  receives  only  partial  acceptance  within  the  parachute  indusuy. 

Provided  that  some  overall  data  for  canopies  during  their  inflation  phase  have  been  obtained  by  conducting 
appropriate  experiments.  e.g.  the  peak  opening  force  coefficient  and  the  dimensionless  filling  time,  together  with  their 
payload  and  canopy  mass  ratios,  Wolfs”  1973  single  degree  of  freedom  canopy  model  solves  the  necessary 
momentum  equations  and  satisfactorily  predicts  some  observed  phenomena,  such  ar,  the  effects  of  altitude  on  parachute 
inflation  time.  The  mass  ratios  are  defined  here  as  the  respective  payload  or  canopy  mass,  divided  by  the  mass  of  air 
displaced  by  a  fully-inflated  representative  spherical  canopy. 
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Different  methods  were  adopted  by  Roberts5-11,  by  Klimas?'*1115  ard  by  others  in  order  to  identify  the 
pressure  distribution  in  the  flow  field  associated  with  the  inflating  canopy.  These  arc  discussed  in  greater  detail  in 
Section  7.1.  In  1981  Purvis5-15  developed  an  analytical  model  based  on  a  simplifeation  of  the  changing  canopy 
shape  during  the  inflation  process.  This  model  needs  no  experimental  inputs  at  all.  It  could  be  used  to  predict  the 
trajectory  of  an  inflating  canopy. 

In  his  model  the  canopy  is  modelled  as  a  right  circular  cylinder  whose  radius  is  free  to  increase  with  time.  For  the 
inflation  of  this  cylinder  immersed  in  an  inviscid  and  incompressible  fluid  an  expression  for  the  time  rate  of  change  of 
axial  momentum  over  the  surface  of  the  expanding  cylinder  is  first  established  and  then  solved.  First-order  effects  only 
arc  considered.  In  reference  5.16  Purvis  made  comparisons  for  both  im porous  and  porous  canopies  between  the  results 
obtained  with  this  model  and  experimentally-obtained  data,  inflating  under  both  infinite  mass  and  finite  mass 
conditions. 

The  real  significance  of  analytical  inflation  models,  such  as  th;s  one  developed  by  Purvis,  is  that  they  reveal  which 
are  the  gross  parameters  governing  the  inflation  process,  indicating  what  may  be  the  consequences  of  their  independent 
variation. 
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6.  EXPERIMENTS  TO  DETERMINE  PARACHUTE  AERODYNAMIC  CHARACTERISTICS 

At  first  sight,  the  proper  manner  in  which  to  conduct  experiments  in  order  to  establish  the  aerodynamic  characteristics 
of  parachutes  would  appear  to  be  through  the  use  of  full-scale  proto  ypes  deploying,  inflating  and  descending  through 
their  natural  environment  Indeed,  some  experiments  are  conducted  in  this  way.  Many  others  are  performed  in  wind 
tunnels  or  similar  environments,  both  because  these  can  readily  be  controlled  and  because  under  these  circumstances  it 
is  often  much  easier  to  obtain  the  required  data  by  instrumenting  the  support  required  for  the  model  parachute  than  it  is 
to  determine  them  from  Eight  tests.  Wind  tunnels  are  most  often  used  for  both  static  and  dynamic  tests  on  model 
parachutes,  as  a  means  of  enabling  the  flow  around  the  parachutes  to  be  visualised  and  for  measuring  as  a  fuction  of 
canopy  attitude  the  aerodynamic  forces  which  are  developed  Other  facilities  which  have  been  used  for  special-purpose 
tests  include  water  tunnels,  in  which  air  is  replaced  with  water  as  the  test  medium  in  which  the  parachute  is  immersed. 

6.1  EXPERIMENTS  CONDUCTED  ON  FULL-SCALE  PARACHUTES  IN  AIR 

As  described  in  Section  23.3,  the  drag  area  QjS.  of  a  descending  parachute  can  be  determined  from  a  knowledge  of 
the  weight  of  the  parachute  system,  the  density  of  the  air  through  which  it  descends  and  the  descent  velocity  V„.  The 
latter  can  be  crudely  estimated  by  observation.  James*1,  for  example,  hung  a  61-m  long  axial  cord  below  the 
payload  and  observed  the  time  which  elapsed  between  the  two  ends  of  this  cord  striking  the  ground.  However, 
Drake*1,  by  analysing  kindtbeodolite  data  over  a  500ft  descent,  used  a  more  sophisticated  technique  with  which  to 
estimate  the  descent  velocity. 

The  drag  area  so  obtained  over  the  time  of  descent  would  be  an  average  value.  Since  the  canopy  nominal  surface 
area  S.  is  known,  the  average  drag  coefficient  for  the  system  is  determined.  For  stable  parachutes  this  measurement 
accords  well  with  wind  tunnel  evaluations,  but  as  the  angle  of  attack  of  an  unstable  parachute  varies  continuously 
during  its  descent,  for  such  a  parachute  there  may  well  be  a  substantial  discrepancy  between  the  average  value  of  C0 
obtained  during  its  free  descent  through  the  air  and  the  corresponding  wind  tunnel  evaluation  of  Q,  at  zero  angle  of 
attack. 

Depending  on  their  application,  parachutes  being  tested  may  be  allowed  to  inflate  freely,  having  been  dropped  from 
tethered  balloons,  they  may  be  ejected  from  aircraft,  using  specially-designed  test  vehicles  such  as  those  which  have 
been  described  by  Key  and  Barker*1  or  they  could  be  launched  from  ground-based  test  vehicles  such  as  the  British 
compressed-air  launcher,  recently  being  used  at  the  Royal  Aircraft  Establishment.  Such  facilities  can  be  instrumented 
to  telemeter  appropriate  data  from  the  parachute  system  to  the  ground.  A  test  vehicle  described  by  Barker  and 
Nosworthy**  was  designed  to  be  dropped  from  tethered  balloons  which  fly  at  altitudes  up  to  1  000  m.  This  test 
vehicle  fell  freely  and  as  it  did  so,  after  pre-set  time  intervals  the  lest  parachutes  were  deployed  from  it  Lingard  used 
this  test  vehicle  to  obtain  required  the  canopy  inflation  dimensionless  time  signature  described  in  Section  5.3.4. 

6.2  EXPERIMENTS  IN  A  CONTROLLED  ENVIRONMENT 

Although  when  the  weights  of  both  canopy  and  payload  are  known  the  drag  coefficients  for  stable  parachute  canopies 
can  be  determined  by  measuring  their  rate  of  descent  through  the  atmosphere,  to  acquire  more  sophisticated 
aerodynamic  data  the  instantaneous  angle  of  attack  of  the  canopy  will  also  be  required.  To  deduce  this  angle  both  the 
magnitude  and  the  direction  of  the  relative  airflow  must  be  obtained.  The  most  satisfactory  way  of  determining  these 
data  is  to  fly  the  canopy  in  a  controlled  environment  in  which  the  relative  velocity  of  the  fluid  win  be  known.  Thus 
the  main  reason  for  conducting  measurements  on  scale  models  of  parachutes  in  wind  tunnels  is  in  order  to  provide 
such  a  controlled  environment.  It  is  highly  desirable  to  use  wind  tunnels  rather  than  flight  tests  when  obtaining  drag 
data  for  unstable  parachute  canopies  end  when  determining  canopy  static  stability  characteristics.  For  such 
measurements  static  tests  can  be  performed  on  rigidly-mounted  canopies.  Similarly,  some  form  of  controlled 
environment  is  preferable  when  canopy  inflation  behaviour  or  dynamic  stability  characteristics  are  sought,  but  for 
these  purposes  dynamic  testing  is  necessary,  on  models  which  can  move  through  the  controlled  environment  with  a 
limited  amount  of  freedom. 

6.3  REQUIREMENTS  FOR  MODEL  TESTS  CONDUCTED  IN  A  CONTROLLED 
ENVIRONMENT 

To  model  the  airflow  round  a  descending  parachute  faithfully  when  it  is  flying  in  a  controlled  environment  such  as 
that  in  a  wind  tunnel  two  important  requirements  must  be  satisfied: 

(i).  the  shape  of  the  model,  including  the  means  of  fixing  it  to  any  force-measuring  apparatus  and  to  the  wails  of 
the  controlled  environment,  together  with  the  canopy  porosity  and  its  flexibility,  must  be  truly 
representative  of  the  prototype  full-scale  parachute.  The  scaling  of  canopy  flexibility,  as  Lceuhas 
indicated,  is  particularly  difficult  to  achieve  and  this  can  lead  to  problems  in  data  interpretation, 
particularly  for  inflation  loads.  This  shape  requirement  also  includes  ensuring  that  any  blockage 
constraint  which  is  caused  by  the  presence  of  the  controlled  environment  walls  is  minimal; 
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(it),  the  Reynolds  number  of  the  test  programme  and  also  its  Mach  number,  where  it  is  applicable,  must  also  be 
representative  of  the  full-scale  parachute.  As  discussed  in  Sections  222  and  2.2.3  these  parameters 
need  not  necessarily  be  equal  to  those  in  full-scale  flight  but  in  the  model  tests  any  differences  in  their 
magnitudes  must  be  considered  carefully. 

Inevitably,  when  conducting  model  tests  the  shape  requirement  will,  to  some  extent,  be  compromised.  Early 
experiments  were  often  performed  on  rigid  rather  than  flexible  models  of  canopies  and  in  some  cases  the  tests 
conducted  on  these  canopies  were  not  performed  at  sufficiently  high  Reynolds  numbers  to  avoid  laminar  boundary 
layer  separation.  Often  in  static  tests  models  are  mounted  on  axial  stings  which  can  both  limit  the  movement  of  the 
model  canopies  and  may  develop  drag  forces  on  their  own  account  Because  parachute  canopies  are  bluff  rather  than 
streamlined  bodies  any  effects  on  aerodynamic  characteristics  caused  by  blockage  constraint  can  be  of  considerable 
significance. 

6.4  WIND  TUNNEL  TESTS  ON  MODEL  PARACHUTES 

The  scale  effect  and  blockage  problems  referred  to  above  were  mentioned  as  long  ago  as  1946,  in  Block's**  brief 
report  Early  German  wind  tunnel  tests  which  Munson*7  described  are  principally  concerned  with  the  establishment 
of  the  proper  dimensionless  parameters  influencing  the  aerodynamic  forces  developed  on  parachute  canopies,  also  with 
determining  both  static  stability  requirements  and  the  opening-shock  forces.  Heinrich's  wind  tunnel  tests,  originally 
conducted  in  Germany  and  later  in  the  United  States,  have  been  described  in  a  variety  of  reports,  such  as  references  1.8, 
2.2, 6.8  and  6.9. 

Later  German  research  has  been  considered  by  both  Doherr,  in  references  6.10  and  2.10.  and  by  Saliaris*11.  In 
this  experimental  work  the  techniques  which  were  adopted  for  static  tests  can  be  considered  as  a  development  of  those 
which  Heinrich  had  earlier  implemented  in  the  United  States.  Experimental  methods  to  determine  parachute  dynamic 
stability  characteristics  were  also  developed  in  Germany. 

Although  much  of  the  more  recent  British  parachute  testing  has  been  performed  in  the  free  air,  in  reference  l.S 
Dennis  refers  to  some  wind  tunnel  testing  in  the  United  Kingdom.  Other  recent  British  experimental  work  has  been 
described  by  Shen  and  Cockrell  in  reference  2.11. 

6.4.1  Flow  Visualisation  Around  Model  Parachutes 

Wool  tufts  fixed  to  detect  the  onset  of  flow  separation  from  model  parachute  canopies  and  smoxe  employed  as  a  flow 
tracer  are  the  most  commonly-used  wind  tunnel  techniques  for  flow  visualisation  around  parachute  canopies,  though 
because  of  the  rapid  dissipation  of  the  smoke  the  latter  is  not  a  very  appropriate  technique  at  Reynolds  numbers  which 
approach  full  scale. 

Techniques  for  using  neutrally-buoyant  helium-filled  soap  bubbles  as  flow  tracers  around  model  parachute  canopies 
have  been  described  by  Pounder*11.  Klimas  A  Rogers^1,  Lingard*'  and  by  Shen  &  Cockrell111. 

One  of  the  major  advantages  of  testing  parachutes  under  water  rather  than  in  air  is  that  this  makes  possible  the  use 
as  flow  tracers  of  either  small,  near  neutrally-buoyant,  polystyrene  beads  and  this  technique  has  been  described  by 
Lingard51, or  of  hydrogen  bubbles  generated  by  local  electrolysis  at  fine  wire  cathodes  immersed  in  the  water.  This 
latter  technique  has  been  outlined  by  Cockrell,  Huntley  and  Ayres*1*. 

When  testing  model  parachute  canopies  under  water  scaling  problems  can  arise.  In  particular,  as  Cockrell.  Harwood 
and  Shen*11  have  discussed,  the  nominal  porosity  Xof  a  parachute  tested  under  water  can  differ  appreciably  from  its 
value  when  determined  in  air. 

6.4.2  Measurement  of  Steady  Aerodynamic  Forces  and  Moments 

In  order  to  measure  steady  aerodynamic  forces  and  moments  developed  on  model  parachutes,  generally  the  models  are 
rigidly  fixed  in  their  test  media.  By  using  strain  gauges  or  other  appropriate  transducers  the  required  aerodynamic 
reactions  on  the  supporting  structure  can  be  determined.  Such  a  measuring  technique  has  been  well  described  both  by 
Heinrich  A  Haak11  and  by  Doherr140. 

6.4.3  Unsteady  Aerodynamic  Measurements 

Unsteady  aerodynamic  measurements  have  been  made  on  model  parachute  canopies  in  order  to  determine  both  the 
opening  loads  during  inflation  and  the  added  masses  which  are  developed  on  unsteadily-moving  fully-inflated  canopies. 

Measurements  of  IntcratttvcBaraa  Owing  Inflation  ante  Finite  Mm  Condi  liens  ifmodei 
parachute  canopies  which  are  rigidly  fixed  in  wind  tunnels  are  inflated,  this  process  will  occur  at  a  constant  wind 
velocity  relative  to  the  model.  In  Section  S.1.I  it  is  explained  that  infmite-mass  measurements  made  under  these 
conditions  of  the  interactive  force  F,  developed  between  the  parachute  canopy  and  the  payload  are  unrepresentative  of 
full-scale  canopy  inflation.  However,  both  Heinrich  A  Noreen*4*  and  Lingard14  have  shown  that  by  mounting 
the  canopy  in  the  wind  tunnel  so  that  as  it  inflates  it  is  five  to  move,  then  during  its  inflation  the  surrounding  air  flow 
will  decelerate  relative  to  the  canopy  and  the  finite -mass  iitflatkxi  process  cut  be  modelled. 
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Fig.  6.1a  FINITE  MASS  WIND  TUNNEL  MEASUREMENT  OF  THE 
INTERACTIVE  FORCE  F  FOR  AN  INFLATING  CANOPY 


By  appropriately  modelling  the  system  shown  in  Fig.  6.1a,  its  ensuing  behaviour  can  be  designed  to  be 
representative  of  that  around  a  vertically-descending  payload  in  the  free  air,  decelerating  under  the  influence  of  a 
parachute  undergoing  a  finite  mass  inflation. 

For  the  purposes  of  analysis  Fig.  6.1b  shows  the  same  system  broken  into  a  number  of  separate  free-bedy 
diagrams. 
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The  equations  of  motion  for  each  of  these  sub-systems  arc  then  as  follows: 

F,  -  T,  =  l,wV 

(T,  -  TJr.  =  laV/r 

T,  +  l,W8  •  T,  a  l,wV 

(T«  -  T,)rk  *  I.V/rk 

T,  -  T,  =  IjwV 

(T<  -  T,)rt  =  IfV/r, 

T,  -  l«wg  -T,  s  I4WV 

T,  -  m^j  s  m,V 


equation  0): 
equation  (ii); 
equation  (iii); 
equation  (iv); 
equation  (v); 
equation  (vi); 
equation  (vii); 
equation  (viii). 


The  solution  to  this  set  of  equations  is: 

V(w(l,+1,+1,+1,)  +  ljrl+ljr}+ljr*  +  m,} 


Fi  *  +  m. 


(6.1) 


This  equation  is  of  the  same  form  as  that  for  the  vertically-descending  payload  shown  in  Fig.  6.1b  where  in  both 
equations  F,  represents  the  required  interactive  force  between  the  payload  and  the  parachute. 

Thus,  with  an  appropriate  choice  of  the  constants  for  the  wind  tunnel  system  equation  6.1  can  be  made  to  represent: 

Vm,  =  F,  +  D,  -  m,g  (6.2) 

and  so  the  interactive  force  can  be  determined. 

6A3.2  Experiment  to  Determine  the  Added  Maanfor  Unfleadily-MoyingJullyrlnflalctl  Parachute  Canopies  in 
order  to  determine  these  added  mass  coefficients  experimentally,  Cockrell,  Shen,  Harwood  and  Baxter*  3  measured 
the  aerodynamic  farces  which  were  developed  on  fully-inflated  model  parachute  canopies  when  the  latter  were  subjected 
to  a  known  periodic  motion.  The  aerodynamic  forces  were  calculated  from  output  signals  transmitted  from  strain 
gauges  which  were  attached  to  the  models'  support  sling. 

6.5  AERODYNAMIC  MEASUREMENTS  MADE  ON  PARACHUTES  IN  OTHER  FACILITIES 
THAN  IN  WIND  TUNNELS 

From  time  to  time  aerodynamic  measurements  on  parachute  canopies  have  been  made  in  other  facilities  than  wind 
tunnels.  For  example,  at  some  time  prior  to  1967  measurements  which  Colboume*17  has  described  were  made  of 
the  drag  developed  by  a  4.6  m.  (IS  ft)  flat  circular  parachute  when  it  was  caused  to  descend  freely  inside  the  107  m. 
(350  ft.)  high  cooling  tower  of  an  electricity  power  station.  This  tower  had  a  base  diameter  of  100  m.  (325  ft),  a 
diameter  at  the  apex  of  66  m.  (218  ft.)  and  a  throat  diameter  of  62  m.  (205  ft.).  Its  varying  diameter  was  a  source  of 
some  difficulty  in  that  it  caused  a  corresponding  variation  in  descent  velocity.  This  facility  was  only  available  over  a 
very  limited  period  for  experimental  purposes  and,  probably  in  consequence.  Col  bourne's  report  on  the  measurements 
which  he  made  is  somewhat  inconclusive. 

As  described  in  reference  6.15,  in  order  to  minimise  the  ratio  of  the  inertia  forces  developed  on  the  canopy  supports 
to  the  added  masses  developed  on  parachute  canopies  when  they  move  unsteadily,  measurements  of  the  aerodynamic 
characteristics  for  fully-inflated  parachute  canopies  have  been  made  under  water  rather  than  in  the  air.  Water  is  a 
suitable  medium  because  the  aerodynamic  forces  are  proportional  to  the  fluid  density  and  that  of  waiter  is  some  800 
times  the  density  of  air,  whereas  the  inertia  forces,  being  proportional  to  the  density  of  the  canopy  supports,  are  of  the 
same  order  of  magnitude  whichever  medium  is  adopted.  In  the  experiments  described  the  canopy  models  were  towed 
through  a  61.0  m.  long  ship  tank,  having  been  suspended  from  the  ship-towing  carriage. 

6.6  BLOCKAGE  CAUSED  BY  MODEL  PARACHUTE  CANOPIES 

The  flow  past  any  body  immersed  in  a  stream  of  fluid  is  subjected  to  blockage  constraint,  caused  by  restraint  of  the 
fluid's  free  lateral  displacement  This  stream  of  fluid  might  be  constrained  either  by  the  test  environment's  solid 
walls,  as  it  is  when  the  body  is  held  in  the  closed  working  section  of  a  wind  tunnel  or  in  a  less  confined  situation 
such  as  would  occur  if  the  model  were  immersed  in  an  open  jet  of  fluid.  In  the  former  state  the  boundary  condition 
imposed  on  the  flow  is  that  at  the  solid  boundary  there  can  be  no  transverse  velocity  component:  in  the  latter  state  the 
corresponding  boundary  condition  is  that  the  pressure  along  the  boundary  must  be  approximately  constant  and  equal  to 
the  ambient  pressure.  In  these  two  states,  the  corrections  which  are  to  be  applied  to  measured  pressures,  aerodynamic 
forces  and  moments  are  of  opposite  sign.  Maskcll441  shows  that  the  dominant  effect  of  blockage  constraint  is  a 
simple  increase  in  the  fluid's  free-strearo  velocity,  in  part  related  to  the  volume  distribution  of  the  body  itself,  termed 
solid  blockage,  and  in  part  related  to  the  dirpiacement  effect  of  the  wake,  termed  wake  blockage. 

Conventional  parachute  canopies  are  bluff  bodies,  that  is  bodies  for  which  the  surrounding  flow  is  dominated  by 
large  regions  of  flow  separation.  Unless  the  blockage  area  ratio,  defined  as  the  ratio  of  the  cross-sectional  area  of  the 


body  S  to  toe  cross-sectional  area  of  the  fluid  stream  at  the  body  location  A,  is  small  (say  under  15%),  the  variation  in 
the  flow  pattern  around  the  body  and  its  downstream  wake  which  is  caused  by  blockage,  is  extreme.  In  the  published 
literature,  e.g.  Blockage  Corrections  for  Bhtf  Bodies  in  Confined  Flows,  reference  6.19,  the  emphasis  is  on  correction 
of  drag  measurements  for  blockage  constraint  and  because  Muff  bodies  are  inefficient  generators  of  lift,  little  has  been 
written  about  the  correction  of  aerodynamic  force  components  developed  on  them  which  are  normal  to  the  direction  of 
the  relative  flow.  An  estimate  from  reference  6.19  for  this  bluff  body  across-flow  force  coefficient  blockage  correction 
AC,  is; 

(ACc )/Cc  -  -  0.2*/A  (6.3) 

and  for  a  parachute  canopy  presenting  as  large  a  blockage  ratio  in  a  wind  tunnel  (or  similar  facility)  as  15%,  this 
gives  a  negligible  •  3%  across-flow  force  coefficient  blockage  correction . 

Reference  6.19  states  that  for  bluff  bodies  such  as  parachute  canopies  blockage  areas  of  5%  and  under  are  to  be 
prefened;  blockage  areas  in  excess  of  10%  are  not  to  be  recommended.  For  blockage  ratios  under  10%  corrections  to 
the  drag  of  conventional  parachutes  can  readily  be  made  by  one  of  four  methods,  the  choice  of  which  depends  on  the 
location  of  flow  separation  from  the  canopy  and  the  amount  of  experimental  data  which  are  available.  This  reference 
also  contains  data  on  blockage  corrections  for  both  mean  pressure  measurements  and  for  fluctuating  quantities,  such  as 
for  Strouhal  numbers. 

Estimates  of  blockage  constraint  for  gliding  parachutes  should  be  based  on  the  more  extensive  literature  available 
concerning  blockage  corrections  for  model  wings  in  wind  tunnel  tests,  e.g.  Lift-interference  and  Blockage  Corrections 
for  Two-dimensional  Subsonic  Flow  in  Ventilated  and  Closed  Wind  Tunnels,  reference  6.20. 

6.6.1  Maskell's  Bluff-body  Blockage  Constraint  Method 

For  a  conventional  parachute  canopy,  from  which  the  flow  separates  at  or  ahead  of  its  maximum  cross-sectional  area, 
the  estimated  conection  to  drag,  AQ>,  consequentional  on  wind  tunnel  blockage  which  was  obtained  by  Maskellu* 
can  be  expressed  as; 

bCJC0  =  -  2.77  CD*/A  .  (6.4) 

Thus,  for  a  model  parachute  canopy  subjected  to  a  5%  blockage  area  ratio  in  a  wind  tunnel,  if  the  measured  drag 
coefficient  were  0.80  the  corrected  drag  coefficient  would  be  0.71. 

6.6.2  Cowdrey's  Alternative  Expression  for  Bluff-body  Blockage  Constraint 

When  values  of  drag  coefficient  are  not  readily  available  Cowdrey's121  alternative  expression  for  blockage  constraint 
can  be  used: 

ACb/Cb  =  •  3.15  */A  .  (63) 

Since  in  this  relationship  Cowdrey  was  primarily  concerned  with  very  large  bluff  bodies  such  as  buildings,  equation 
63  tends  to  overestimate  the  blockage  conection  necessary  for  parachute  canopies.  When  the  blockage  area  ratio  is 
only  5%  the  drag  coefficient  for  the  wind  tunnel  model  must  be  as  high  as  13  in  order  to  acheive  the  same  drag 
correction  from  equation  63  as  is  given  by  equation  6.4.  If  the  measured  drag  coefficient  were  0.80.  then  by  using 
equation  63  the  corrected  drag  coefficient  would  be  0.67,  instead  of  0.71  which  is  given  in  the  paragraph  above. 

6.63  The  Quasi-Streamlined  Flow  Method  for  Bluff-body  Blockage  Constraint 
If  separation  were  to  occur  downstream  of  the  canopy's  maximum  cross-sectional  area,  a  relationship  originally 
developed  by  Gainer,  Rogers,  Acum  and  MaskeU“*  gives  for  the  blockage  constraint  of  a  conventional  parachute 
canopy  mounted  centrally  in  a  wind  tunnel  which  has  a  rectangular  working  section  of  aspect  ratio  A: 


ACn/Co  =  -  034  X  (*/A)**  -  030  C„*/A. 

(6.6) 

where 

X  *  0.72(1+A)/A  . 

(6.7) 

For  a  measured  drag  coefficient  of  0.80  on  a  conventional  parachute  canopy  with  a  5%  blockage  area  ratio  in  a 
1.4:1  aspect  ratio  wind  tunnel  working  section,  the  corrected  drag  coefficient  determined  by  using  equation'  6.6  and 
6.7  is  0,77. 

6.6,4  Bluff-body  Blockage  Constraint  Determined  from  the  Working  Section  Wall  Pressure 
.  Distribution 

If  the  static  pressure  can  be  measured  on  the  wind  tunnel  wall  at  the  sections  at  which  the  blockage  is  a  maximum  and 
at  a  second  section  downstream  of  the  first  a  further  method,  described  in  reference  6. 19  and  originally  developed  by 
Garner  et  al**  can  be  used  to  estimate  the  blockage  correction. 


6.6.5  Estimation  of  Blockage  Constraint  in  Wind  Tunnel  Tests  on  Parachute  Canopies 
In  order  to  minimise  blockage  effects,  aerodynamic  tests  on  model  parachutes  must  be  performed  in  large  wind 
tunnels.  For  much  of  the  published  experimental  dam,  not  only  is  it  not  known  if  any  blockage  corrections  were 
made  by  the  originators,  but  as  wind  tunnel  blockage  area  ratios  may  not  now  be  readily  available,  making 
contemporary  corrections  is  often  not  possible.  Values  of  drag  coefficient  obtained  from  different  wind  tunnel  tests 
can  vary  substantially  from  one  another  and  although  in  published  data  the  sign  of  d(yda  may  be  reliable,  its 
numerical  value  is  much  less  certain. 

In  determining  the  aerodynamic  characteristics  of  parachute  clusters  Braun  and  Walcott*”,  using  the  3.7  m.  (12 
ft.)  diameter  vertical  wind  tunnel  at  the  Wright-Patterson  Air  Force  Base,  Ohio  found  that  their  blockage  factors 
approached  19%.  Earlier  Amerson*1*,  using  the  7.3  m.  (24  ft)  diameter  low-speed  wind  tunnel  at  the  Royal 
Aircraft  Establishment  Famborough,  UJC,  performed  tests  on  clusters  of  from  one  to  five  parachute  canopies  with 
only  a  4%  maximum  blockage  area  ratio.  Heinrich  and  Noreen*1'  tested  clusters  of  from  one  to  four  solid  flat 
circular  and  ring  slot  parachute  canopies  in  the  1.5  m.  x  1.5  m.  (5  ft  x  5  ft)  open  working  section  wind  tunnel  at  the 
University  of  Minnesota  at  Minneapolis,  experiencing  blockage  factors  of  up  to  22%.  Using  the  2.1  m.  x  3.0  m.  (7 
ft  x  10  ft)  working  section  of  the  Vought  Corporation  wind  tunnel,  clusters  of  between  one  and  eight  0.4  m.  (16  in.) 
nominal  diameter  flat  circular  canopies  were  tested  by  Baca*29.  In  these  latter  tests  the  blockage  factor  varied  from 
1.8%  to  14%.  In  his  determination  of  the  canopy  drag  coefficient  for  the  fixe  stream  dynamic  pressure  Baca  estimated 
the  effect  of  ’his  blockage  factor  by  determining  the  dynamic  pressure  variation  along  the  wind  tunnel  ceiling  outside 
of  the  parachutes'  boundary  layer  and  using  the  •.'slue  which  was  obtained  in  the  plane  closest  to  that  containg  the 
canopy  skin  hem.  His  method  is  an  approximation  io  that  outlined  in  Section  6.6.4  above. 

Shen  and  Cockrell111  measured  the  (frag  of  cross-shaped  parachute  canopies  both  in  a  wind  tunnel  with  a  working 
section  of  1.14  m  x  0.84  m  and  in  a  water  tank  which  had  a  cross-sectional  area  of  3.66  m  x  1.83  m.  Hie  blockage 
area  ratios,  based  on  the  projected  areas  of  the  canopies,  were  between  7%  and  8%  in  the  wind  tunnel  and  of  the  order 
of  a  negligible  1%  in  the  ship  tank.  Corresponding  drag  coefficient  characteristics  have  been  drawn  in  Fig.  6.2.  The 
upper  curve  shows  the  values  of  drag  coefficient  which  were  obtained  in  the  wind  tunnel.  The  corrected  wind  tunnel 
results  which  are  shown  by  the  open  symbols  in  Fig.  62  were  obtained  by  using  equation  6.4  for  the  correction  to 
drag,  ACd„  These  compare  very  well  with  the  drag  coefficient  measurements  independently  obtained  in  the  ship  tank 
with  the  much  larger  cross-sectional  area  using  the  same  model  parachute  canopies. 

- \ 

Drag  Coefficient,  Cq 


Canopy  Porosity,  in  cu.ft./sq.ft/sec. 

At  10  inches  of  Water  Pressure 

Fig.  6.2  BLOCKAGE  CORRECTIONS  APPLIED  TO  CROSS-SHAPED 
PARACHUTE  CANOPIES  -  Shea  sad  Cockrell,  ref.  2.11 
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7.  METHODS  OF  ANALYSIS  FOR  FLOW  AROUND  PARACHUTE  CANOPIES 

For  a  variety  of  reasons  analytic  solutions  are  required  to  a  number  of  parachute  aerodynamic  problems.  Because 
of  the  interactive  relationship  between  the  canopy  shape  and  the  aerodynamic  forces  developed  on  it,  even  when  it  is 
possible  to  make  experimental  measurements  they  may  well  be  insufficient  for  predictive  purposes.  Often  the 
unsteady  nature  of  the  flow  around  bluff  parachute  canopies  causes  difficulties  in  both  the  planning  and  the  execution 
of  experimental  programmes  and  even  where  many  aerodynamic  characteristics  have  been  described  in  Section  2  as 
steady,  they  are  in  fact  the  average  values  of  time-varying  quantities  which  fluctuate  with  the  unsteady  wake  pattern 
downstream.  For  these  and  other  reasons  reliable  experimental  solutions  for  parachute  aerodynamic  characteristics  can 
be  difficult  to  obtain  and  analytical  solutions  are  sought,  not  only  for  predictive  purposes  but  also  in  order  to  obtain  a 
better  understanding  of  the  fluid  flow  processes  involved. 

A  thorough  analysis  of  the  whole  flow  field  around  the  parachute  is  possible  only  by  solution  of  the  Navier- 
Stokes  equations.  Though  certain  low  Reynolds  number  flow  problems  do  respond  to  the  computer  solution  of  these 
equations  the  flow  round  a  parachute  canopy  and  payload  is  at  a  high  Reynolds  number,  implying  that  convection  in 
the  flow  field  is  of  more  significance  than  diffusion  and  that  viscosity  is  a  relatively  unimportant  fluid  property.  In 
the  immediate  future  parachute  aerodynamic  characteristics  are  most  unlikely  to  be  obtained  from  full  solutions  to  the 
Navier-Stokes  equations  and  other  methods  must  be  sought  instead.  However,  because  the  three-dimensional  flow 
field  around  a  bluff  body  such  as  a  parachute  canopy  contains  large-scale  structures  which  arise  from  free-shear  layers 
brought  about  by  flow  separation,  its  solution  is  far  from  staightforward.. 

Currently  there  are  two  possible  lines  of  attack,  in  both  of  which  a  largely  irrotational  flow  field  is  assumed.  In 
the  first  of  these  approaches,  in  spite  of  some  evident  disadvantages  in  this  assumption,  the  entire  flow  field  outside  of 
any  near-surface  singularities  is  considered  to  be  irrotational.  This  is  the  method  that  Ibrahim,  Klimas  and  Roberts  all 
adopted,  described  below  in  Section  72.  In  the  second,  described  in  Section  73,  an  identifiable  region  of  vorticity  is 
used  to  model  the  characteristics  of  the  wake  which  is  shed  by  the  parachute. 

Irrotational  flow  is  an  essential  ingredient  of  any  analytical  model  which  is  developed  in  order  to  determine  the 
characteristics  of  the  flow  around  an  arbitrarily-shaped  body,  for  it  is  only  if  the  flow  is  irrotational  that  it  possesses  a 
velocity  potential  and  hence  can  be  termed  potential  flow.  Potential  flows  can  be  steady  or  unsteady.  In  steady 
potential  flows  streamlines  can  be  drawn  orthogonal  to  the  equipotential  lines.  Vorticity  is  a  measure  of  fluid  element 
rotation  thus,  in  a  region  in  which  a  fluid  is  considered  to  be  irrotational,  there  can  be  no  vorticity.  Since  vorticity  is 
necessarily  present  in  fluid  regions  in  which  there  are  3hear  stresses,  within  boundary  layers  and  wakes  the  flow  is  not 
potential.  Thus,  it  is  an  idealisation  to  conceive  of  the  whole  flow  field  in  which  a  given  body  is  immersed  as  being 
potential.  To  contrast  this  idealisation  with  the  actual  flow  field,  in  which  vorticity  is  present  in  certain  high  shear 
regions,  the  potential  flow  field  is  also  referred  to  as  ideal  fluid  flow,  or  sometimes  as  perfect  fluid  flow.  It  is 
customary,  though  not  essential,  to  consider  that  a  characteristic  property  of  an  ideal  fluid  is  its  incompressibility. 

Certain  features  of  unsteady  flows  around  bluff  bodies  such  as  parachute  canopies  can  be  predicted  with  remarkable 
accuracy  by  using  steady  flow  models  but  with  such  methods  a  cautious  approach  is  necessary.  In  reference  7.1  it  is 
stated  that  "they  may  amount  to  little  more  than  the  creation  of  a  highly  idealised  flows,  some  of  whose  features 
coincide  with  the  corresponding  features  of  the  real  flow.  The  extent  of  the  approximations  inherent  in  such  models 
might  only  be  revealed  by  the  comparison  with  experiment  of  other  features  of  the  flow,  such  as  the  pressure 
distribution  over  the  downstream  surface  of  the  bluff  body.  It  is  at  this  point  that  we  are  handicapped  by  the  fact  that 
experimental  techniques  are,  at  this  moment,  lagging  behind  the  advance  of  theory”. 

7.1  RELEVANCE  OF  POTENTIAL  FLUID  FLOW  SOLUTIONS  TO  PARACHUTE 
AERODYNAMICS 

In  steady  potential  fluid  flow  all  or  a  part  of  the  immersed  body's  impervious  boundary  is  considered  to  be  one  of  a 
family  of  streamlines  which  represent  the  flow.  Since  these  streamlines  and  their  associated  equipotential  lines  form 
an  orthogonal  net,  in  a  fluid  which  is  wholly  irrotational  the  flow  pattern  which  is  developed  around  a  symmetrically- 
shaped  body  is  itself  symmetrical.  This  leads  to  a  symmetrical  pressure  distribution  and  in  consequence  in  steady 
flows,  to  zero  normal  pressure  drag,  otform  drag,  developing  on  the  immersed  body.  By  considering  the  momentum 
of  the  entire  flow  field  surrounding  the  immersed  body  it  can  readily  be  shown  that,  provided  that  the  body’s 
dimension ;  are  finite  and  it  is  considered  to  be  immersed  in  a  steady,  frictionless  incompressible  fluid  which  is 
entirely  free  of  vorticity,  then  regardless  of  its  shape  no  net  aerodynamic  force  can  be  developed  on  that  immersed 
body.  This  is  a  statement  of  the  D'Alembert  Paradox. 

In  view  of  this  paradox  and  considering  the  aerodynamic  characteristics  which  are  required  from  an  analytical  model 
of  a  parachute,  the  apparent  lack  of  relevance  of  potential  flow  solutions  to  parachute  aerodynamic  problems  must 
certainly  be  considemL  In  spite  of  what  D’Alembert  declared  to  be  a  paradox,  lift  is  certainly  developed  on  some 
bodies  which  are  immersed  in  potential  flows  and  similarly,  trail ing-vortex  drag  component  forces  can  be  generated  on 
models  which  have  a  finite  span.  However,  the  major  part  of  the  drag  developed  on  parachute  canopies  is  form  drag 
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and  if  they  are  immersed  in  steady  potential  fluid  flows  which  are  wholly  irrotational,  it  is  not  possible  to  determine 
this  form  drag  from  solutions  to  potential  flow  problems.  This  is  the  evident  disadvantage  of  the  first  approach. 

Through  free-streamline  theory,  originating  with  Helmholtz7'1  and  Kirchhoff71,  which  postulated  that  the 
pressure  is  constant  on  streamlines  which  extend  to  infinity  from  the  bluff  body  and  which  bound  the  wake  region,  a 
non-zero  drag  force  on  such  an  immersed  bluff  body  can  be  established.  Following  Fage  and  Johansen’s7  4 
assumption  that  the  fluid  pressure  varies  along  these  boundary  streamlines,  a  plausible  value  for  the  bluff  body  drag 
can  be  obtained.  Although  in  free-streamline  theory  the  assumption  is  made  that  flows  are  steady,  this  approach 
serves  as  a  useful  introduction  to  the  vortex-sheet  methods  of  analysing  unsteady  flows  around  bluff  bodies,  described 
below  in  Section  7.3. 

7.2  THF  IRROTATIONAL  FLOW  FIELD  APPROACH 

7.2.1  Ibrahim’s  Solution  for  the  Added  Mass  of  Fully-Inflated  Parachute  Canopies 

By  idealising  their  shapes  into  thin-walled,  rigid  cup  shapes  and  applying  conformal  transformation  techniques 
Ibrahim44,  in  his  1963  doctoral  thesis,  developed  a  potential  solution  for  flow  around  fully-inflated  parachute 
canopies.  His  reason  for  making  this  analysis  was  in  order  to  determine  the  added  mass  coefficients  associated  with 
the  unsteady  oscillatory  motion  of  the  parachute  during  its  descent.  Lester4-15,  whose  work  was  broadly 
contemporary  with  that  of  Ibrahim,  noted  that  “the  theoretical  concept  of  added  mass  with  regard  to  motion  of  bodies 
in  an  ideal  fluid  is  not  necessarily  representative  of  the  physical  phenomena  which  occur  in  a  real  fluid”.  Later 
experimental  measurements  tend  to  justify  this  earlier  opinion. 

7.2.2  Klimas’  Parachute  Canopy  Method 

Only  certain  representative  shapes  of  immersed  bodies  respond  well  to  the  methods  of  conformal  transformation 
and  a  more  flexible  technique  is  required  with  which  to  represent  axi-symmetrical  canopies  with  arbitrary  cross- 
sectional  shapes.  Milne-Thomson4-5  modelled  two-dimensional  thin  aerofoil  sections,  representing  them  by  a  line 
vortex  sheet  and  in  1972  Klimas5-15  followed  his  example,  modelling  an  axi-symmetric  parachute  by  a  system  of 
vortex  rings  which  covered  the  canopy.  The  modelling  of  axi-symmetric  shapes  by  vortex  rings  is  discussed  in 
reference  7.1.  This  representation  enabled  Klimas  to  include  the  effects  of  canopy  porosity  in  his  model. 

Klimas’  initial  objective  was  that  of  determining  the  pressure  field  round  an  inflating  canopy.  In  seeking  to  meet 
|  this  aim  he  needed  good  experimental  data  with  which  to  validate  his  model  and  this  was  not  easy  to  find.  In  1977, 

like  Ibrahim,  Klimas4  7  used  a  development  of  his  earlier  mode!  as  the  basis  for  determining  the  added  mass 
coefficients  in  unsteady  motion,  remarking  that  “no  obstacle  exists  to  extension  of  the  approach  to  include  the  non- 
stcady  canopy  geometries  (ot  the  inflation  process)”.  In  a  third  paper  published  in  1979  Klimas514  adopted  his 
earlier  vortex  sheet  canopy  model  as  the  means  by  which  he  determined  the  pressure  field  round  an  inflating  parachute 
|  canopy.  The  shapes  adopted  by  the  inflating  canopy  were  assumed  and  the  process  of  inflation  was  permitted  to 

continue  until  the  axial  aerodynamic  force,  evaluated  by  integration  of  the  pressure  distribution,  was  equal  to  a  value 
which  had  been  indepently  determined  or  had  been  assumed. 

An  objection,  ascribed  to  Roberts  and  Reddy55  in  Section  3.3,  to  this  technique  which  Klimas  adopted  is  the 
necessity  for  the  inflating  canopy  shape  to  be  an  input  to  the  determination  of  the  pressure  field,  rather  than  the 
canopy  shape  developing  as  a  significant  output  to  the  analytical  procedure  which  is  adopted. 

7.2.3  Roberts’  Inflating  Canopy  Method 

In  order  to  determine  inflation  times  and  inflation  loads  for  parachutes,  Roberts511  determined  the  unsteady 
pressure  distribution  over  an  axisymmetric,  impervious,  inflating  shell.  In  order  that  its  potential  flow  pressure 
distribution  could  be  developed  by  conformal  transformation  techniques  from  that  in  a  right-angled  comer  a  paraboloid 
was  chosen  as  the  general  shape  of  this  shell.  At  any  instant  in  time  the  dynamics  of  the  canopy  inflation  process 
could  be  determined  from  the  payload  mass  and  knowledge  of  the  drag  force  developed  on  the  canopy,  the  latter  being 
given  immediately  from  the  known  pressure  distribution  which  developed  over  it 

Its  aerodynamic  analysis  is  developed  by  considering  the  existence  of  a  starting  vortex  ring  which  forms  and  grows 
in  a  location  adjacent  to  the  skirt  of  the  shell.  It  is  through  the  existence  of  this  vortex  ring  that  the  Kutta-Joukowski 
condition  is  satisfied  at  the  shell  skirt.  From  observations  of  parachute  characteristics  it  is  evident  that,  as  the  shell 
becomes  nearly  fully  inflated,  the  vortex  ring  at  the  shell  skirt  become  unstable,  drifting  downstream  from  the  shell 
and  forming  the  axisymmetric  wake  in  the  canopy.  In  his  model  Roberts  showed  how  the  position  of  this  vortex  ring 
could  be  determined. 

In  Roberts’  method,  unlike  the  one  which  Klimas  adopted,  the  particular  shape  of  the  inflating  canopy  followed 
from  the  determination  of  the  pressure  field  which  developed  around  if.  Essentially,  it  is  a  vortex  sheet  method  of 
calculation  applied  to  a  canopy  shape  which  varies  with  time. 

7.3  VORTEX  SHEET  METHODS  OF  REPRESENTING  THE  WAKES  SHED  BV 
PARACHUTES 

The  representation  within  an  otherwise  irrotational  fluid  flow  of  an  identifiable  region  of  vorticity  has  a  long 
history.  This  is  the  method  which  Prandtl7-5  adopted  to  model  the  boundary  layer.  Lanchestcr7*  similarly 
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modelled  the  vortex  system  v'hich  forms  behind  an  aircraft  wing  and  von  Kirmdn7,7  used  vortex  sheet  methods  to 
model  the  wake  established  behind  a  bluff  body .. 

The  first  representation  by  discrete  vortices  of  a  two-dimensional  continuous  vortex  sheet  was  made  by 
Roscnhead™.  He  showed  that  this  vortex  sheet  was  unstable,  rolling  up  smoothly  into  concentrated  clusters  of 
vortices.  An  extensive  review  of  more  recent  two-dimensional  vortex  sheet  representations  has  been  made  by 
Clements  and  Maull”. 

The  simulation  by  vortex  rings  of  the  wakes  shed  by  axi-symmetrical  bodies  poses  an  additional  problem  to  the 
two-dimensional  representation;  that  of  the  self-induced  velocity  field  induced  by  the  curvature  of  the  vortices.  In  order 
to  model  the  asymmetry  of  the  real  fluid  flow  around  a  symmetrical  parachute  canopy,  Shirayama  and 
Kuwahara71*71'  have  simulated  the  wake  by  using  segmented  vortex  rings,  thereby  forming  a  series  of  vortex 
sticks.  In  a  recent  review  paper  Strickland711  has  described  a  number  of  contemporary  applications  of  vortex-sheet 
analysis  methods  for  parachute  canopies. 

By  using  vortex-sheet  methods  it  is  possible  to  determine  the  pressure  distribution  over  an  immersed  bluff  bodv 
under  any  known  instantaneous  velocity  and  acceleration.  This  can  then  be  integrated  to  obtain  the  component 
aerodynamic  forces  and  moments  developed  in  any  required  direction.  Work  by  Rehbach7 11  has  shown  that  by  such 
a  method  normal  force  coefficients  developed  on  inclined  flat  plates  can  be  determined. 

Canopy  porosity  can  be  simulated  in  any  panel  of  the  bound  vortex  rings.  Although  compressibility  effects  have 
been  included  in  steady  flow  vortex  methods  Strickland  remarked  that  at  this  instant  of  time  their  inclusion  into  the 
unsteady  method  necessary  to  simulate  the  a  bluff  body  wake  appeared  to  be  premature. 

7.3.1  Vortex  Sheet  Methods  Applied  to  Bodies  of  Various  Geometries  by  Meyer  and  Purvis 

Preliminary  results  from  a  model  capable  of  predicting  unsteady,  viscous,  incompressible  flows  about  parachute 

canopies  have  been  published  by  Meyer  and  Purvis714.  At  the  time  at  which  they  made  this  presentation  the  results 
from  their  model  were  limited  to  the  two-dimensional  pressure  distribution  around  a  circular  cylinder,  which  they 
compared  with  experimental  results.  Ultimately,  they  plan  to  extend  their  technique  so  that  it  can  be  applied  to 
predictions  of  the  pressure  distribution  associated  with  an  inflating  axi-symmetric  parachute  canopy. 

7.3.2  Flow  Around  Discs  by  de  Bernardlnis,  Graham  aud  Parker 

Assuming  potential  fluid  flow  around  a  disc  and  representing  that  disc  by  an  appropriate  distribution  of  bound 
vortex  rings  while  the  vortex  sheets  that  it  sheds  are  simulated  by  discrete  vortex  rings,  de  Bemardinis,  Graham  and 
Parke711  obtained  a  soluuon  to  flow  behaviour  in  the  neighbourhood  of  a  sintisoidally-oscillaring  disc.  The 
method  which  they  adopted  is  essentially  similar  to  that  used  by  Roberts  in  reference  3.11.  The  Kutta-Joukowski 
condition  was  applied  at  the  region  of  flow  separation  from  the  disc.  From  this  region  a  free  vortex  sheet,  represented 
by  a  number  of  vortex  rings,  was  converted  with  die  local  fluid  velocity. 

The  flow  was  determined  at  a  number  of  different  Keulegan -Carpenter  numbers,  K  =  V^T/R,  where  V„,  is 
the  amplitude  of  velocity  oscillation  for  the  disc,  R  is  the  disc  radius  and  T  the  period  of  oscillation.  These 
calculations  were  found  to  be  in  reasonably  good  agreement  with  results  from  flow  visualisation  experiments.  The 
unsteady  pressure  field  acting  on  the  disc  was  also  determined.  It  was  concluded  that  the  vortex-sheet  method  which 
they  adopted  accurately  predicted  the  dominant  features  of  the  flow  around  the  oscillating  disc. 

7.3.3  McCoy  and  Wcrme’s  Axisymmetric  Vortex  Lattice  Method  Applied  to  Parachute  Shapes 

Applications  to  axisymmetric  parachute  canopy  shapes  are  shown  in  the  paper  describing  McCoy  and 

Werme’s71*  vortex  lattice  method,  still  in  its  early  stages  of  development  By  assuming  that  the  velocity- 
dependent  drag  coefficient  was  equal  to  its  steady-state  value  in  a  decelerating  flow  and  considering  the  remaining  part 
of  the  total  aerodynamic  force  to  be  acceleration-dependent,  the  instantaneous  added  mass  coefficients  described  in 
Section  4.2  were  calculated.  McCoy  and  Werme  found  that  although  the  initial  value  of  these  coefficients  compared 
well  with  an  inviscid  flow  calculation  made  in  the  absence  of  vortex  sheets,  over  a  period  of  time  the  added  mass 
coefficients  rapidly  decreased  in  value,  eventually  becoming  negative. 

7.4  SUMMARY  OF  REQUIREMENTS  FOR  THE  DETERMINATION  OF  PARACHUTE 
AERODYNAMIC  CHARACTERISTICS  BY  VORTEX  SHEET  METHODS 

Currently,  analytical  methods  are  required  by  which  the  aerodynamic  characteristics  of  parachute  canopies  can  be 
determined  under  the  following  flow  conditions: 

(0.  in  steady  flows. 

The  current  state  of  the  an  is  that  parachute  drag  coefficients  appear  to  be  calculable  by  vortex-sheet 
methods  and  it  is  probable  that  by  similar  techniques  both  normal  and  axial  aerodynamic  force  , 

coefficients  could  be  calculated.  Although  determinations  have  been  made  by  Rehbach711  on  ; 

rectangular  plates,  no  publications  are  known  in  which  this  has  been  shown  for  axisymmetric  bodies  | 

like  parachute  canopies,  set  at  known  angles  of  attack  in  steady  flows.  | 


62 


(ii).  in  unsteady  flows,  for  fully-inflated  parachute  canopies. 

The  problem  here  is  to  determine  the  added  mass  coefficients  analytically.  Ibrahim44  and  others  have 
achieved  this  for  wholly  itrotational  flows,  although  Lester4,45  has  questioned  the  engineering 
significance  of  the  results  which  they  obtained.  As  McCoy  and  Wem;'  •*  have  indicated  such 
determinations  can  be  made  by  vortex-sheet  methods  but  as  yet  it  is  still  too  early  for  much  published 
data  to  be  available. 

(iii).  in  unsteady  flows,  for  inflating  parachute  canopies. 

Here,  reliable  methods  are  still  required  with  which  to  couple  the  aerodynamic  forces  obtained  by 
vortex-sheet  methods  to  a  structural  analysis  technique.  Roberts511  has  indicated  that  such  a 
coupling  is  possible:  Meyer  and  Purvis7,14  are  attempting  its  solution,  but,  as  yet,  there  are  no 
indications  that  a  solution  has  been  obtained  to  this  problem. 

REFERENCES 

7. 1  Thwaites,  B.  (Editor).  Incompressible  Aerodynamics.  Clarendon  Press,  Oxford.  I960 

7.2  Helmholtz,  H.von.  Obcr  Discontinuierliche  FKlssigkeitsbewegungen.  ‘Mber’  K.  Akad.  Wiss.,  Berlin,  22, 

215,  1868. 

7.3  Kirchhoff,  G.  Zur  Theorie  Freier  Fllissigkeitsstrahlen.  J.  Heine  Ajigew.  Math.  2Q,  289, 1869. 

7.4  Page,  A.  and  Johansen,  F.C.  On  the  Flow  of  Air  Behind  an  Inclined  Flat  Plate  of  Infinite  Span.  U.K. 

Aeronautical  Research  Council  R.  &  M.  2127, 1927 

7.5  Prandtl,  L.  Ober  Fliissigkeitsbewegung  bei  sehr  kleiner  Reibung.  Verh.  3.  int.  Math.-Kongr.,  Heidelburg 

484,  1904. 

7.6  Lanchestcrp.W.  Aerodonetics.  Constable,  London.  1908. 

7.7  Kdrmfa,  T.von  FIfissigkeits-  und  Luftwiderstand.  Fhys.Z  12, 49, 191 1. 

7.8  RosenheadJ-.  Formation  of  vortices  from  a  surface  of  discontinuity.  Proceedings  of  the  Royal  Society.  A, 

134.  pp.170-192, 1931. 

7.9  Clements,  R.R.  and  Maull,  DJ.  The  Representation  of  Sheets  of  Vorticity  by  Discrete  Vortices.  Progress 

of  Aerospace  Science,  ifi,  2,  pp.129-146,  Pergamon  Press,  1975. 

7.10  Shirayama,  S.  and  Kuwahara,K.  Computation  of  Flow  Past  a  Parachute  by  a  Three-Dimensional  Vortex 

Method. 

AIAA-CP-86-0350.  Proceedings  of  the  AIAA  24th  Aerospace  Sciences  Meeting,  Reno,  January  1986. 

7. 1 1  Shirayama,  S.  and  Kuwahara  JC.  Superconductors  and  Fluid  Dynamics.  Proceedings  of  the  First  Nobcyama 

Workshop, 1985,  Edited  by  Kuwahara,  Mendez  and  Orszag.  Springer- Verlag,  1985. 

7.12  StricklandJ.H.  On  the  Utilization  cf  Vortex  Methods  for  Parachute  Aerodynamic  Predictions.  AIAA-CP-86- 

2455.  Proceedings  of  the  AIAA  9th  Aerodynamic  Decelerator  and  Balloon  Technology  Conference, 
Albuquerque.  October  1986. 

7.13  Rehbach.C.  Numerical  Calculation  of  Three-Dimensional  Unsteady  Flows  with  Vottex  Sheets.  AIAA-CP- 

78-111.  Proceedings  of  the  16th  Aerospace  Sciences  Meeting,  Huntsville,  January  1978. 

7.14  Meyer,  J.  and  Purvis,  J.W.  Vortex  Lattice  Theory  Applied  to  Parachute  Canopy  Configurations.  AIAA-84- 

0795.  Proceedings  of  the  AIAA  8th  Aerodynamic  Decelerator  and  Bal'oon  Technology  Conference, 
HyannisApril  1984. 

7.15  deBemardinis.B.,  GrahamJ.MJL  end  Parker JCH.  Oscillatory  Flow  Around  Disks  and  Through  Orifices. 

Journal  of  Fluid  Mechanics,  102.  po.279-299. 1981. 

7.16  McCoy JH.H.  and  Wcrme.T.D.  Axisymmetric  Vortex  Lattice  Method  Applied  to  Parachute  Shapes.  AIAA- 

CP-86-2456.  Proceedings  of  the  AIAA  9th  Aerodynamic  Decelerator  and  Balloon  Technology 
Conference,  Albuquerque,  Ocobcr  1986. 


t 

I 


8.  EXTRA-TERRESTRIAL  APPLICATIONS  OF  PARACHUTES 


During  the  ten  years  which  began  in  1976  intensive  activities  in  both  the  United  States  of  America  and  the  Soviet 
Union  have  been'directed  towards  exploring  the  characteristics  of  three  other  planets  in  the  solar  system,  Mars,  Venus 
and  Jupiter.  In  separate  missions  during  July  and  August  1976  two  Viking  instrument  packages  launched  from  the 
United  States  were  soft-landed  on  to  the  surface  of  the  planet  Mars.  Parachutes  were  used  on  both  of  these  operations 
as  a  significant  part  of  their  deceleration  phase.  The  entry  capsules  entered  the  Martian  atmosphere  at  a  Mach  number 
of  about  2.0,  at  which  stage  the  parachutes  were  deployed  in  order  to  decelerate  the  system.  They  were  then 
jettisoned  as  retiopropulsion  was  responsible  for  the  payloads'  terminal  descent  As  the  atmospheric  characteristics 
on  the  planet  Mars  might  well  support  life  as  we  understand  it  it  was  important  to  sterilise  the  entire  descent 
systems,  including  the  parachute  canopies.  This  was  done  for  200  hours  at  a  temperature  of  about  140  deg.C.  (280 
deg.F). 

A  mission  whose  purpose  is  to  bring  back  to  earth  a  selected  S  kg.  set  of  sample  materials  from  the  surface  of 
Mars  had  been  planned  by  the  United  States.  This  was  to  have  taken  piace  in  the  late  1980’s,  but  budgetary  and  other 
considerations  have  delayed  this  intended  mission  until  the  1990's. 

In  September  1977  two  United  States  Pioneer  probes  began  their  journey  to  the  planet  Venus,  having  been 
scheduled  to  arrive  there  in  December  1978.  For  partial  trajectory  control  over  altitudes  from  67  -  44  km  in  the 
Venus  atmosphere  the  largest  of  these  two  probes  was  designed  to  employ  a  parachute.  Though  this  mission  was 
sufficiently  successful  for  some  data  to  be  acquired,  it  is  not  known  by  the  writer  to  what  extent  this  larger  probe  and 
its  parachute  deederator  materially  contributed  to  this  success. 

Reference  8.1  radicals  that  in  March  1982  two  Soviet  probes  which  had  been  launched  in  1981  safely  landed  on 
Venus.  These  probes,  which  were  called  Veneras  13  and  14,  had  parachutes  which  were  jettisoned  about  45  km  (28 
miles)  above  that  planet's  surface.  Since  the  Venus  atmosphere  is  hostile,  having  an  atmospheric  pressure  some  90 
times  that  on  e=vtlt  and  an  atmospheric  temperature  which  is  around  450  deg  C.  (850  deg.  F.)  the  surface  survival 
time  would  be  small.  However,  it  was  long  enough  for  chemical  analyses  to  be  made  of  the  soil  and  for  landing  sites 
to  be  photographed,  the  relevant  data  being  radioed  back  to  the  earth.  By  late  1982  a  total  of  seven  successful  Soviet 
landings  on  Venus  had  been  reported  but  it  is  not  known  for  how  many  of  these  missions  parachutes  were  employed 
to  decelerate  the  probes. 

A  Galileo  spacecraft  designed  to  enter  the  atmosphere  around  the  planet  Jupiter  in  August  1988  was  also  planned. 
Although  scheduled  for  an  American  launch  in  May  1986  this  research  programme  has  also  suffered  considerable 
delays. 

Over  this  ten  year  period  there  has  been  a  very  considerable  increase  in  knowledge  concerning  the  physical 
characteristics  of  the  planets  which  were  being  explored.  Many  writers  report  that  at  the  time  the  Venus  probes  were 
launched  the  atmosphere  around  that  planet  was  much  better  understood  than  was  the  atmosphere  around  Mars  in  the 
late  1970’s. 

For  ail  of  these  planets  the  adopted  parachute  designs  closely  resemble  those  which  have  already  been  used  for 
spacecraft  recovery  in  the  earth's  atmosphere. 

8.1  ATMOSPHERIC  CHARACTERISTICS  ON  MARS,  VENUS  AND  JUPITER 

According  to  reference  1.9,  little  is  known  with  any  certitude  about  the  characteristics  of  the  atmosphere  around 
Mars.  At  the  surface  it  is  certainly  very  cold  and  it  has  only  about  0.01  times  the  density  of  the  earth’s  atmosphere. 
The  pressure  there  has  been  obtained  from  die  Mariner  space  research  programme.  Darnell,  Henning  and 
Lundstrom*-2  report  data  front  Mariner  IV  which  gives  pressure  at  the  surface  of  Mars  as  probably  between  5  and  10 
millibars.  Moog.  Bendura,  Timmons  and  Lau,J  add  that  even  with  the  data  from  the  Mariner  9  mission  large 
uncertainties  still  remain  about  the  Martian  atmospheric  densities  and  scale  heights.  For  landings  on  Mars  high  entry 
velocities  are  necessary,  i.’ference  8.2  quoting  3.7  to  4.9  m/sec.  (presumably  what  was  intended  was  3.7  to  4.9 
km/sec.)  (12  000  to  16  000  ft/sec.)  ?.t  altitudes  of  4.6  to  6.1  km  (15  000  to  20  000ft.).  The  corresponding  entry 
Mach  Number  was  believed  to  be  a  liule  above  1.0.  Reference  83  states  that  the  parachute  must  be  capable  of 
operating  over  a  Mach  number  range  from  as  high  as  2.0  to  a  low  subsonic  Mach  number  and  perform  without 
damage  in  a  range  of  dynamic  pressures  from  24  -  479  N/m1  (03  to  10  Ibf/ft*).  In  this  reference  an  entry  dynamic 
pressure  of  239  -  31 1  N/mJ  (5.0  to  6.5  Ibf/ft1)  is  claimed.  In  a  much  earlier  paper  Heinrich*4,  quoting  NASA*3 
and  General  Electric  Company*4  data,  estimated  the  first  stage  of  entry  into  the  Martian  atmosphere  to  be  at  a 
velocity  of  0.93  km/sec.  (3  C56  ft/sec.)  at  an  altitude  of  7.35  km.  (24  000ft)  and  a  Mach  number  of  4  0.  Even 
allowing  for  the  undoubted  rapid  refinement  of  data  on  the  atmosphere  around  Mars  it  is  clear  that  there  are  very  wide 
uncertainty  bands  present. 

On  the  planet  Venus  reference  1.9  indicates  that  the  atmospheric  conditions  are  reasonably  well  established.  The 
atmospheric  dcruHy  is  about  100  times  that  of  the  earth’s  atmosphere  and  the  surface  temperature  is  dose  to  480 
deg.C  (900  deg.  F.)  Because  of  this  hostile  environment,  life  was  not  considered  to  be  possible  and  thus  a 
biologically  clean  exploratory  system  >ai  not  considered  essential. 
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A  source  of  relevant  ~nd  up-to-date  information  about  the  properties  cf  the  atmosphere  around  these  three  planets 
is  the  Journal  of  Geophysical  Research.  For  example,  data  on  the  Martian  atmosp'iere  are  contained  in  references  8.7 
and  8.8.  Reference  8.9  is  also  to  a  model  of  the  atmosphere  around  Mars,  while  reference  8.10  is  to  the  atmosphere 
around  Jupiter. 

8.2  MISSION  REQUIREMENTS  FOR  PARACHUTES 

On  the  Viking  Mars  mission  reference  1.9  explains  that  the  capsules  were  designed  to  enter  the  Martian 
atmosphere  at  about  245  km  (800  000  ft).  At  about  6 5  km  (21 000  ft),  above  the  surface  and  at  a  velocity  of  about 
365  m/sec.  (1 200  ft  /sec.)  a  disc-gap-band  parachute  opened.  This  parachute  was  disconnected  at  an  altitude  of  about 

1.2  km  (4  000  ft.)  when  the  payload  velocity  was  about  60  m/sec.  (200  ft  /sec.). 

On  the  Pioneer  Venus  mission,  reference  1.9  states  that  the  planet’s  atmosphere  was  entered  at  about  67  km  (220 
000  ft),  the  300  kg  (670  lb)  probe  decelerating  to  a  Mach  number  of  about  0.8.  At  this  altitude  the  dynamic  pressure 
was  3  300  S/m1  (69  Ibf/ft3).  Using  a  guide  surface  pilot  parachute  the  main  conical  ribbon  parachute  was  then 
deployed.  The  prime  function  of  this  5  m  (16.2  ft)  nominal  diameter  canopy,  which  had  a  drag  coefficient  of  0.52, 
was  to  stabilise  the  probe  through  the  Venus  cloud  cover  so  that  scientific  examination  of  the  atmosphere  could  be 
carried  out  By  47  km  (155  000  ft)  the  velocity  of  the  parachute  and  its  payload  was  so  low  that  the  parachute  was 
jettisoned,  the  impact  of  the  probe  on  the  surface  of  Venus  occurring  37  minutes  later.  The  requirement  for  this 
parachute  was  thus  much  more  limited  than  was  that  for  the  Viking  mission  to  Mars.  However,  as  the  density  of  the 
atmosphere  on  Venus  is  large  and  the  atmospheric  temperature  close  to  the  planet  very  high,  at  a  lower  altitude  than 
47  km  a  parachute  is  neither  needed  nor  would  it  have  been  a  practicable  proposition. 

Corridan,  Givens  and  Kepleyui  indicate  that  the  purpose  of  the  Galileo  mission  is  to  explore  the  planet  Jupiter 
and  its  satellites  by  indirect  measurements,  made  from  an  orbiting  vehicle,  as  well  as  more  direct  by  atmospheric 
measurements,  made  from  an  entry  probe.  This  probe  is  designed  to  enter  the  atmosphere  of  Jupiter  at  47  km/sec. 
(29  miles/sec.),  then  to  be  slowed  by  its  blunt  forebody  to  a  transonic  velocity.  At  a  Mach  number  of  between  0.91 
and  1.01  and  a  corresponding  dynamic  pressure  of  between  4  850  and  7  650  N/to2  (102  and  160  lbf/ft1)  a  20  deg. 
conical  ribbon  pilot  parachute  of  1.1  m  (3.74  fi)  nominal  diameter  is  to  deploy.  When  the  entry  Mach  number 
decreases  to  between  0.87  to  0.97  the  conical  ribbon  main  parachute  of  3.8  m  (12.48  ft)  nominal  diameter  is  to  be 
deployed.  The  purposes  of  this  latter  are  to  separate  the  instrumented  descent  module  fom  its  heat  shield  and  to 
provide  drag  for  a  controlled  descent  through  the  atmosphere.  Further  details  of  the  Galileo  mission  to  the  planet 
Jupiter  are  given  in  references  8.12  and  8.13. 

8.3  HEINRICH’S  1966  ANALYSIS  OF  EXTRA-TERRESTRIAL  PARACHUTE 
AERODYNAMICS 

In  1966,  when  considering  the  behaviour  of  parachutes  descending  in  a  Martian  environment,  Heinrich'4 
stated  that  the  most  significant  performance  characteristics  of  a  parachute  system  were: 

(i) .  its  rate  of  descent; 

(ii) .  its  dynamic  stability  characteristics  in  pitch; 

(iii) .  its  opening  time  and 

(iv) .  its  opening  shock  load. 

These  characteristics  depend  on  the  weight  cf  the  payload  and  on  the  parachutes's  aerodynamic  characteristics. 
These  latter  are  functions  of  the  canopy  shape  and  attitude,  its  Reynolds  number  and  its  Mach  number.  Heinrich 
argued  that  the  most  significant  way  in  which  the  parachute  canopy  shape  would  be  altered  when  it  descended  through 
the  atmospheric  environment  of  a  planet  other  than  the  Earth  is  through  the  changes  which  would  occur  in  the 
effective  porosity  of  the  canopy. 

Working  from  earlier  research113  which  he  had  undertaken  on  the  porosity  of  parachute  canopies,  Heinrich 
argued  that  at  a  given  dimensionless  pressure  ratio  across  the  canopy  the  variation  in  the  effective  porosity  of  the 
canopy  from  that  at  sea  level  in  the  earth’s  atmosphere  could  be  expressed  as  a  function  of  two  variables,  the  ratio  of 
the  density  of  the  fluid  in  which  the  canopy  is  immersed  to  that  of  air  at  sea  level  and  the  Reynolds  number  at  which 
the  parachute  descends.  Heinrich  defined  the  dimensionless  pressure  ratio  across  the  canopy  fitbric  as  the  ratio  of  the 
actual  pressure  difference  across  the  canopy  to  that  which  would  establish  sonic  flow  in  the  fabric  interstices.  Because 
the  fluid  density  in  the  Martian  environment  is  much  lower  than  it  is  above  the  earth’s  surface,  he  established  that  the 
effective  porosity  of  a  parachute  canopy  descending  on  the  planet  Mars  would  be  only  between  15%  and  30%  of  its 
value  when  descending  above  the  eanh’s  surface. 

Fbr  the  type  of  parachute  canopy  envisaged  for  high  Mach  number  inflation  in  the  Martian  environment  he  argued 
that  the  variation  of  its  aerodynamic  characteristics  with  Reynolds  number  would  be  weak  and  that  the  Mach  number 
was  the  dominating  dimensioniess  parameter.  He  therefore  recommended  that  in  model  rests  the  Mach  number  should 
be  made  equal  to  that  of  the  prototype,  defining  the  modelling  conditions  as  follows:- 
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(i) .  the  model  Mach  number  should  bethct  of  the  prototype; 

(ii) .  the  fluid  density  for  the  model  tests  should  be  that  of  the  fluid  in  which  the  prototype  is  immersed; 

(iii) .  the  dimensionless  pressure  ratio,  defined  above,  across  the  canopy  fabric  for  the  model  tests  should  be  the 

same  as  that  for  the  prototype. 

By  establishing  these  three  modelling  conditions,  Heinrich  sought  to  establish  values  of  aerodynamic  coefficients 
on  model  canopies  in  the  earth's  atmosphere  which  would  be  identical  to  those  on  full-scale  prototype  canopies  in  the 
Martian  atmosphere.  Then,  by  making  the  fluid  density  for  the  model  canopy  equal  to  that  for  the  full-scale 
prototype,  the  rates  of  descent  VD  could  be  related  through  the  respective  sizes  of  canopy,  the  known  gravitational 
accelerations  in  both  environments  and  the  known  masses  of  the  payloads  which  would  be  attached  to  the  canopies. 

Next,  he  considered  a  simplified  expression  for  the  motion  of  a  parachute  and  it:  payload  when  oscillating  in 
pitch.  Using  a  relationship  which  was  similar  in  form  to  equation  3.20,  but  taking  the  origin  of  the  parachute  and 
payload  system  to  be  at  its  centroid;  by  neglecting  the  gravitational  moment  of  the  canopy  about  the  origin  compared 
with  that  of  the  payload  at  bD  from  the  centroid,  he  obtained: 

mb’D’q  =  '4pVB*(*/4)  D*C.  (8.1) 

and  substituting  for  the  descent  velocity  V0 

VJ  =  ■*  /  I'/jPCM  D*CdJ  (8.2) 

produced: 

q  =  (C.  /C#)(g/b*D)  .  (8.3) 

With  identical  aerodynamic  characteristics  in  the  two  environments  Heinrich  thus  established  that  for  q  to  be  the 
same  in  each,  then: 


(8/D)t.r»  =  (g/D)«-  -  (8.4) 

Since  the  ratio  of  the  gravitational  accelerations,  g^/gem*  ~  0.38,  this  established  that 

IWDbrtk  =  0-38  .  (8.5) 

Thus,  the  size  of  the  model  would  need  to  be  2,63  times  the  size  of  the  prototype.  Similarity  in  the  rate  of 
descent  requirement,  equation  8 2,  then  led  to  the  ratio  of  payload  masses: 

*nMu»/Bllara  -  0.38  •  (8,6) 

Thus,  in  order  to  achieve  the  same  asrodynamic  coefficients  in  the  two  different  environments  Heinrich  sought  to 
make  tests  on  a  model  which  was  of  the  same  shape  as  the  prototype.  To  ensure  the  shapes  were  wholly  identical  he 
considered  how  to  make  the  canopy  porosity  on  the  Earth  similar  to  that  on  Mars.  Although  Heinrich’s  model 
canopy  was  tested  at  full-scale  Mach  number  its  Reynolds  number  was  dl  imilar  to  that  of  the  prototype. 

If  model  aerodynamic  tests  in  the  earth's  atmosphere  are  performed  under  conditions  in  which  the  air  density  is 
equal  to  the  fluid  density  in  the  atmosphere  around  Mars,  equation  8.2  then  indicates  that  the  rate  of  descent  for  the 
model  will  be  equal  that  d  the  prototype.  This  entails  testing  the  model  canopy  in  the  earth’s  atmosphere  at  an 
altitude  of  about  32  km  (105  000  ft)  above  the  earth’s  surface.  Under  these  circumstances,  in  the  earth's  atmosphere 
the  Proude  number  F,  shown  in  equation  2.39  to  equal  V£/D.g,  will  be  about  one-seventh  of  its  value  for  the 
prototype  descending  in  the  Martian  atmosphere. 

Although  the  model  test  in  the  earth's  atmosphere  may  meet  all  of  Heinrich’s  specifications,  because  of 
differences  m  Reynolds  number  and  of  canopy  porosity  in  there  tests  there  is  no  possibility  of  ensuring  total  dynamic 
similarity.  However,  the  prototype  parachute  in  the  Martian  environment  has  to  be  deployed  at  high  Mach  numbers 
and  in  order  to  meet  this  requirement  highly  porous  ribbon  canopies  are  required.  For  these  there  will  be  minimal 
variation  of  aerodynamic  characteristics  with  Reynolds  number.  Since  a  given  parachute  canopy  is  less  porous  in  the 
Martian  atmosphere  than  it  is  in  the  atmosphere  surrounding  the  earth  and  as  has  been  discussed  in  Section  22.1.4, 
reducing  die  porosity  decreases  the  canopy  stability  in  pitch,  it  is  necessary  to  test  in  the  earth’s  atmosphere,  but  at 
altitudes  of  some  32  km  (105  000  ft)  where  the  atmospheric  density  is  of  the  order  of  that  on  Mars,  parachutes 
intended  for  use  within  the  Martian  atmosphere.  If  under  these  conditions  the  ribbon  canopies  exhibit  the  required 
stability  characteristics  in  pitch,  then  in  the  Martian  environment  they  should  be  adequately  stable. 


8.4  TESTS  ON  PARACHUTES  FOR  EXTRA-TERRESTRIAL  APPLICATIONS  CONDUCTED 

IN  1HE  EARTH’S  ATMOSPHERE 

Reference  1.9  describes  the  three  experimental  programmes  which  were  devised  to  establish  the  characteristics  of 
the  parachutes  for  Viking  Man  landings.  Initially,  the ‘Planetary  Entry  Parachute  Programme  was  established  to 
select  the  most  appropriate  canopy  shape.  Disk-gap-band,  ringsail  and  cross  parachutes  were  deployed  at  altitudes  in 
excess  of  30  km  (100  000  ft)  and  at  Mach  numben  from  1.0  to  2.8  in  these  rocket  and  balloon-launched  tests.  Asa 
consequence  the  disk-gap-band  parachute  was  selected. 

In  the  subsequent  Low  Altitude  Drop  Test  Programme  parachute  opening  loads  and  stresses  at  up  to  1.5  times  the 
predicted  design  loads  were  investigated  as  the  canopies  were  deployed  from  a  B-57  aircraft  flying  at  IS  km  (SO  000  ft) 
altitude.  Finally,  in  order  to  test  the  transonic  interference  effects  produced  by  the  large  and  blunt  forebody  and  to 
check  the  stability  characteristics,  Moog,  Bendura,  Timmons  and  Lauu  report  that  balloon-launched  tests  were 
conducted  in  which  a  simulated  full- scale  Viking  vehicle  attached  to  a  16  metre  (S3  ft)  nominal  diameter  parachute 
was  tested  over  a  Mach  number  range  from  0.47  to  2.18. 

From  reference  1.9  it  is  seen  that  tests  in  the  earth’s  atmosphere  of  the  Pioneer- Venus  probe  vehicle  and  parachute 
occurred  in  two  stages.  In  the  first  of  these  a  bomb-shaped  test  vehicle  was  dropped  flora  an  F-4  aircraft  flying  at  12 
km  (40  000  ft)  altitude.  In  the  second  stage  a  simulated  probe  test  vehicle  with  a  mass  of  304  kg  (670  lb)  together 
with  its  entire  parachute  system  was  released  from  a  balloon  at  27.5  km  (90  000  ft)  and  the  parachute  was  deployed 
satisfactorily. 

Tests  in  the  earth’s  atmosphere  devised  to  simulate  the  Galileo-Jupiter  mission  arc  reported  by  Corridan,  Givens 
and  Kepley*-11.  Parachute  deployment  was  required  at  16.5  km  (54  391  ft),  a  dynamic  pressure  of  6  000  N/m1 
(125  lbflft*)  and  a  Mach  number  of  0.92.  At  transonic  speeds  problems  were  encountered  because  of  poor  parachute 
performance  in  the  wake  of  the  blunt  forebody  and  some  wind  runnel  deployment  tests  were  initiated  in  the  NASA 
Langley  transonic  dynamics  wind  tunneL  In  these  experiments  freon  12  used  as  the  test  medium  and  one-quader  and 
one-half  scale  models  of  the  pilot  parachute  were  tested.  Since  the  wind  tunnel  working  section  had  a  cross-sectional 
area  of  23  raJ  (248  ft1)  the  blockage  area  ratios,  described  in  Sectioa  6.6,  for  these  tests  were  a  negligible  0.4%  for 
the  one-quarter  scale  model  and  between  1%  and  3ft  for  the  one-half  scale  modeL 
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9.FURTHER  AERODYNAMIC  RESEARCH  INTO  PARACHUTES 


Knacke111  has  suggested  three  important  milestones  that  have  occurred  in  parachute  development: 

(i) .  Irvin’s  free-fall  parachute  jump  in  April  1919  as  a  member  of  Hoffman's  U.S.  Army  team.  After  freefalling 

for  a  short  distance  Irvin  pulled  a  ripcord  and  opened  his  parachute  pack  which  had  no  static  line 
connection  between  the  pack  and  the  aircraft; 

(ii) .  the  work  by  Madelung’s  team  at  the  Flugtechnisches  Institut,  Stuttgart  leading  to  the  development  of  ribbon 

parachutes  in  about  1934.  These  were  necessary  in  order  to  increase  the  stability  and  reduce  the 
inflation  loads  on  parachutes  required  for  the  in-flight  and  landing  deceleration  of  aircraft; 

(iii) .  manoeuverable  gliding  parachutes.  These  have  a  complex  history  of  development,  achieving  designs  capable 

of  contemporary  commercial  development  through  Jalbert’s  1961  ram-air  design. 

Parachute  aerodynamics  has  correspondingly  developed  from  its  initial  rile,  that  of  providing  a  service  to  aid  the 
understanding  and  extension  of  full-scale  flight  trials,  conducted  on  descending  parachutes  in  the  atmosphere.  Instead, 
it  has  become  the  means  whereby  appropriate  wind  tunnel  tests,  together  with  the  analytical  prediction  of  parachute 
performance  and  stability  characteristics,  can  supplement  without  supplanting  these  flight  trials. 

9.1  AERODYNAMIC  PROBLEMS  IN  FULL-SCALE  FLIGHT  TESTINC 

As  has  been  indicated  in  Section  6. 1 ,  the  testing  of  parachutes  during  free  flight  in  e  atmosphere  appears  to  be 
the  most  obvious  experimental  procedure.  However,  because  the  atmospheric  environment  is  uncontrolled  it  can  be  a 
difficult  medium  in  which  to  make  aerodynamic  measurements  which,  having  been  made,  can  be  as  difficult  to 
interpret.  Determination  of  the  parachute  drag  is  still  a  fundamental  problem  but  this  is  given  by  equation  234  once 
the  parachute’s  rate  of  descent  is  known.  Near  the  ground  a  mean  rate  of  descent  can  be  determined,  either  crudely  by  a 
timing  process  or  by  more  sophisticated  kindtheodolite  methods,  but  both  of  these  techniques  are  limited  to  relatively 
low  altitudes. 

Measurcinen'  of  the  manner  in  which  the  canopy’s  aerodynamic  coefficients  vary  with  the  angle  of  attack  is 
limited  by  the  difficulty  in  determining  this  angle  in  flight.  It  can  only  be  estimated  when  the  instantaneous  direction 
of  the  relative  airflow  is  known  with  certainty.  Since  the  parachute  camopy  is  bluff  the  flow  around  it  is  strongly 
influenced  both  by  its  own  shape  and  that  of  its  payload.  The  only  known  method  of  estimating  the  relative  airflow 
is  from  frame-by-frame  study  of  cind-ftlm  records  and  as  these  are  unlikely  to  include  any  means  of  flow  tracing,  at 
best  such  a  technique  can  only  be  a  very  approximate  process. 

Much  of  the  parachute's  flight  performance  and  stability  analysis  must  be  determined  by  visual  inspection, 
supported  by  photographic  records.  Because  of  this  limitation  it  is  necessary  to  judge  the  stability  or  the  instability 
in  pitch  of  a  descending  parachute  solely  in  terms  of  the  angle  through  which  the  descending  parachute  oscillates. 
During  its  descent  the  variation  of  such  art  angle  can  be  measured  through  the  use  of  gyroscopically-con trolled 
instruments.  Although  this  angle's  amplitude  is  of  significance  if  the  parachute  is  unstable  in  pitch,  for  a  stable 
parachute,  as  has  been  indicated  in  Section  23.2,  its  amplitude  may  be  as  much  a  function  of  the  local  atmospheric 
instabilities  as  it  is  of  the  degree  of  the  parachute's  stability  in  pitch.  The  widespread  practice  of  quoting  average 
angles  of  oscillation,  particularly  for  stable  parachutes,  can  therefore  be  misleading. 

Wherever  possible,  canopy  inflation  tests  are  best  carried  out  in  the  atmospheric  environment  using  full-scale 
parachutes.  Depending  on  the  parachute  application,  different  designs  of  launchers  are  used  to  cany  these  parachutes 
up  to  their  deployment  altitude.  As  the  parachutes  deploy  the  required  aerodynamic  data  is  relayed  back  to  the  ground 
control  station. 

9.2  AERODYNAMIC  PROBLEMS  IN  WIND  TUNNEL  TESTING 

Although  the  ultimate  criterion  of  a  parachute's  performance  is  its  behaviour  in  die  atmospheric  environment,  for 
a  better  appreciation  of  its  aerodynamic  characteristics  it  is  clearly  necessary  to  supplement  performance  data  gathered 
there  with  information  from  appropriate  wind  tunnel  tests.  In  order  to  perform  these  tests  the  difficulties  in  achieving 
geometric  similarity  between  the  wind  tunnel  model  and  the  full-scale  prototype  parachute  canopy,  which  have  been 
described  in  Section  63,  must  be  overcome.  These  include  blockage  constraints  around  the  bluff  parachute  canopy 
imposed  by  the  wind  tunnel  walls  and  which.  Section  6.6  indicates,  necessitate  testing  in  a  facility  whose  working 
section  cross-sectional  area  is  some  twenty  times  the  parachute  canopy’s  projected  area.  Parachutes  required  for  re¬ 
entry  vehicles  and  for  extra-terrestrial  applications  must  be  tested  in  wind  tunnels  operating  at  full-scale  Mach  number 
and,  as  has  been  outlined  in  Section  8.4,  this  is  particularly  important  for  parachutes  operating  in  the  transonic  Mach 
number  range. 

Although  procedures  by  which  canopy  inflation  loads  under  finite  mass  conditions  can  be  determined  in  wind 
tunnels  have  been  outlined  in  Section  6.43.1,  they  are  only  appropriate  for  large  scale  wind  tunnel  models,  because 
the  dimensionless  inflation  force/time  signature  relationship  described  in  Section  53.4  is  dependent  on  the  model 
shape  and  the  smaller  the  wind  tunnel  mode!  the  more  difficult  it  is  to  model  its  shape  faithfully.  Leeu  has 
confirmed  that  when  determining  inflation  loads  the  shape  of  the  canopy  includes  its  flexibility  and  that  canopy 


flexibility  has  a  marked  effect  on  the  values  obtained  for  both  the  dimensionless  opening  time  and  dimensionless 
opening  force. 

The  only  solution  to  these  wind  tunnel  problems  is  to  use  costly  test  facilities  that  will  be  adequate  for  the 
proposed  experimental  programme. 

9.3  AERODYNAMIC  PROBLEMS  IN  THE  ANALYTICAL  DETERMINATION  OF  A 
PARACHUTE’S  CHARACTERISTICS 

Following  the  work  of  Rosenhead’ *  and  of  others,  as  has  been  described  in  Section  7.3,  a  promising  start  to  the 
development  of  a  theoretical  model  capable  of  representing  the  continuous  vortex  sheet  which  is  shed  by  parachute 
canopies  has  been  made.  Srickland711  has  made  clear  that  applications  of  this  technique  to  parachutes  are,  as  yet, 
in  their  infancy.  Currently,  analytical  solutions  are  sought  for 

(i) .  the  steady-state  pressure  distribution  over  an  imporous  canopy  at  zero  angle  of  attack; 

(ii) .  the  effect  of  canopy  porosity  on  these  last  solutions; 

(iii) .  the  steady-state  pressure  distribution  as  the  angle  of  attack  is  varied,  making  the  canopy  asymmetric  to  the 

flow; 

(iv) .  the  corresponding  pressure  distributions  when  the  flow  is  unsteady; 

(v) .  the  effect  on  these  pressure  distribution  of  the  canopy  shape  varying  with  time,  as  occurs  in  the  inflation 

process;  and 

(vi) .  compressibility  effects  on  all  of  these  solutions. 

There  would  appear  to  be  excellent  prospects  of  achieving  such  solutions  in  the  not-too-distant  future. 

9.4  EFFECTS  OF  GROWTH  IN  COMPUTER  POWER 

The  single  most  significant  factor  in  the  contemporary  development  of  parachute  aerodynamics  is  the  growth  in 
computer  power,  made  readily  available  through  developments  in  both  mainframe  and  mini-computers.  The  power 
that  is  now  available  makes  possible  the  replacement  of  semi-empirical  and  strongly  experimentally-orientated 
methods  of  analysis.  Not  only  does  the  computer  readily  access  a  vast  amount  of  experimental  data  so  that  these  can 
be  compared  for  design  optimisation  studies  but  it  enables  solutions  to  be  made  of  arrays  of  non-linear  equations 
where  linearisation  or  special-case  solutions  were  all  that  used  to  be  available.  It  also  demands  a  new  type  of 
experimental  approach,  since  analysts  now  require  the  answers  to  different  kinds  of  questions.  Not  only  does  the 
computer  make  it  possible  to  achieve  complete  solutions  to  sophisticated  problems,  but  through  its  data  storage 
capabilities  and  its  ability  to  achieve  numerical  solutions  to  differential  equations,  it  can  identify  the  major 
independent  panmeters  in  any  investigation  and  determine  how  required  dependent  parameters  vary  as  a  function  of  any 
one  of  them.  As  a  direct  consequence  of  these  computer-based  parameter  identification  and  optimisation  techniques 
industrially-relevant  predictive  studies  ate  now  made. 

9.5  REWARDING  FIELDS  IN  CONTEMPORARY  EXPERIMENTAL  RESEARCH 

For  much  contemporary  wind  tunnel  research  Section  9.2  has  indicated  that  large,  high-speed  wind  tunnels, 
available  only  at  national  or  international  levels,  may  be  necessary.  However,  other  problems  can  be  identified  which 
are  capable  of  study  in  more  locally-available  experimental  facilities. 

9.5.1  The  Aerodynamics  of  High  Performance  Gliding  Parachutes 

Lingard*1  states  that  the  optimum  performance  for  conventional  ram-air  gliding  parachutes  occurs  at  an  aspect 
ratio  of  about  2:1  when  the  ratio  of  lift  to  drag  is  about  3:1.  By  developing  a  mathematical  model  for  the 
aerodynamics  of  ram-air  parachutes  he  has  identified  the  factors  which  limit  this  gliding  parachute’s  performance. 

To  obtain  a  higher  lift  to  drag  ratio  it  is  first  necessary  to  close  the  leading  edge  of  the  parachute,  leaving  a  single 
open  cell  on  the  under  surface,  litis  leads  to  a  substantial  reduction  in  drag  but,  because  the  internal  pressure  of  the 
wing  is  then  somewhat  less  than  the  stagnation  pressure,  the  leading  edge  of  the  parachute  will  now  collapse  inwards 
and,  in  order  to  maintain  the  inflated  shape,  it  is  necessary  to  sweep  the  leading  edge  backwards.  Where  this  has  been 
done,  the  resulting  swept-wing  closed-cell  concept  for  the  ram-air  gliding  parachute  develops  lift  to  drag  ratios  ratios 
which  exceed  5:1.  Ratios  of  up  to  6:1  have  been  predicted,  doubling  the  glide  ratio  from  the  maximum  currently 
possible  with  the  contemporary  unswept  ram-air  parachutes. 

In  experiments  conducted  in  the  7.3  m.  (24  A.)  diameter  low-speed  wind  tunnel  at  the  Royal  Aircraft 
Establishment,  Fam  borough,  U.K.,  lift  to  drag  ratios  in  excess  of  5:1  have  been  measured  on  such  gliding  parachutes. 
Currently,  experimental  research  into  high  performance  gliding  parachutes  is  continuing  in  order  to  obtain  this 
performance  improvement  without  impairment  of  either  the  parachute’s  inflation  characteristics  or  its  stability  and 
controllability. 

9.5.2  Tbe  Aerodynamics  of  Routing  Parachutes 

An  important  application  for  rotating  parachutes  which  was  mentioned  in  Section  23.4  is  to  provide  a  facility  by 
which  submunitions  can  rotate  about  one  of  their  body  axes  during  their  initial  deceleration  and  subsequent  steady 
descent  Ibrahimu  has  described  the  search  pattern  which  is  traced  by  certain  forms  of  submunition  during  such  a 
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descent.  Pepper*-2  shows  how  these  canopies  are  designed  to  generate  me  necessary  rotation,  a  component  of  the 
total  aerodynamic  force  providing  the  torque  which  results  in  the  canopy  autorotation.  He  also  describes  other  benefits 
which  can  follow  from  canopy  rotation,  these  are  increased  canopy  drag  and  improved  stability  characteristics  in 
pitch.  For  the  rotating  parachute  application  which  he  describes,  that  of  recovery  of  high  performance  re-entry 
vehicles,  high  rates  of  rotation  are  desirable  and  with  these  the  drag  coefficient  developed  on  the  autorotating  parachute 
described  was  some  four  or  five  times  the  value  that  it  would  have  had  if  this  parachute  did  not  rotate. 

DohetT,  MQnscher  and  Saliaris*J  defined  a  rotor  quality  number  R<,  in  terms  of  the  canopy  drag  coefficient  Q, 
and  the  rotor  coefficient  Q. 


Rq  = 

cfVcD 

(9.1) 

where 

P 

it 

f  D./Vd 

(9.2) 

and  f  is  the  autorotation  frequency  in  revolutions  per  second,  D„  is  the  canopy  nominal  diameter  and  VD  its  rate 
of  descent,  so  that  C,  =  1/J, I  being  the  advance  ratio,  as  normally  defined  for  propellors.  The  significance  of  this 
rotor  quality  number  has  been  discussed  by  Doherr  and  Synofak*'*,  who  describe  a  series  oi  wind  tunnel  experiments 
devised  to  measure  for  a  rotating  guide  surface  parachute.  The  methods  of  performance  evaluation  which  they 
describe  are  recommended  for  more  general  application  in  wind  tunnel  tests  on  rotating  parachutes. 

9.5.3  Experiments  to  Further  the  Application  of  Vortex  Sheet  Theories  to  Parachutes 

Writing  of  the  development  of  vortex  sheet  theories,  a  quotation  included  in  the  introduction  to  Section  7  stated: 
“It  is  at  this  point  that  we  are  handicapped  by  the  fact  that  experimental  techniques  are,  at  this  moment,  lagging 
behind  the  advance  of  theory”.  Given  the  desire  and  the  facilities  to  perform  these  necessary  experiments,  what  facts 
need  to  be  discerned  from  then? 

Much  of  the  experimental  work  performed  on  parachutes  in  wind  tunnels  has  been  with  a  view  to  determining  the 
mean  values  of  the  aerodynamic  coefficients  which  are  developed.  Any  observation  that  these  coefficients  might 
fluctuate  in  magnitude  has  been  considered  to  arise  from  extraneous  factors,  such  as  poor  wind  tunnel  design  or  the 
blockage  constraint  imposed  by  the  model. 

A  deeper  understanding  is  now  required  of  the  nature  of  the  flow  around  bluff  bodies  in  general  and  around 
parachute  canopies  in  particular.  A  fundamental  question  is:  when  the  canopy  is  set  symmetrically  in  the  flow,  what 
is  the  frequency  at  which  it  sheds  vortices  and  what  determines  this  frequency?  Cain  it  be  disassociated  from  the 
stiffness  of  the  sting  support?  Is  it  a  function  of  the  properties  of  the  canopy  fabric?  Is  the  oscillation  in  resonance 
with  a  much  lower  amplitude  driving  oscillation?  Is  it  eynolds  number  dependent? 

If  this  flow  in  the  wake  is  periodic,  then  are  the  aerodynamic  forces  developed  on  the  canopy  correspondingly 
periodic  and  if  they  are,  with  what  amplitude  do  they  vary? 

Since  the  purpose  for  which  this  information  is  required  is  the  construction  of  a  vortex  sheet  model,  what  is  the 
simplest  form  in  which  these  new  physical  insights  can  be  expressed? 
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10.  POSTSCRIPT 

There  continue  to  he  a  number  of  fundamental  and  challenging  problems  in  parachute  aerodynamics,  the  reasons 
for  which  are  very  similar  to  what  they  were  when  W.D.Brown  published  ' Parachutes' 14  in  1951,  declaring  as  his 
aim  that  of  selecting  ‘the  principal  aerodynamic  characteristics  of  parachutes  and  the  various  known  factors  which 
affect  these  characteristics’.  The  bluff  body  of  the  conventional  parachute  canopy  still  possesses*  a  non-rigid  structure 
and  has  a  mass  which  is  of  the  order  of  mass  of  the  air  which  it  displaces.  However,  Brown  could  not  have 
appreciated  how  the  passage  of  35  years  would  bring  such  a  diversity  of  parachute  application  or,  accompanying  this 
diversity,  the  necessity  for  a  deeper  fundamental  understanding  of  physical  principles.  One  of  the  many  benefit  of 
trying  to  bring  threads  together  in  this  AGARDograph  has  been  a  much  deeper  appreciation  of  the  developments 
which  have  occurred  in  such  a  short  period  of  time. 


Twenty  years  later,  in  a  review  on  aerodynamic  decelerators  written  from  the  Sandia  Laboratories  at  Albuquerque, 
Pepper  and  Maydew1-1  remarked  that  Brown  had  written  the  only  bode  on  parachute  technology  that  was  known  to 
them.  At  that  time  extensive  ribbon  parachute  development  work  was  taking  place  at  die  Sandia  Laboratories. 
Although  this  AGARDograph  might  not  have  done  justice  to  that  particular  activity,  it  does  contain  more  than  a 
dozen  references  to  significant  aerodynamic  research  performed  subsequendy  by  individuals  working  at  the  Sandia 
Laboratories. 

What  led  to  these  developments  has  been  both  the  circumstances  and  the  individuals  whose  contributions  to  the 
subject  have  been  demanded  by  these  circumstances.  Knacke141  states  that  'it  was  World  War  II,  its  forebearings  and 
its  aftermath,  that  started  the  widespread  applicadon  of  parachutes  for  the  air  drop  of  troops  and  supplies,  for  the 
retardation  of  ordnance,  the  in-flight  atd  landing  deceleration  of  aircraft  and  the  recovery  of  missiles,  drones  and 
spacecraft’  and  although  Brown  wrote  with  direct  experience  of  that  war  he  could  not  have  foretold  all  of  this 
aftermath. 


Brown  mentions  Heinrich  and  Knacke,  the  first  of  whom  he  describes  as  a  German  technician,  who  designed 
during  World  War  II  the  ‘mushroom’  or  ‘beret’  parachute,  primarily  for  dropping  heavy  bombs  and  sea  mines.  He 
referred  to  Knacke  as  a  German  scientist  who  appears  to  have  invented  the  "TaschengurT.  In  this  AGARDograph 
there  are  nine  separate  references  to  Heinrich’s  work  and  the  contribution  made  to  parachute  research  and  development 
by  the  Minneapolis  postgraduate  school  which  he  established,  has  been  outstanding.  Similarly,  to  rough  his  research 
and  leaching.  Knacke  has  made  a  series  of  memorable  contributions  to  this  subject  Both  have  been  honoured  by  the 
AI AA  for  the  outstanding  pans  they  have  played  in  the  development  of  parachutes. 

Over  these  years  the  most  significant  contribution  to  the  dissemination  of  information  on  parachute  technology 
has  been  the  AIAA  Aerodynamic  Decelerator  Conferences,  held  in  the  United  States  every  two  and  a  half  years.  Most 
active  participants  in  parachute  aerodynamics  throughout  the  world  have  been  present  and  have  contributed  to  these 
Conferences.  Pepper  and  Maydew  refer  to  the  first  and  second  ones,  held  in  1966  and  1969.  Since  that  time  these 
conferences  have  grown  in  significance.  In  October  1986  the  9th  AIAA  Aerodynamic  Decelerator  and  Balloon 
Technology  Conference  took  place  in  Albuquerque,  New  Mexico. 
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14.  Abstract 


This  AG  ARDograph  discusses  the  principal  aerodynamic  characteristics  of  parachutes  and 
the  factors  which  affect  those  characteristics.  It  takes  into  account  many  of  the  publications 
that  were  summarised  in  the  1 963  and  1 97JJUnitcd  States  Air  Force  Parachute  Design 
Guides,  the  proceedings  of  biennial  American  Institute  of  Aeronautics  and  Astronautics 
Aerodynamic  Dccclerntor  Conferences,  the  Helmut  G.Heinrich  Dccclcrator  Systems 
Engineering  Short  Courses  held  in  1 983  and  1 985  and  *  fuc  Parachute  Recovery  System 
Design  Manual’,  snortly  to  be  published  by  the  United  States  Naval  Weapon  Center. 

It  is  anticipated  that  its  main  readers  will  be  recent  engineering  graduates  entering  research 
establishments,  parachute  companies  or  related  industries  so  some  appreciation  of  tasic 
mechanics,  the  principles  of  computing  and  elementary  fluid  mechanics  on  the  part  of  the 
reader  has  been  assumed. 

4  Its  content  includes  Steady-State  and  Unsteady  Aerodynamics,  Parachute  Deployment  and 
Inflation,  Experimental  Investigations,  Methods  of  Analysis,  Extra-Terrestrial  Parachute 
Applications  and  Some  Suggestions  for  Future  Research. 

This  AGARDograph  has  been  produced  at  the  request  of  the  Fluid  Dynamics  Panel  of 
AGARD. 


